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PREFACE 


This NASA Conference Publication contains the proceedings of the NASA 
Conference on A'lvanced Technology Airfoil Research held at Langley Research 
Center on March 7-9, 1978, which have unlimited distribution. Conference 
oochairmen were Alfred Sessow, NASA Headquarters, and Robert P. Bower, Langley 
Research Center. Hono;ary cochairmen were Ira H. Abbott, NASA Headquarters 
(retired) , and Richard T. Nhitoomb, Langley Research Center . 

The conference i«s planned to provide a con^rehensive review of all NASA 
airfoil research, conducted both in-house and under grant and contract. In 
addition, a broad spectrum of airfoil research outside of NASA was reviewed. 

A total of 64 technical papers were presented at 12 sessions. Six workshops 
were also held to discuss progress, further immediate and long-range research 
needs, and important unresolved issues. A roundtable discussion summarised 
the technical sessions and workshops. 

This volume contains papers presented at technical sessions covering the 
following subjeotst 

(1) Airfoil Analysis and Design of Single-Element Airfoils 

(2) Airfoil Analysis and Design of Multielement Airfoils 

(3) Airfoil Analysis and Design Topics 

(4) Research Facilities and Test Techniques 

(5) Facilities and Test Technique Topics 

(6) Unsteady Aerodynamics 

The major thrusts of the technical sessions were in three areast devel- 
opment of computational aerodynamic codes for airfoil analysis and design, 
development of experimental facilities and test techniques, and all types of 
airfoil applications. The conference proceedings are presented in two volumes* 
Volume I is unclassified with unlimited distribution and Volume II is unclassi- 
fied but with limited distribution. 

The included pap>ers are largely as submitted as camera-ready copy. Only 
minor editorial revisions have been made and a title page and abstract have 
been added. 

Use of trade names or names of manufacturers in this report does not 
constitute an official endorsement of such products or manufacturers, either 
expressed or implied, by the Netional Aeronautics and Space Administration. 


P. K. Pierpont, Conference Organiser 
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A PARAMETRIC EXPERIMENTAL STUDY OF THE 30 

SLOTTED-WALL BOUNDARY CONDITION 

Joel L. Everhart* and Richard W. Barnwell 
NASA Langley Research Center 

SUMMARY 

An experimental study of slotted upper and lower walls In a two-dimensional 
transonic wind tunnel with solid sidewalls Is reported. Results are presented 
for several slot spacings and slot openness ratios# The experlinental data are 
pressure measurements which were made on an airfoil model and on a sidewall near 
one of the slotted walls. The slotted-wall boundary-condition coefficient, 
which relates the pressure and streamline curvature near the wall. Is determined 
from the wall pressure measurements. The measured wall -Induced Interference Is 
correlated with the experimental values for the boundary-condition coefficient. 
This correlation Is compared with theory. 


INTRODUCTION 

Slotted walls have been used to relieve blockage effects In transonic wind 
tunnels for three decades. In this paper an experimental study of slotted 
walls with different slot spacings and openness ratios in a two-dimensional 
transonic tunnel with solid sidewalls Is reported. Pressure measurements were 
made both on a model In the tunnel and near one of the slotted walls so that 
the wall -Induced Interference on the model and the nature of the flow near the 
slotted wall could be determined and correlated. These results can be used In 
the design and estimation of Interference of other wind tunnels with slotted 
walls. 

There was a practical Incentive for the present test. The amount of 
experimental data on the flow near slotted walls Is scarce (only three experl- 
mental data points (refs. 1, 2, and 3) are known to have been published prior 
to this test), and the data that have been published do not agree with either 
of the basic theoretical nmdels (refs. 4 anti 5). It should be noted the 
theoretical model presented In reference 3 Is, In effect, an empirical modifl- 
cation to the method of reference 4# It should also be noted that the theoretv 
cal model of reference 5 Is simply a correction to that of reference 1. An 
analysis of the theoretical and experimental results (ref. 6) has shown that 
the design of two-dimensional tunnels with little or no blockage and streamline 
curvature effects Is feasible If the experimental results are correct but 
unfeasible If the generally-accepted theoretical model of Davis and Moore 
(ref. 4) Is correct. It was concluded that a parametric study Involving wall- 
geometry parameters, Mach number and angle of attack should be made. The 
Reynolds number could not be varied Independently since an atmospheric facility 
was used. 


♦Former Graduate Research Scholar Assistant with Joint Institute for 
Advancement of Flight Sciences, George Washington University. 
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There are three major steps* in the estimation of interference effects 
caused by slotted walls. First, the general form of the boundary condition 
for flow near the wall must be determined. Second, the unknown coefficients 
in this boundary condition must be evaluated in terms of the slotted wall 
geometry. Third, the relationship between the wall-induced interference and 
the values of the boundary-condition coefficients must be determined. In this 
paper • results which affect each of these three major steps will be presented. 
In the case of the first step» the Ideal s1otted*wa11 boundary condition Is 
used. For the second and third steps, correlations of parametric experimental 
data for the tunnel -geometry to boundary-condition-coefficient relationship 
and the boundary-condition-coefficient to interference-effects relationships 
are presented. 


SYMBOLS 


A empirical boundary-condition coefficient 

a slot spacing 

B porous wal’ boundary-condition coefficient 

Cj^ normal -force coefficient 

Cm slope of the normal force curve 

o 

Cjn pitching-moment coefficient 

Cp pressure coefficient 

Cp* critical pressure coefficient 

S»PLENUM P''*®ssure coefficient in the plenum chamber 
c airfoil chord 

h tunnel semi -height 

K slotted-wall boundary-condition coefficient 

k slotted-wall correlation parameter, k » |k 

^PLENUM number based on plenum pressure 

Ppi^ENUM Plenum pressure 

R Reynolds number 

X distance along tunnel center line 
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upstream location used In evaluation of equation (7) 
y distance nnrrual to tunnel center line 

a angle of attack 

6 slot width 

B flow angle 

0Q value of 0 at x^ 

SLOTTEDoUALL BOUNDARY CONDITION 

The terms used In the slotted-wall boundary condition are Illustrated In 
figure 1* The usual form of the boundary condition Is derived from the Ideal 
slot condition* which states that the pressure In the slots Is equal co the 
plenum pressure. The resulting boundary condition Is 


a2 



S. plenum 


( 1 ) 


where Cp and Cp^pj^gjjym are the pressure coefficients In the tunnel near 

the wall and In the plenum* respectively. The quantities 6* a* and e are 
the slot width* the slot spacing* and the flow angle In the tunnel near the 
wall. Note that the quantity da/6 Is the nondlmenslonal cross-flow velocity 
at the slot. The coordinate x Is the distance along the tunnel axis* and 
K Is the slotted-wall boundary-condition coefficient. 


When the cross flow at the slot Is small* equation (1) can be approximated 
as 


C 


P 



S.plenum 


( 2 ) 


DETERMINATION OF BOUNDARY-CONDITION COEFFICIENT 

The experimental procedure which was used to measure the boundary-condition 
coefficient Is described In this section. Since the flow In the tunnel near 
the slotted walls Is Irrotatlonal and the disturbances are small* the small - 
disturbance Irrotatlonal condition 


dd 

n 


1 

7dy 


(3) 
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can be used to rewrite equation (1) as 


J df. 

Cp -i-4-t K- 




PLENUM 


( 4 ) 


*. ** y measured perpendicular to the slotted top and 

bottoin walls of the tunnel. Finite-difference approx linatlons for tne pressure 

JSfcelSl!''L 4 . 5 P 2"J derivative aCn/ay are made and are evaluated with 
K! !f !5 obtained from pairs of orifices drilled In the solid sidewall near 
the slotted top wall. The location of these orifices Is Indicated In floure 2 
The expressions for these quantities are ^ 


and 




ay 


^p(x»y^) ** Cp(x,y2) 


yi -yg 


(5) 


( 6 ) 


where y^ and yg are values of the y coordinate for two orifice rows. The 
flow angle e Is obtained from an Integration of the Irrotatlonal condition as 

X 


e = 


-1 


/• 

/sy^*** 


+ a. 


(7) 


AA coefficient K Is obtained from a comparison of 

distributions of Cp and aCp/ay. 


EXPERIMENTAL FACILITY AND TEST CONDITIONS 

The experimental facility Is the Langley 6- by 19-1nch transonic tunnel. 

This Is a two-dimensional blowdown tunnel In which transonic airfoils are tested. 
The model used for the present study was an NACA 0012 airfoil with a chord of 
15.24 cm (6 In.). Pressures on the airfoil surface were measured with strain- 
gage pressure transducers with a range of ±10.34 N/cm^ (±15 psi). The accuracy 
of these gages Is 1 percent of full scale. 

Three rows of pressure orifices have been drilled In one of the solid 
sidewalls at distances of 2.54 cm (1 In.), 3.81 cm (1.5 In.), and 5.08 cm 
(2 In.) from the slotted upper wall. Two of the orifice rows extend upstream 
beyond the point where the slots start st that the Integration to determine 
the flow angle e can start at a point where 0 Is known. The sidewall 
pressures were measured with variable capacitance precision transducers which 


have reading errors no larger than .5 percent of reading. The transducers 
were referenced to the plenum pressure In order to keep the reading small and 
hence minimize measurement error. Near the model, transducers with a range 
of +5.17 N^m" (+7.5 psi) were used. At other orif ices where the measurements 
wer¥ less critical, transducers with a range of +10.34 N/cm* (£15 psi) were 
used. In all, there are 52 orifices on the sidewall where pres«5Mre measi'^e- 
ments were made. 

It should be noted that the pressure measurements on the tunnel sidewall 
and on the model were all made at the same time. In references 1, 2, and 3, the 
flow angle was measured one po.int at a time. This noncoincident data was then 
differentiated numerically to determine the derivative d#/dx. 

The test conditions are shown In figure 3. Tests were conducted for 
slotted walls with openness ratios of .01, .025, .05, .075, .10, and .15 and 
slot spacings of 3.81 cm (1.5 In.), 7.62 cm (3 In.), and 15.24 cm (6 In.). 

These slot spacings correspond to 4, 2, and 1 slots, respectively. In the 
6>1nch-w1de slotted wall. Some data were also obtained In a closed tunnel. 

It should be noted that some combinations of the openness ratio and slot 
spacing values listed above were not tested. It should also be noted 
some of the combinations were tested In two different conflgurat” .r 
example, some of the single-slot combinations were tested v<; . . slot 

In the middle and with 1/2 slot at each sidewall. The bas*. »hI 1 configurations 
are shown In figure 3. 

The wall configurations were tested at Mach numbers of .5, .6, .7, .8, 

.825, .85, .875, and .9. Each wall was tested with the airfoil at several 
different angles of attack. Since the airfoil model Is symmetric, results 
for the upper and lower slotted walls were obtained with the present experi- 
mental facility, which Is only Instrumented near the upper slotted wall, by 
placing the airfoil at positive and negative angles of attack. 


EVALUATION OF BOUNDARY-CONDITION COEFFICIENT 

An example of the evaluation of the coefficient K Is depicted In 
figure 4. It Is assumed that the boundary condition !s the combined slotted 
wall, porous-wall condition with an added constant term. This boundary 
condition Is written as 


a2 ^2 

S — W ' S.plenum 


dC 

-Ka + B8 + A 


( 8 ) 


The constant A Is added to account for the fact that the plenum pressure 
Is slightly lower than the average pressure In the tunnel even when no model 
Is present. 

Results for the distributions of Cp, 3Cp/ay, and 8 are plotted 
against tunnel station. It can be seen that the distribution of Cp Is 
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similar to that of dCp/fly but unlike that of Q» For example* the local 

maximum and minimum on the curve for $ do not occur near similar features 
on the curve for Cp, but the local maxlmums on the curves for Cp and 

dCp/8y occur at the same place. It Is concluded that, for all means and 

purposes, the coefficient B Is zero, A good match of the curves for C 

and 3Cp/dy Is obtained when the coefficients K and A have the values 

1.7 and ,022, respectively. 

Results for the value of the slotted-wall coefficient K for several 
openness ratios and slot spacings are depicted In figure 5. All of the present 
data shown are for a Mach number based on plenum pressure of ,7. Also shown 
on the figure are the three previous experimental data points (refs. 1, 2, and 
3) and the results of two theoretical models (refs. 4 and 5), It can be seen 
that the four present data points for the tir.*o*slot configuration vary con- 
sistently with openness ratio. It can also be seen that the coefficient K 
appears to be dependent on the slot-spacing to tunnel-height ratio. 


DETERMINATION OF INTERFERENCE EFFECTS 
Blockage Interference 

Blockage Interference Is a result of either the underexpansion or the 
overexpansion of the fluid about the model as compared to free-air flow. This 
effect Is directly related to the wall characteristics (I.e., openness ratio, 
number of slots, slot arrangement, etc.). In figure 6 the effect of the open- 
ness ratio on the shock location Is shown at non-lifting conditions for a wall 
with 1 slot. The Mach number Is about .833. The openness ratio has been 
varied from .01 to .05. On this airfoil at this Mach number, a 30-percent 
chord change In shock location Is Indicated with a change of 4 percent In 
openness ratio. 

The effect of the number of slots and their arrangement In the wind 
tunnel Is a constraint which. In general, has been overlooked. In figure 7 
the effect of slot arrangement on shock location for an openness ratio of .05 
and Mach number of about .83 Is shown. The open circles represent a wall with 
one slot In the middle of the tunnel. The closed circles represent a wall 
with one-half slot located at each sidewall. Open squares are data obtelned 
from a two-slot wall with each slot located a distance of one-half the slot 
spacing from the nearest wall, while the closed squares represent a two-slot 
wall with one slot In the center and one-half slot on each sidewall. The 
differsnee in the shock wave location obtained with the two one-slot walls and 
the difference obtained with the two two-slot walls Is about 2 percent In each 
case. Differences of this magnitude can be expected from either Instrumenta- 
tion error or the Inability to repeat test conditions. These data Indicate 
that shock location Is relatively Insensitive to slot arrangement. However, 
there Is a larger difference In shock wave location between the one-slot and 
two-slot results. This difference Is between 5 and 10 percent for the case 
shown. 
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Lift Interference 
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Streamline Curvature Interference 
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CONCLUDING REMARKS 
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Figure 1.- Slotted-wall boundary condition 
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TRANSONIC ASSESSMENT OF TWO-DIMENSIONAL WIND 

TUNNEL WALL INTERFERENCE USING MEASURED WALL PRESSURES 

William P Kemp, Jr. 

NASA Langley Research Center 


SUMMARY 


i Is described for assessing wall Interference In two-dimensional 

P««Hsm and accuracy than Is achieved with methods 
Interference theory. Measured pressure distributions 
* nonlinear transonic form of the 

fJJfJrJ' accurate computational reproduction of the 

velocity perturbations are then extracted 
both corrections to the tunnel test conditions and a measure of the 
adequacy of these corrections to account for tunnel Interference. Appll- 

lllustr t d"* ® transonic tunnel with variable porosity walls Is 


INTRODUCTION 


The methods used most frequently for predicting wall Interference In 
tunnels draw their basic principles from the classical wall 

ir?!; ® theoretical representation of the tunnel flow Is 

II?®!'' superposition of two-dimensional potential fields due to 
jTOdel lift, thickness, drag, and the wall reaction to each, and the wall 

®*'® ®®t^sf1ed uniformly along surfaces extending from 
upstream to downstream Infinity. This approach has been satisfactory In 
general when applied to solid wall wind tunnels at purely subsonic speeds. 

nAA^Al!Ii^. 4 ^^®-i?^'^®[!^ transonic test sections having either slotted or 
perforated walls, boundary condition expressions were developed In linear 
form to repi^sent these walls In the classically based procedures. Subsequent 
exp®r1ence has shown that the resulting wall Interference predictions are^ 
often Inadequate for correlating airfoil test results from different facilities. 
4 ® P®P®<^» several areas from which Inaccuracies arise will be 

loentlfled and a new procedure will be described which should lead to signifi- 
cantly improved accuracy for wall Interference assessment. 
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SYMBOLS 


Cp pressure coefficient 

c airfoil chord 

c„ airfoil section normal-force coefficient 

h tunnel height 

M Mach number 

u* V total perturbation velocity components In x and y d1 recti ons, 
respectively; normalized by tunnel reference velocity 

Um* wall -Induced perturbation velocity components In x and y directions, 
respectively; normalized by tunnel reference velocity 

Vp velocity component normal to wall 

X. y coordinates In streamwise and crossflow directions, respectively 
a angle of attack 

T open area ratio of perforated wall 

SOURCES OF INACCURACY 


Several problem areas which contribute to Inaccuracy In wall Interference 
prediction, especially In transonic test sections, are pointed out In 
figure 1. Although the primary reason for using a transonic test section 
In an airfoil test facility Is to minimize the exaggeration of wall Inter- 
ference by the supercritical flows occurring at high subsonic speeds, the 
limitations of the linearized potential flow representation used In the 
classical wall Interference theory are violated In this speed range. 

The I mposi tion of wall boundary conditions uniformly on Infinitely long 
surfaces in the classical approach overlooks such real tunnel geometry 
features as the finite length of a ventilated wall, the possible existence of 
diffuser entry flaps, a choked diffuser entrance, or the presence of a wake 
survey rake and Its support. An additional source of error Is the boundary 
layer on the tunnel side walls. The pressure field around the test airfoil 
causes variations In the displacement thickness of this wall boundary layer, 
thereby violating the two-dimensional flow assumption of the classical theory. 
The resulting effects Include not only localized three-dimensional flow 
distortions but also the one-dimensional distortion of the effective tunnel 
stream tube area. 
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FtATURES OF IMPROVED METHOD 

Kis » ffi 

Li Kd 

procedure described In det^l in ® ^^J®«>*|zed potential flow 

In referencll s co«a?ed J floSJ^rSJth i 5^ the results alven 

approach as Implemlnted In referl^e 4 Ihe'^I^Irf^JLfT.*'* ® 
obtained at the University of Southamnfnn analyzed were 

course of a self-streami Inina wafl^sSdS the 

Center. The test Mach number was low^anminh^ift^l sponsored by Langley Research 

?SoW 

large wall Interference effects. * expected to produce unusually 


centeM1nrS^the*’h?I!rJlIn^^^‘'I® ^ distributions along the tunnel 
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by the method using measured boundary conditions Is probably caused by the 

To extend this cdoablllty to supercritical flows, the measured boundary 
conditions «re iSlJosed on « prtentlnl flow Je'llM 

r=;5orn?SiiorofT??:iU'^i«e‘^^^^^ 

Slln‘’iri!!Ie"ecMSS5urreSiri«lSdinipCT^?^ 

£rf3 “prn??;^ «?o ; oJx rsi'4r,s 

J?e direlopid « polyLil.ls co.spat1ble vrtth all ’"‘ij?."® 

prescrited flow direitlons at the two upstream 
of airfoil lift and drag coefficients are used to constrain the 
?lrcul.tlin and the trailing edge thickness by * /KSl*’ 

scheme at the mesh interval bracketing the airfoil leading edge. 

The wind tunnel flow calculated by this method combines the realism 
1ntroJS«d bj tre*per^!lental boundarj conditions w'‘b ‘be accurate treatm^^ 
A*f ennov*rr*^f iral flows mdde possible bv the nonlinaer transonic potential 
eJuiro""Uralsr“«rvrr/?hl wal1-1nduced pe^ 
identified part of the nonlinear tunnel flow solution, 

“nlrSr eJ^esiStlS^ ^ 

the model -induced perturbation. 

Tft facilitate both the boundary condition transfer and the subsequent 

rxterd'Ui'nrt°hr5;™rn-anJ'’?aX?i%S^^^^^^^ 

:"ccru^*%rthrSe?Sdli «u«e iSSn^^^ co"r45p"ding to*?he prescribed 
drag coefficient. 
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Because at transonic speeds the perturbation solution becomes strongly 
coupled to the Mach number of the uniform flow that Is used as a perturbation 
base, the Mach number of this uniform base flow (that Is, the far field Mach 
number) is corrected In the free air computation to account for wall Interfer- 
ence. The streaiiiWlse component of the wall-induced velocity Is calculated at a 
prescribed chordwise location on the model after each Iteration of the free 
air flow solution and Is used to update the far field Mach number for the next 
Iteration. After convergence, the local Mach number at this location In the 
free air flow matches that In the tunnel flow because the wall-induced 
contribution simply has been shifted from the perturbation In the tunnel 
case to the uniform base flow In the free air case. This location, therefore, 
will be called the match point. Because of the velocity Jump boundary 
condition, this match occurs on both the upper and lower model surfaces. It 
follows that 1f the wall-induced blockage velocity were uniform over the model 
chord, the entire Mach number distribution (or static pressure distribution) 
over the model In the free air flow would match that In the tunnel flow. 

It Is not necessary to Iteratively correct the far field flow direction 
to account for the wall -Induced upwash because the transonic small disturbance 
equation used In the present method treats such crossflow perturbations 
linearly. Instead, one might Imagine rotating the entire free-air flow 
solution to align the local flow direction at the match point with that In 
the tunnel flow. Then If the wall Induced upwash velocity were uniform over 
the model chord, the entire model shape would match that In the tunnel flow 
and the far field velocity vector would define the free-air Mach number and 
angle of attack to which the data measured In the tunnel should be applied. 

In general, however, the wall-induced blockage and upwash velocities will 
not be found uniform over the model chord and the amount of nonuniformity Is 
Indicative of the residual error In the data after correcting for mean values 
of the wall-induced Increments In Mach number and angle of attack. Although 
at low speeds dav.a corrections sometimes are applied as linearly super- 
imposed first order effects of gradients in wall-induced velocity, these 
corrections overlook such nonlinear phenomena as shock movement and are, 
therefore, of questionable validity at supercritical speeds. 


ILLUSTRATION OF TRANSONIC ASSESSMENT 


Data obtained during tests of a 10-percent thick supercritical airfoil 
In the Lockheed Compressible Flow Facility and discussed in reference 7 have 
been assessed for wall Interference using the present method. This facility 
has perforated top and bottom walls with variable porosity control and was 
equipped with pressure measurement rails to determine the static pressure 
distributions near the top and bottom walls. Because the pressure rails did 
not extend upstream Into the solid wall nozzle, the upstream boundary of the 
tunnel flow computation region was located within the perforated region of 
the real tunnel test section and consequently the flow directions at the 
two upstream corners were not accurately known. To analyze the test results 
from this facility, therefore, the flow angle at these points was input as 
zero and the calculated tunnel flow was then rotated to achieve a best 


alignment of the model surface streamlines In the calculated flow with the 
known airfoil shape and angle of attack setting. The flow direction was 
thereby resolved within about 0.02°. Fxperlence has shown that the 
calculated airfoil surface streamline shape Is sufficiently reliable for this 
approach If the angle of attack Is small and the airfoil pressure data are 
of good quality. 

Airfoil pressure distributions measured at a Mach number of 0.8 with 
two different wall porosities are compared In figure 3. In each case the 
angle of attack In the experiment was set to obtain a normal force coefficient 
equal to the airfoil design value. The angle of attack settings required with 
the two porosities differed by more than 1°, and, even at the same Cn, 
the large differences observed between the two pressure distributions Indicate 
changes In wall Interference which cannot be accounted for simply by adjusting 
angle of attack. 

Wall Interference results obtained by applying the present method to the 
0.013 wall porosity case are shown In figure 4 for two prescribed locations 
of the match point used In updating the far field Mach number. Again, the 
results are given as distributions along the tunnel center line of the 
blockage and upwash components of the wall*1nduced velocity normalized by 
the tunnel reference velocity. Observe that although the blockage velocities 
used for the far field update differed by only about 0.005 between the two 
match point locations, the resulting wall -Induced velocity distributions are 
very sensitive to this change, particularly In the vicinity of the shock 
terminating the upper surface supercritical region. Observe also that the 
discontinuity in wall-induced velocity across the shock Is essentially 
eliminated by choosing a match point location just ahead of the shock. This 
Is because the shock strength Is equalized In the two flows by matching the 
local Mach number entering the shock. 

One might be disturbed at this point by the apparent dependence of the 
wall Interference results on an arbitrarily chosen match point location. 

Note, however, that this dependence arises from the gradient In wall-induced 
velocity and would vanish If the wall properties were such as to produce a 
uniform wall-induced velocity over the model. In such a case, the free-air 
condition to which the tunnel data properly should be applied would be 
Indicated by the present method regardless of match point location. In the 
more general case, the nonuniformities In wall-induced velocity yield an 
uncorrectable residual error and the arbitrary match point location provides 
a means of minimizing this residual error by seeking a most nearly uniform 
wall -Induced perturbation, thereby Identifying that free-air condition to 
which the tunnel data are most nearly applicable. The experience gained to 
date In applying the method to supercritical flow cases has Indicated that 
locating the match point just ahead of the shock Is a generally effective 
way to minimize the nonuniformities. 

The wall-induced velocity distributions In the vicinity of the model 
with four different wall porosities are compared In figure 5. For each 
case, representative single values of wall-induced blockage and downwash 
would be used to correct the Mach number and angle of attack. Best accuracy 
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of the resulting data requires minimizing not the magnitude of the wall-induced 
velocities but rather their nonuniformity over the model. The nonuniformities 
in the results shown in figure 5 include localized irregularities in the super- 
critical flow region which are most pronounced for the minimum porosity wall. 
With increasing porosity, the localized irregularities diminish and essentially 
vanish at 0.06 porosity, only to reappear with a somewhat altered nature at 
0.08 porosity. The localized Irregularities are superimposed on a general 
gradient which, in the case of the blockage component, is also largest for the 
minimum porosity and nearly vanishes at 0.06 porosity. The general gradient 
of the wall-induced downwash, however, crosses zero between the 0.013 and 0.04 
porosity cases and is clearly negative at 0.06 porosity. The results shown in 
figure 5 indicate, therefore, that none of the wall porosities produce data 
which are completely correctable for wall interference. The residual error is 
probably minimized with a porosity of slightly less than 0.06. 


It is tempting to speculate that the previously noted localized 
irregularities of wall-induced velocity in the supercritical region might 
arise from the same wall properties that govern the wave reflection charac- 
teristics of the wall at higher Mach nundi>ers. To date no cases have been 
analyzed at sufficiently high Mach nun^er to either support or refute this 
speculation. 


In order to compare the results shown in figure 5 with predictions from 
the classical approach, values are needed of the wall porosity parameter 
which appears as a coefficient in the perforated wall boundary condition. 
Inasmuch as the tunnel flow computation in the present method is an accurate 
reproduction of the actual tunnel flow, it can be used to provide a wall 
crossflow calibration from the actual test data. The procedure is illustrated 
in figure 6 for wall porosities of 0.013 and 0.08 although it was performed 
for all porosities. The variation of crossflow velocity v (positive inward) 
with longitudinal perturbation velocity u along both the top and bottom 
walls is plotted for each porosity. The slope of this variation is taken 
as the wall porosity parameter. In determining the slope, the portions of 
the walls nearest the model, which generally Included tne largest magnitude 
of velocity perturbations, were considered most significant. The curly 
variations near the origin in figure 6 occurred far upstream of the model 
and were ignored. 

The porosity parameters thus determined were used with the charts of 
reference 4 to predict the magnitude and gradient of the wall-induced 
»elocities for each wall porosity. These predictions are shown in figure 7 
as dashed lines and are compared in this figure with the results of the 
present method reproduced from figure 5. Considering first the wall-induced 
downwash results, the gradients predicted friMi reference 4 agree well with 
the general gradients obtained by ignoring the localized irregularities 
near the shock in the results of the present method. The downwash magnitude 
agrees reasonably well between the two methods for the 0.04 porosity case 
but significant discrepancies are apparent for the higher porosities. 
Additional information on the wall-induced downwash magnitude can be obtained 
from the angle of attack settings during the experiment. It is reasoned 
that because the airfoil normal force coefficient was the same for all four 
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wall porosfty cases* the change in angle of attack setting between any two 
porosities should reflect with fair accuracy the change In wall-induced 
downwash. Accordingly, the negative of the actual angle of attack settings 
(In radians) Is Indicated In figure 7 by the t1^.k marks. Although the agree- 
ment In absolute level with the results of the present method Is only 
coincidental, the Increments due to changes In wall porosity corroborate 
the results of the present method rather than those from the classical 
predictions. 

The magnitudes of the wall -Induced blockage velocity that were predicted 
using reference 4 are shown In figure 7 to be smaller In general and much less 
sensitive to changes In wall porosity than those assessed by the present 
methcd. The large sensitivity of the wall-induced blockage to wall porosity 
changes Indicated by the results of the present method at low porosities Is 
corroborated by the discussion In reference 7 of the same experimental data. 

In reference 7 It was shown that the airfoil pressure distribution measured 
with 0.02 porosity walls at a reference Mach number of 0.796 but reduced as 
though the Mach number were 0.817 agreed very well with that measured with 
0.04 porosity walls at a reference Mach number of 0.817. This Implies that 
the normalized blockage velocity u^ was 0.023 more positive for a wall 
porosity of 0.02 than for a porosity of 0.04; this agrees much better with 
the trend of the present method on figure 7 than with that of the predictions 
using reference 4. 

Although the reasons for discrepancy In the particular cases shown are 
not known. It should be pointed out that the present method Is responsive 
to several phenomena which are overlooked by the classical approach. These 
Include the effect of discrepancies In the empty tunnel calibration, the 
Interaction between the model and the tunnel reference pressure source, and 
the change In effective wall boundary layer thickness due to model -Induced 
flow through the wall perforations as well as the sources of Inaccuracy 
discussed In a previous section. 

In addition, an anomaly In the method of reference 4 with respect to the 
wake blockage prediction should be noted. In the classical approach, the 
test section walls are extended an Infinite distance upstream from the model 
to the reference f4ow. For perforated walls, the establishment of the 
reference flow is dominated by a pressure balance condition, even In the limit 
of zero porosity, rather than by the mass conservation condition which 
dominates the solid wall case. As a result, the predicted perforated wall 
wake blockage Increases negatively with decreasing porosity and reaches a limit 
at zero porosity which Is equal In magnitude but opposite In sign to that 
for solid walls. This result Is clearly Inappropriate for real perforated 
tunnels In which the reference pressure Is measured much closer to the model. 
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CONCLUDING REMARKS 


The method described In this paper for assessing the wall Interference 
in transonic two-dln^nslonal tunnel tests embodies two significant departures 
f|^m the methods generally used for wall Interference prediction. The use 
of experimentally measured boundary conditions provides the means of accounting ' 
realistically for a number of tunnel Interference phenomena which are either 
overlooked or represented less accurately by other methods. The use of a 
nonlinear transonic potential equation to represent the tunnel flow assures 
applicability of the method to supercritical test conditions. The method yields 
corrections to the test Mach number and angle of attack and also yields a 
measure of the adequacy of nhese corrections to account for the Interference 
existing during the test. The Improved accuracy and realism Illustrated In 
this paper should lead to Increased confidence In the use of airfoil test 
facilities, particularly at supercritical speeds. 
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Figure 2.- Comparison of wall-induced velocities calculated at tunnel 
center line using measured pressures with predictions by clasale«i 
approach. NACA 0012-64 airfoil In University of Southampton tunnel 
with straight solid walls; c/h • 0,9; M ■ 0,1, 



Figure 3»- Effect of wall porosity on airfoil pressure distribution. 
10-percent thick supercritical airfoil In Lockheed CFFj c/h « 0.25; 
M ■> 0.8; design c^. 



Figure 4.- 


Effect of far-fleld Mach number update 
velocities. T » 0.013; M *> 0.8; 


on calculated wall-induced 
design c„. 
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RESEARCH ON SELF-CORRECTING WIND TUNNELS* 

R. J. Vidal and J. C. Erickson, Jr. 
Calapan Corporation 


SUMMARY 


«ur^sS;.if3S«r‘f T” perf0""«> at conditions where the “ Ihl 

conditHL ‘^Th! procedure is described for these 

»II5*e. * results Of an analysis of the flow in the auxiliary suction 
y tem and test section illustrate the trade-offs available In the nf 

self-correcting wind tunnel test sections and in^SSl 


INTRODUCTION 


research is in progress at the Calspan Corporation to de- 
pa seif-correcting wind tunnel and to demonstrate that interferencA ^ 

iolM°^ measuring the cpmponents of the disturbance veloJity^at^disSete^" 

^®fi"ary control surfaces, or interfaces, in the flow field with- 
in the tunnel. A theoretical formulation for the flow «Aid ! e.u “ 

t^raiVSstJbi^ct s::; 

pSvS?^-"‘« 

Snsl*5”?*Vh*l^®”"°°"^ 'aad“s?ed‘m?irthS‘M2“jS^^ 

condition for unconfined flow. In\his way the 
best theor^ical and experimental features are combined to minimize wall*inter 
ference. The flow chart for the self-corr ecting tunnel schte 

» Research sponsored by the Office of Naval" Research and the Air Force Office 
of Scientific Research under Contract No. N00014-72-C-0102, with supolemen- 

SHE =SSsS:= 
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figure 1. 


which^hara*"solL’‘ML?®®" “ two-dimensional airfoil model, 

field was euporcritlcal[!'LftL'£L':rtr«ur:.r«^^^^^ S“" 

trat“thS‘wairi!!Jrf°'’*‘'‘ "Viewed‘h«e to ill”. 

to trade-offs which can be made in t.st*“«K:M“S?^ral:/“"“' 


SYMBOLS 


fnci Customary Units 

and calculations were made in U.S. Customary Units. 


The measurements 


c airfoil chord length, m (in.) 

®l,^d,^m ®itfoil section lift, drag, and pitching-moment coefficients 

Cp pressure coefficient 

M Mach number 

U stream velocity, m/sec (ft/sec) 

(ft/^ecf^** streamwise components of disturbance velocity, m/sec 


X 

€t 

A 


streamwise coordinate with origin at airfoil leading edge, m (in.) 
airfoil angle of attack 


(1-Mi 


) 


1/2 


Subscripts : 

c calculated value 

ni measured value 

oo free-stream conditions 


\ 
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THE CALSPAN SELF-CORRECTING WIND TUNNEL 


Tho Calspun 0>30 n (1 ft.) Solf-Corroctinu Wind Tunnoi 9 nn^i 

1 = a continuous-flow facility/ Tho tost sLtlon irtio TLLSil „U nor- 
ous top and bottom walls of 22.5% open area. Tho plenum chambers behind tho 
porous wails have boon Jlyldod Into 18 sog.,onts^ ironh^wfandT™ tw 
bottom, and each segment Is connected to a pressure and a suction source 
through individual control valves. The pressure source is the tunnel stilline 

auxiliary compreLi? dlL^iig 

iJeli?;; diffuser. Six plenum chambers in the imSlXJe 

^ provisions for a distributed porosity which cm be 

varied linearly in the streamwlse direction. ^ 


0012 section are shown in figure 4. The model is an NACA 

12 airfoil section with a 0«1S m (6 in*) chord (6% solid blockaae*) anA 

"iJ’’ “ “* P’^'ssure orifices and with a three-co»poient%orce 

talance supporting a metric section on the tunnel centerline. The test sec- 
tion instrumentation (fig. 4) consists of two static-pressSr; pipes! eacHith 

»uS Krr”d':ia!!: Tf's;! 

strumentation are given in references 3 and 4, * 


EXPERIMENTS 


• experiments were performed with the 0.15 m (6 in.) chord model 

solid blockage) to icq^ire a 

^ ^ regarded as essentially free from wall interfer- 

ence. These data are presented and analyzed in reference S. 


Initial experiments with the model in the 0.30 m (1 ft.) tunnel were made 
to simulate tests in a conventional porous-wall facility This simuiai-ion 

t1o:la°fh“^‘"‘."““ “ tmlfo™ «!!:i p^slu^HumSu! "■ 

Jir.™ ‘ftt section. The model was then installed and tested with 

the same valve settings. The lift, drag and pitching-moment results obtained 
at Moo= 0.725 are compared with the 2.44 m (8 ft.) tunnel data in fieure 5 

Comparison shows that there are wall -interference effects’on 
the model forces Md moment in this simulation. The effects on lift are mod- 
erate and on pitching-moment are large, and each is qualitatively consistent 
interference for lift coefficients greater than 0.1. The 
effects on drag are appreciable and indicate open- jet interference for sub- 

uJ d«! “r-a !f P"OUs-;al/theo«uSr“nd o^JortL- 

tal data (refs. 6 and 7) are applied to these cases, they would predict wall 
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Intorforonco Th« data 

tMa probably rafiocta tSHaat tLri°t"' 

categorize these results within ‘oversimplification to nttemn? <-n 

latarfaronco. Tba latarfa^aa" i„\luda“^* f?rr'“" or oJ™”Jot 

w““?a?ti!' ff »'>“*k-wavfSaftloa‘M'‘' "“l™- 


Iterations for Subcritical Walls 


cxporiments wero porformeH in 4>u^ 

the model at <t . 2 * and M_. o 725 ”°'*® '>)' iterating with 

Static pressures ai i 4 * 1 . ®*«wise disturbance velocities i 1 « 

functional relationships Jo recorded and used in evaliatiM'tir®^ 

o‘si.:;roTHS;^^ 

points, where each^nSbe^lJdicaterthf*'?!™ ‘"/‘S®’'® 5 by the numbered data 

‘be equivalent daw point iS tL I?m T'/a'*'’' ™® ‘““-I «Ung?r 
It can be seen that after three lt«i.ad? 4 ^”® simulated conventional facility 

interference effects in that tL liA inTAT eliminated thJ wall 

!^‘‘fl data, within the scatter in that 2.44 m (8 ftO 

l^ormally, we would have moment agrees 
the evaluation of the functi^iarreJ^tioJsMnf the iterations bfcauL 
vergence, but we were at the limite nf ^ips did not indicate complete cnn- 
model and could not achieve the next iter2tf°^ °"® ^1®"“® section above the 


Operating Procedure for Supercritical 


Walls 


n T/m ■^Jr^’^^ments were performed with 

0.30 m fl ft "i tlinnrwl „ 4 . w with (,„o U.ia 

d* :.^^“""®1 ®t Moo“ 0.8S 


- V* tunnel at M^o n ar >» „ io .. tnora model in the 

cohtrol . be^talnf 
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streamwise disturbance velocities. We were able to 
!? desired distributions upstream of the model quarter chord but 
flow downstream of the quarter chord The result in0 
flow field is illustrated by the schlieren ph??oJrJJh i^figuS 8 

shwk ^**® contained multiple 

SHOCK waves in the remainder of the test section. 

t*»e desired flow-field control by modifying our 
trol at J ^**® **®®* technique is to establish con- 

irinSeLfttaTeJ'^nSw ‘^5 "® ®“*>c^iticsl» and then sequentially 

tniltt distributions of the streamwise dinurtaJle ^^ScuLs. 

In this particular case, we established flow-field control by settine the first 

wSrSutTOrt^thn^* J^^’^aasad the Mach number to 0.85, readjusted the 

shiin ™ o !wf Obtalnol th. flow field 

snow in tigitre 9* The shock waves are close to the correct oosition 

evid^ced by the 2.44 m (8 ft.) tunnel data, and the walls LnuJwcktLal 

from^th«*^^?^«.?^ '^®^® ®®anating from the flow-angle probL The data 

ftout 1 05 S“* "®" *■'0 wall was 

pu):L‘s^iioS rftijs^srof ''-® 

CHARACTERISTICS OF THE AUXILIARY SUCTION SYSTEM 

The recurring inability to achieve full control of the flow field at all 

«‘® «<* i” the LSLar? section smL 

lights and conclusions of the analysis are given here* 

analysis consists of writing a pressure balance for the 
auxiliary suction system, which is shown in figure 3. The analysis includes 
the operating characteristics of the auxiliary compressor, the losses in the 

iit^th^ diffSeS’'®th^°" pressure loss arising from the compressor dis^haJge 

presswe drop across the porous walls, and the required 
wconfined-flow pressure at the wall locations due to the model. When suitable 
enpneenng appox^ations for these pressure terms are made, using, in part 
data measured in the 0.30 m (1 ft.) tunnel, the predicted limitrSS aJaiable 
control agree reasonably well with experimental Observations a'^ailable 


coi.id°h! !!5 ^ analysis was that the recompression penalty 

^ reduced considerably. This could be accomplished by introducing an 
area change in the tunnel diffuser at the location where the flow from the 
auxiliary blower is vented into the tunnel circuit. The best arrangement is 

«>'« would bnu«ion“n the 
«“ 'l»"8b modification promises a considerable Jm- 
provement in the circuit performance and has been carried out. 
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The analysis also illustrates the three major trade-offs available in the 
design and application of self-correcting wind tunnels with porous walls; 
namely, compression ratio, porosity, and model size. The results in reference 
8 show that the tunnel performance could be inq>roved by increasing the compres- 
sion ratio from 1.43 to 2.0, but the in^rovement would be available only at 
Mach numbers above about 0.75 because the wall perforations would be choked at 
Mach numbers below that value. That restriction could be relaxed considerably 
by using a higher porosity, although the shock-wave reflection characteristics 
may not be as favorable. Decreasing the model size from 6% blockage to a 
smaller value would decrease the magnitude of the disturbance velocities at the 
walls and would inqprove the tunnel performance. The latter alternative, de- 
creasing the model size, has been selected for the Calspan facility and itera- 
tion experiments have begun at higher Mach numbers with a model having 4% 
solid blockage. 


CONCLUDING REM/^RKS 


A two-dimensional self-correcting wind tunnel has been built and tested 
at the Calspan Corporation to demonstrate that wall-interference effects can 
be minimized by actively controlling the flow in the vicinity of the tunnel 
walls. Experiments were made with a 6%-blockage airfoil model, using the 
facility to simulate a conventional wind, tunnel, to show that the wall-inter- 
ference effects can be large in such tunnels. Several cases have been iterated 
successfully in the self-correcting. mode. In one supercritical flow over the 
airfoil, we iterated three times and were able to reproduce the correct shock 
position, eliminate wall-interference effects on lift and drag, and decrease 
them to 5% to 7% on pitching moment. Experiments with supercritical walls 
showed that it is necessary to establish wall control initially when the flow 
at the walls is subcritical, and then sequentially to increase the Mach number 
and readjust the wall control until the desired test condition is achieved. 

We have devised a method for analyzing the performance of a self-correcting 
test section with porous walls. This method makes it possible to examine the 
trade-offs between compression ratio, wall porosity, and model size. It should 
be a useful tool in future design studies of self-correcting wind tunnels. 
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Figure 2.- Cutaway sketch of Calspan two-dimensional self-correcting 

wind tunnel. 
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Figure 6.- Simulated conventional wind tunnel. 
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Figure 8.- Schlieren view of the flow field with incorrect wall control. 
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ANALYTICAL DESIGN OF A CONTOURED WIND-TUNNEL 
LINER FOR SUPERCRITICAL TESTING 

Perry A. Newman and E. Clay Anderson* 
NASA Langley Research Center 


SUMMARY 

The present analytical design procedure Is being developed In order to 
determine the shape of a contoured nonporous wind-tunnel liner for use In the 
Ames 12-foot-pressure-w1nd-tunnel test of a large-chord, laminar flow control 
(LFC), swept-wing panel which has a supercritical airfoil section. This 
procedure Is applicable to the two-dimensional streamlined tunnel problem and 
a first check on Its validity would be a comparison of the calculated tunnel - 
wall shape with that found experimentally. Results for such a comparison are 
given and the favorable agreement Is encouraging. 


INTRODUCTION 

The analytical design procedure discussed herein Is being developed In 
order to meet the special requirements of a laminar flow control (LFC) test. 

A transonic test condition Is needed In order to establish the compatibility 
of an active LFC wing suction system with the current high-performance super- 
critical airfoil technology. This test must be done in a wind tunnel which 
has low levels of stream turbulence and acoustic noise so that the laminar 
flow Is not unduly disturbed. Conventional slotted- or porous-walled transonic 
tunnels are, therefore. Inappropriate In this regard. Transport aircraft 
presently envisioned for LFC applications have moderately swept wings; thus, 
the laminar boundary- layer crossflow stability must be Investigated at the 
appropriate flight crossflow Reynolds number. This requirement, together with 
the physical-size limitations set by slot/duct construction and the required 
roughness-height Reynolds number for laminar flow, results in a large-chord, 
swept-wing panel. In order to produce a transonic wind-tunnel flow which 
simulates free-air flow about an Infinite yawed wing, one must contour all 
the nonporous bounding walls. The ratio of total tunnel height to chord for 
this test Is about 1 for the Ames 12-foot pressure wind-tunnel. 

The design procedure Is based on a simple Idea and several existing 
computational tools which make It feasible. Basically, one determines bounding 
streamlines in the desired flow, makes all required blockage corrections, and 


*E. Clay Anderson, a consultant, has performed the boundary-layer analyses 
for this study under NASA contract NASI -14517. 
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then builds the resulting contoured ncnporous tunnel liner. At supercritical 
flow conditions, one must use an appropriate transonic analysis In order to 
determine the desired flow-field and bounding streamlines. In order to account 
for the blockage due to viscous effects on the liner, a boundary-layer analysis 
must be made along the bounding streamlines with (local) edge conditions 
determined by the local flow properties. The bounding streamlines are then 
displaced to relieve the viscous blockage. This displaced shape Is the 
geometric nonporous contour to be built. The liner Is restricted to a given 
airfoil configuration at a single design-point condition. Currently available 
two-dimensional transonic analysis codes (refs. 1 and 2) and a laminar- turbulent 
Jboundary- layer code (ref. 3) have been used. The procedure, however, has a 
general utility which Is restricted by our current ability to calculate the 
viscous transonic flow field about arbitrary configurations. 

This design procedure is, of course, directly applicable to the two- 
dimensional streamlined tunnel problem. One should be able to predict the 
tunnel -wall shape required to simulate free-air flow about an airfoil. A first 
check of Its validity Is provided by a direct comparison of the analytically 
determined tunnel wall shape with one determined experimentally In a stream- 
lined tunnel experiment. Barnwell and Everhart 1 have recently completed such 
an experiment for a symmetric airfoil at zero lift In the Langley 6- by 
19-Inch transonic tunnel. In that test, the slotted walls of the tunnel were 
replaced by nonporous flexible plates which were adjusted using an experimental/ 
analytical procedure until a simulated free-air flow over the airfoil section 
was obtained. Thus, the tunnel -wall displacement required to produce equivalent 
free-air flow Is determined as a result of the experiment. Their Iterative 
experimental /analytical technique Is similar to that presented by Goodyer 
(ref. 4). 

This paper should be viewed as a progress report In the continuing develop- 
ment of the design procedure for the LFC swept-wing test application. Some 
(preliminary) test-section contours which were generated to check the numerical 
procedure are given In the next section. The following section summarizes an 
application for a two-dimensional streamlined tunnel experiment and compares 
present analytical predictions and experimental results for tunnel-wall shapes. 


SYMBOLS 

c airfoil chord length, cm (Inch) 

M free- stream Mach number 


U streamlining procedure (unpublished) has been developed by Richard W. 
Barnwell and Joel L. Everhart of Langley Research Center. The present analyt- 
ical results are compared with the wall shape experimentally determined by 
this procedure. 
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^Re,“ free-stream Reynolds number based on airfoil chord 

X streamwise tunnel coordinate measured from tost section origin 

indicated {Q>) in figure 4, cm (inch) 

y lateral tunnel coordinate normal to tunnel span, measured from test 

section origin indicated (©) in figure 4, cm (inch) 

a angle of attack, deg 

«*(x,y) local boundary- layer displacement-thickness correction, measured 
normal to tunnel wall, cm (inch) 


LFC SWEPT-WING PANEL TEST LINER 

Current plans for the LFC swept-wing panel test are to utilize the Ames 
12-foot pressure wind-tunnel to test a large-chord model, which has a cambered 
supercritical airfoil section, at lift. Recall (as shown in fig. 1) that for 
flow over a lifting yawed wing the streamline which divides around the leading 
edge of the airfoil develops a different spanwise displacement on the two 
wing surfaces. The resulting displacement at the trailing edge persists down- 
stream and one must contour the end walls to properly simulate the inviscid 
crossflow. 

The LFC airfoil section design and its aerodynamic characteristics are 
for a given set of wing-normal conditions. (See ref. 5, for example.) However, 
as indicated in the last section, considerations of crossflow Reynolds number, 
unit Reynolds number, and slot/duct construction bear on the selection of a 
model chord length and test conditions. These are certainly not independent of 
the tunnel operating characteristics and it appears that the required conditions 
may well be near the tunnel power limitations. All of these aspects result 
in a ratio of tunnel height to chord of approximately 1. The contraction ratio 
of the tunnel is 25 to 1 and an appreciable increase could easily result in 
unsatisfactory diffuser performance. A liner of ''octagonal'' cross-section 
results in less additional contraction from the existing circular cross-section 
than would result from a "rectangular" liner. Figure 2 is a view of the tunnel 
showing an "octagonal" cross-section liner and the large-ciiord swept-wing panel. 

All walls are contoured in order to simulate free-air flow about an 
infinite-aspect-ratio yawed wing at the design condition. A preliminary design 
of an "octagonal" liner shape has been made, considering only the inviscid 
streamline shapes, in order to check the numerical procedure. This was done 
using an early LFC supercritical airfoil design of chord length 2.44 m (8 feet), 
swept 35°, at a tunnel Mach number of 0.89. Schematic views of this contoured 
liner shape near the LFC model are shown in figure 3. The section views 
clearly show the development of the end-wall steps which persist downstream of 
the trailing edge. It is also seen that some of the streamlines fall outside 
the existing tunnel in places for this first try. A turbulent boundary layer 
will exist on all tunnel walls and the local viscous blockage correction for it 
has not been included in this shape. This correction will move the liner even 
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further outside the existing tunnel In places. The sizing will ^ . 

Iteratively In order to minimize the additional contraction. 5"^ ^ 44 .<rtn 

fSre 3 shows the distorted flow which will have to be reshaped 
section In the liner In order to fair Into the existing diffuser. This fairing 
will be shaped to minimize the degradation of the existing tunnel pet forman . 


2-0 STREAMLINED TUNNEI WALL SHAPE 

The liner design procedure being developed for the LFC 
applicable to the prediction of wall shapes for two-dimensional streamlined 
tunnels. Details concerning application of the present approach to a two- 
dimensional tunnel are given In reference 6. The results from that study are 
summarized In this section. 

In the Barnwell -Everhart experiment the slotted lateral 
Langley 6- by 19-Inch transonic tunnel were 

walU%s Indicated In fig. 4) In order to Investigate the streamlined tunnel 
concept for minimizing wall Interference In transonic tests. 
their experimental/analytical procedure iterates 

until consistency is obtained between near-field and tar-field solutions. The 
wind tunnel determines the near-field solution for a Qj^en wall shape (I.e., 
exoerlmentally "solves" a direct, viscous, boundary-value problem). The 
Inverse far-fleld solution corresponding to uniform parallel flow at if'tinity 
Is determined numerically using the experimental 

to predict a new wall shape. Thus, the tunnel-wall ^JaP® i% detemlned as a 
result of the experiment. Their procedure has been used to tijj<i Jhe wall p 
for only one model and test condition. However, the first ^^ep In their pro 
cedure Is to determine a wall shape which produces uniform pressures In the 
test section for the empty tunnel (i.e., no model). 

For the empty tunnel the lateral wall shape required to produce uniform 

test-section pressures Is a weighted 6Mrv)°"m1nus 

to all of the tunnel wall boundary layers* That is> the local o (x^yh minus 

that at the upstream flexible-wall/rlgld-wall junction. ^ ^ J 

the perimeter of the tunnel cross section to obtain a 

which is then applied at the two adjustable walls. Comparison of the analyti 
cally determined wall shape and that found experimentally by ^ 

Everhart is shown in figure 5. The tunnel was operating at a [tech number of 
o! 902 and a unH Reyn^ number of 2.24 x 105/cm 3.41 x l06 based on jnodel 
chord). The maximum difference In the wall shapes Is not Inconsistent with the 
pressure measurement accuracies. (See ref. 6.) 

The model test was conducted at a Mach number of 0.765 using ® ]5.24-cm 
(6-1nch) chord symmetric NACA 0012 airfoil section at zero lift. This size 
model results in a ratio of tunnel height to chord of approximately 3* There 
is a turbulent boundary layer on all of the walls and the model. As Indicated 
n flSSre i! onrthrilaX lateral tunnel walls were flexible and these were 

deformed to simulate free-air flow about the model °"®. Ill®hp°5hl 

is centered at the location (©) Indicated in figure 4; this is J®ken to be the 
origin for streamwise and lateral coordinates in what follows. Each flexible 
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wall is supported by 11 adjustable screw jacks which are located 10.16 cm 
(4 inches) apart. The streamwise extent of the flexible-walled section is 
from 4.333 chords upstream of the leading edge to 2.833 chords downstream of 
the trailing edge. 

Calculations were made along 48 free-air streamlines about the NACA 0012 
airfoil in order to estimate the tunnel-wall boundary layer in the presence of 
the model. Boundary-layer displacement-thickness distributions across the 
tunnel rigid walls (at several streamwise locations) are shown in figure 6. 
These clearly show a significant local influence due to the pressure field of 
the airfoil. Note that the y/c axis is located at a value of 6vc » 0.014 on 
all plots; the maximum value of 6*/c (shown on the last curve in fig. 6(a)) 
is 0.03. At the lateral wall (largest y/c shown) the displacement thickness 
exhibits a monotonic Increasing behavior. However, along streamlines near the 
airfoil surface, 6*/c undergoes very rapid streamwise changes. There is a 
pronounced local boundary- layer thickening near the leading and trailing edges 
of the model with a thinning near the expansion peak just downstream of the 
leading edge. The distribution in the wake region may not be very realistic 
since the wake is not properly modeled in the transonic airfoil program. 

The flexible wall shape is determined as the sum of the deflection of the 
outermost lateral streamline (i.e., the inviscid compressible blockage correc- 
tion) and the viscous correction due to the local boundary-layer displacement 
thickness on all walls (i.e., flexible and rigid). The method for determining 
this latter 6* correction was indicated when discussing the empty tunnel 
results. Comparisons of analytical and experimental flexible wall shapes are 
shown in figure 7. The agreement is seen to be very good. The outermost 
lateral inviscid streamline is shown as a dashed line so that the relative 
contributions to the wall displacement can be seen. The viscous wall effect 
is not small; the boundary layer on the rigid wall near the model produces a 
rapid streamwise variation in blockage which must be accounted for in the 
flexible wall. One must have sufficient resolution in the wall control at 
lateral regions around the model in order to accommodate this behavior. 


CONCLUDING REMARKS 

The present design procedure has a general utility which is restricted by 
our current ability to calculate transonic flow about arbitrary configurations 
including a proper account of the viscous phenomena. However, continuing rapid 
advances in computational machinery and methods will allow for better 
numerical modeling of these aspects. Comparisons with the Barnwell -Everhart 
streamlined tunnel experiment provide some degree of verification of the 
numerical design procedure which is being developed for the LFC application. 
These comparisons demonstrate that the blockage due to the viscous-boundary 
layer on the rigid, (side) wall cannot be neglected in tunnel streamlining. In 
fact, it is primarily this wall boundary layer near the model which responds 
to the pressure field of the model and produces the rapidly varying streamwise 
character in the lateral wall shape. For more severe conditions of trcnsonic 
flow with strong shockwaves and high lift, it is felt that the local rigid- 
(side-) wall boundary-layer effects will distort the test so that the desired 


simulation is not achieved (i.e., 2-D flow will not be achieved in the tunnel). 
The present design procedure contains an option for direct calculation of 
suction rate distributions required to maintiin attached boundary-layer flow. 
This option provides a means for determining suction requirements needed for 
boundary- layer control near the model/liner junction. 
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Figure 2.- View of the Ames 12-foot pressure wind tunnel showing "octagonal” 
test-section liner and large-chord swept-wing panel for LFC test. 





SECTION A-A SECTION B-B SECTION C-C SECTION D-D 


Figure 3.- Schematic view of test-section liner for an LFC model in the Ames 
12-foot pressure wind tunnel. Liner was contoured for a preliminary LFC 
section neglecting viscous wall effects. 

RIGID WALLS 
SECTION A*A 

V DENOTES WAa JACK LOCATION 

Figure 4.- Schematic of two-dimensional streamlined tunnel (Langley 6- by 
19-inch transonic tunnel) used by Barnwell and Everhart. 
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= 0.765, a = 0°, = 3.05 x 10 

O EXPERIMENT (BARNWELL-EVERHART) 

— CALCUUTED WALL SHAPE 

— CALCULATED STREAMLINE 


Figure Comparison of flexible nonporous wall shapes for an NACA 0012 
airfoil test at zero lift in the Langley 6- by 19-inch transonic tunnel. 





N79-19994 

PRECEDING PAGE BLANK NOT FILIUEO 


34 


l‘‘A 
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6 IN. BY 22 IN. TRANSONIC AIRFOIL TUNNEL 


John D. Lee 

The Aeronautical and Astronautical Research Laboratory 
The Ohio State University 

. iJiNG PAGE BLANK NOT FILMED 


SUMMARY 


Interference in a 6 in. by 22 in. "two-dimensional" wind tunnel has been 
evaluated at Mach numbers up to 1.06 by comparing pressure distributions from 
airfoil models of differing size, Mcdels of the NACA 0012 profile having chords 
of 76, 152, and 305 mm were used in one phase of the evaluation program and 
models of a supercritical profile having chords of 76 and 152 mm were used in 
another. This report documents the confinen»nt interferaice, l.e., blockage, 
doTOiwash and streamline curvature, all of which are quite small on a model 
having a chord of 152 mm and which can, for most applications, be ignored. 
Specifically the corrections are lumped into an attack angle adjustment of -0.16 
degrees per unit lift coefficient on a 152 mm model. 


INTRODUCTION 


The interpretation of wind txmnel data in terms of equivalent "free-air" 
Information has always been hindered by considerations of the constraints of 
the boundaries. For some situations, theoretical analyses, moderated by experi- 
ence have yielded suitable correction procedxires. In the transonic range two 
features have led to serious difficulties: (a) the use of ventilated walls 
Introduced peculiar constraints which are not yet properly understood and 
(b) local shock waves interacting with the boundary layers on the txmnel walls 
will often disrupt those layers and produce extensive lateral disturbances. 

For the past 25 years in which tests have been conducted in the transonic range, 
there has been a multitude of innovations but there is yet no simple, direct 
procedure by which valid "free-air" data may be acquired from a groxmd test 
facility other than by using models vdiich are relatively small in a given 
facility. 

Due to peculiarities in design and operation, each facility has interfer- 
ence effects which are \mllke those in another, with the result that proced- 
ures for testing and for modifying the data in one facility will not, in 
general, apply to another facility. It is necessary to examine three areas; the 
facility design, the data acquisition technique and the procedure for deter- 
mining the nature and degree of wall interference, in order to adjust the data 
to equlvtlent "free-air" information. 
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airfoil t^el was designed for the qjeolflc task of airfoil test- 
ing and Its configuration was based upon previous research with three other 
tra^onlc tunnels. A wide variety of airfoils has been tested to date and a 
series of careful programs has been conducted to define and evaluate Interfer- 

describee results from tests utilizing models of 
the NACA 0012 prof ne to define the Interference from the confining boundaries, 

^ similar evaluation program was conducted using an 
aft-loa^d or supercritical" profile to support the results with the 0012 
mdei. The progiram was based on the premise that final reduced and/or 
corrected" data from similar models differing In chord by factors of two and 
fow i^st agree with each other. As a further test, such data must also agree 
with theory, to an acceptable degree, and with data from other facilities which 
were relatively much larger In comparison to the model size. 

The symbols are defined in an appendix. 


THE 6" X 22" WIND TUNNEL SISTM 


airfoil wind tunnel is shown In figure 1. The 

fWllon-""! T ^he ^ ® two-diSional Sn- 

* are solid while the 6" walls are per- 

tS porosity of about 7% and overlie individual plenum chambers! 

c2 ?? zone downstream of the test section and thus 

tSn pressure changes caused by models in the test sec- 

Plenm configuration leads to the extremely low interfer- 
^ experimental data for two-dimensional airfoils over a wide 
range of operating conditions and model attitudes* 


^ mounted in the test section between two ports in the 22" side 
section, and the angle of attack of the airfoil can be 
^ consists of two solid aluminum bSote 

IsmS continuous first and second derlyatl?es Iran- 

is at tl» nosale salt. A pressure drop deulce and a tmi-staged bell^™ 

^ to 

*^®®J humber is fixed by a choke consisting of an array of bars 
be vfriS Induction/mixing zone. The Mach number can 

Pfti number and/or the diameter of the bars. With a typi- 

ance S mmheT can be varied from 0.3 to 1.1 with a tSer- 

Reynolds number can be varied by changing the total pres- 

fT! flYPrt providing a range In Re/m from 1 to 100 million 

tin? This particular feature Is extremely Impor- 

^?? Maoh/Reynolds number effects on two-dimensional airfoUs! 
The facility operating envelope Is given in figure 2. uxnoxxs. 

t^el Is operated by pre-setting the control valve (fig. 1) to achieve 
reservoir pressure. The plug valve is opened and thf ci^ult pres- 
izes In a few seconds. Thereafter, the pressure drops in proportion to the 
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mass flow from the storage tanks. The storage system provides 140 oublo meters 
of air at pressures up to 18 x 10^ N/m2, Conventional air driers and oil 
separators are used to maintain a very high purity of the airflow. 

Figure 3 shows the pressure history for the facility during a typical test 
run. Data are normally taken Immediately after the pressure peak when a period 
of nearly constant pressure Is available. For the conditions of figure 3, 
tunnel shutdoVfn was delayed for purposes of evaluating the response of the 
overall system, where the good uniformity of pressure ratio Is evident. 


Calibration 

The basic characteristics of the tunnel were determined by means of static 
pressure distributions through the test section. The calibration of the tunnel 
was based i 5 )on the average of seven static pressxn*es at the model-mounting site 
(with no model present) referred to the plenum pressure. Calibration tests 
were also perfomed with one plenum vented In order to evaluate the effects of 
the small imbalance which occurs when testing with a lifting model. 

The vedldlty of applying such an ”enpty- tunnel" calibration when a model 
Is Installed Is based on the assunptlon that the plenums, through the perfora- 
ted boundaries, supply the necessary environment to simulate the far-fleld. 

The assunqptlon has been verified by the correlations of the data. 


Data Acquisition and Processing 

Steady-state pressure measurements are made with a Soanivalve and a trapped- 
volume system. Tubes from model and tunnel taps are ducted to a Scanlvalve as- 
tern, consisting of a set of valves (basically rotary guillotines) preceding a 
regular Scanlvalve, At the data point, the valves are closed, thus trapping 
seuiples of all the pressures at the same tunnel operating pressxire. The pres- 
sures are subsequently read through the Scanlvalve as the tuimel is shutting 
down. Due to tte small but finite volume in the transducer cavity, adjustments 
to the raw data are necessary to accoimt for pressure differences as Ihie Scani- 
valve is cycled between ports. This is accomplished with simple calibration 
procedures, during which the Scanlvalve Is cycled two or more times. The multi- 
cycle calibration is frequently repeated during a test program as a cross check 
on the system and for detection of leaks. 

Data taken during the run as well as those from the Scanlvalve are digiti- 
zed and stored on magnetic tape In the Digital Computational Facility of the 
Aeronautical and Astronautlcal Research Laboratory. The data Is reduced to co- 
efficient form, reviewed on a CRT display, then Integrated to lift and pitching 
moment coefficients, printed and plotted in hard copy form in five minutes from 
the run. 

The wake Is traversed by a motor-driven pitot probe, 25 cm downstream from 
the trailing edge, and the pressure recorded as a continuous trace of deficit 
from the stagnation pressure. Again, this datum is stored during the run cycle. 
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Pressure tape are Incorporated by blind diSlHne"^'"ri,°f coordinate seta, 
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. ebows a ^el In tbe tesi 


INTERFERENCE EVALUATION 


I. Flow Quality 

stat^ pressure mea^SeSs°through\S°Lst*se^ <ietennlned by means of 
are given In figure 6. section; two such sets of data 

Proce<i^^^y S SdJSrfrS°datrsuch S^th the call- 

and B, For such purposes the mi? P^®sented In figures 7 

be noted that a i^lS?v2y ^ ® It should 

of agreement. ^ number will shift the data out 




attao”S2L^?lJe°^?i'’of^harh'"®““*® 7°” ^r° 

enoee between the three modeS ‘l>a email dlffer- 

in the direction to be e^Sd fS^\o»d\if 








been Scui^tenJ m^ai^*S‘'S%™ thJe; 

from two models of a 3upercrltl^aiXn^„^S^!t "“'’S® °°12 Pf“«le and 

ranges In Mach number ^ attack Lke^'^IntP^r ^°th tested over wide 

sons of the data between models with da+n rr* niay be deduced by coinparl- 

retlcal predictions ( l^lSlS vl^c Juf “« ’'I**' 

examples are presented hers In figures 10 thrmi^h o®^v ^ I'®'"' 

the conclusion that the interfereSe on a 15 ? comparisons lead to 

the full range of Mach number and In faf.+ ®°^el Is negligibly small over 

305 mm chord In subcrltlcal flows' ®^®“ ® “O'!®! cf 
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A conoldoraWe number of data were uaed to determine the Interference 
field In flome detail. For example, the load differences plotted In figure 14 
Indicated that the Interference could be separated into Its components due to 
downwaoh and to streamline curvatxrre. Further data were used to establish the 
dependence of both components on the lift coefficient and on the square of the 
chord-to-helght ratio. The Interference field could then be expressed In 
general form as given by figure 15. It should be noted that the scales of 
figures 14 and 15 are higher expanded for analytical purposes. 

Results from the supercritical models were essentially similar to those 
from the models of the 0012. Figure 15 also shovra that the chordwlse Interfer- 
ence load on the supercritical section was the same as that found on the 0012. 

For comparison purposes, the lift Interference Integrates to a correction 
on the attack angle of -0.16® per unit lift coefficient on a model having a 
chord of 152 mm. 


CONCLUDING REMARKS 


Interference generated by the proximity of the walls in the 6 In. by 22 in. 
airfoil tunnel has been evaluated as extremely small by means of data frcm 
three models of the NACA 0012 profile and from two models of a supercritical 
section. Using a 152 mm model chord for reference, the net lift interference 
may be expressed as an attack angle error of -0.16® per unit lift coefficient ige 
and the blockage as negligible. 


APPENDIX 


C 


P 

Po 

P» 

Pc 

Re 

X 


^“dw 


AC. 


sc 




model chord length (cm) 
lift coofflolent 
jiressure coefficient 
wind-tunnel height 
Mach number 
pressure 


stagnation pressure 
free-stream static pressure 
pressure in plenum of wind tunnel 
Rejmolds number 


chordwlse distance on airfoil (cm) 
angle of attack (degrees) 

airfoil loading (Cp (lower surface)) - (c^ (upper surface)) 

angle of downwash induced by wind-tunnel interference (degrees) 

of lift from streamline curvature induced by 
wind-tunnel interference oy 
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Figure 7.- Pressure distributions from the 152 mm 0012 model conq>ared in a 

Prandtl-Glauert form. 



« - DEGREES 


Figure 8.- Integrated lift coefficients for the 0012 models compared in a 

Prandtl-Glauert form. 
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Figure 10.- Pressure distributions from the 0012 im>dels con^ered to show 
Interference et a sub critical condition. (Data from the 152 mm and 
76 mm models fall within the same symbol on thi > scale.) 
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Figure 11.- Pressure distributions from the 0012 models compared to show 
Interference at a supercritical condition. 
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12.— Pressure distribution from the 152 tun 0012 model compared with 
that from a 152 mm model In the CALSPAN tunnel. 



Figure 13.- Pressure distribution from the 0012 models compared at an 
extreme condition In the 6 In. by 22 In. tunnel. 





x/c 

Figure 14.- Loading interference in the 6 in. by 22 in. tunnel on the 305 mm 
0012 model compared with a superposition of downwash and streamline 
curvature. 
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Figure 15.- Generalized loading intcrferi'Mci' . 'mii.iri-d for the various sets 

of models te.sti iL 
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VISUALIZATION OF THE SEPARATION AND 
SUBSEQUENT TRANSITION NEAR 
THE LEADING EDGE OF AIRFOILS* 


Anthony V. Arena and Thomas J. Mueller 
University of Notre Dame 


SUMMARY 


A visual study was performed tislng the low speed smoke wind tunnels of 
the University of Notre Dame with the objective of obtaining a better under- 
standing of the structure of leading edge separation bubbles on airfoils. The 
location of separation, transition and reattachment for a cylindrical nose 
constant-thickness airfoil model were obtained from smoke photographs and 
surface oil flow techniques. These data, together with static pressure distri- 
butions along the leading edge and upper surface of the model, produced the 
Influence of Reynolds number, angle of attack and trailing edge flap angle on 
the size and characteristics of the bubble. Additional visual Insight Into 
the unsteady nature of the separation bubble was provided by high speed 16 mm 
movies. 8mm color movies taken of the surface oil flow supported the findings 
of the high speed movies and clearly showed the formation of a scalloped 
spanwlse separation line at the higher Reynolds ntimber (Re^ > 450,000). Results 
obtained from these experiments have Improved our knowledge of the laminar 
separation bubble and the transition process. 


INTRODUCTION 


The next generation of business and commercial aircraft must not only 
have significant reductions In the noise and exhaust pollutants produced but 
must be fuel-conservative. One of the crucial areas where significant 
advances have been made In the past few years Is In transonic and subsonic air- 
foil design. Although advances have been made, there are several areas which 
require careful Investigation If further Improvements are to be realized, 
especially In developing an airfoil with laminar flow over a substantial 
portion of Its surface. One Important area of concern Is the occurrence and 
behavior of the leading edge separation bubble. The separation bubble plays 
an Important part In determining the behavior of the boundary layer on the 
surface, including the stalling characteristics of the airfoil. In practical 
situations, this bubble may be partly laminar, with transition occurring after 
laminar separation, or completely turbulent (ref. 1). The primary objectives 


* This research was supported by NASA Langley Research Center under 
Grant NSG 1419. 
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of this paper are to describe the unique low speed smoke wind tunnels and vlS" 
uallzatlon techniques developed at Notre Dame and to use these techniques to 
obtain a better understanding of the structure of the leading edge separation 
bubble and the transition process. 


SYMBOLS 

Values are given In both SI and U.S. Customary Units. The measurements 
and calculations were made In U.S. Customary Units. 

c Chord 

C Pressure coefficient ■ (P - P )/4pU ^ 

p 00' ' 2 

P static pressure on model surface 

Freestream static pressure 

cUj) 

Chord Reynolds number ■ — 

R Reattachment location 

5 Separation location 

T Beginning transition location 

Freestream velocity 

X Distance along airfoil chord 

Distance from beginning of Test Section 

cx Angle of attack 

6 Flap angle at centerline 

6 Diffuser angle, deg 

y Absolute viscosity 

p Density 

EXPERIMENTAL APPARATUS AND TECHNIQUE 

A schematic drawing of one of the Indraft low turbulence subsonic wind 
tunnels Is shown In figure 1. The air passes through the twelve anti- 
turbulence screens before entering the 24:1 contraction section. The first 
seven screens are 14 x 18 mesh bronze and the last five are 25 x 35 mesh 
marquisette. Both the Inlet and test sections are square and are mounted on 
rollers to provide an easy means of Interchanging these components. In 
addition, each test section has removable rear panels, allowing ready access 
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to the test section. Downstream of the test section Is the diffuser which Is 
fixed Into the wall of the laboratory and which transitions from the square 
cross-section of the test section to a circular cross-section at the fan. A 
four Inch foam rubber section separates the tost section and the diffuser to 
ensure no vibrations are transmitted to the working section from the fan. The 
fan is driven by an 11,185.5 W (15 hp) variable speed electric motor, located 
on a separate foundation outside the laboratory. The two-dimensional airfoil 
type model used In this Investigation was the cylindrical leading edge constant- 
thickness model with a movable flap shown In figure 2. Thirty-four static 
pressure taps were placed along the lending edge and upper surface of the model 
ns far as x/c ■■ 0.54. 

For smoke visualization, the smoke Is Introduced upstream of the first 
screen from a rake which can be moved and positioned horizontally or vertically 
so smoke lines may be placed anywhere In the test section. The smoke Is gen- 
erated by a device which allows deodorized kerosene to drip on to electrically 
heated plates; the smoke then proceeds to the smoke rake. The smoke rake has 
a filter bag and cooling colls which reduce the smoke temperature to approx- 
imately ambient before passing through the anti-turbulence screens and Into the 
test section. Still photographs of the smoke flow were taken, using a 
Graflex 10 x 13 cm (4 x 5 Inch) camera and several high intensity lights with 
a duration of about 20 micro-seconds. High speed movies were obtained with a 
Wollensak WF-3 Fastex 16 mm camera and several 1000 and 2000 watt quartz 
lights. Surface oil flow studies were also made, using time lapse and 8 mm 
color movie photography. 


DISCUSSION OF RESULTS 


Turbulence Intensity Measurements 

It was necessary to measure the freestream turbulence In the wind tunnel 
to ascertain the effect of turbulence Intensity on the flow over the airfoil 
model In this and future hot-wire studies. Two stations were chosen at which 
hot-wire measurements would be conducted. The first was near the entrance to 
the constant-area tunnel section (}I^ “ 37 cm) and the second point was further 
downstream coincident to the location of the leading edge of the model 
(X^ « 76 cm). The point “ 0 signifies the point where the reduction cone 
ends and the test section begins. The model was removed from the tunnel while 
measurements were taken at thq second station. Four different sets of 
measurements were taken at each station, varying the velocity from 6 m/sec to 
25 m/sec. Three of the four runs utilized turbulence screens directly at the 
test section entrance (X^ » 0) . Results from these experiments are shown in 
figures 3 and 4. Turbulence Intensities ranging from less than 0.2% (no 
turbulence screen) to less than 1.2% at station X^ = 37 cm (i.e. fig. 3). 

These values dropped 0.1% and 0.7%, respectively, at X^ » 76 cm, as shown In 
figure 4. 
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Siatlc Pronaure Dlatrlbutluna 


Static proaaurn coofllclent plota were obtained for nnRlea of attack 
ranp,lnR from 6° to "4°, flap angina of 0*^, +30° and "30° (poaltlvn chonnr, aa 
downward), and chord Rnynolda numbnra from ISO, 000 to 460,000. Typical proa" 
auro copfflclent vnraua percent chord rcaulta arc ahown In figure 5. Thcac 
prcaaurc coofflclont pl<»ta for the ahcar layer flow may be divided into four 
dlatinct roglona: (1) the laminar boundary layer extending along the airfoil 

aurfnce from the leading <‘dgc to the point of acparatlon, characterised by a 
marked decrease in preaaure to a minimum peak, followed by an adverae pressure 
gradient causing the flow to aeparate; (2) the constant-pressure laminar 
shear layer, existing from the points of separation to transition; (3) the 
turbulent mixing region, between the onset of transition and the point of 
rcattachment, characterized by a quick rise in pressure; (4) the turbulent 
boundary layer stretching from reattachment on downstream. 


Smoke Pictures 

For the smoke flow visualization portion of the experiments, several 
photographs were taken for each condition that pressure plots were available. 
The purpose was to freeze the action of the separated flow and identify points 
of separation and transition. Examples of these smoke photographs are shown In 
figures 6 and 7. Interestingly enough, even though several good photographs 
were taken at the same ambient and flow conditions, slightly different transi- 
tion locations were obtained. High speed movies taken with the Wollensak 
Fastex camera confirmed this unsteady phenomenon in that the laminar portion of 
the separated flow was observed to grow and shrink in length during small time 
intervals. The still photographs show that the separation point is always 
within 1% of the top of the cylindrical leading edge. The angle at which the 
flow departs from the surface Is between 9° and 22°, as measured from the 
photographs. 


Oil Surface Flow 

Before starting up the wind tunnel, the angle of attack and flap angle 
were chosen and a line of oil (approximately 1 cm x 10 cm) was placed on the 
center of the airfoil in a chord-wise direction. Once the tunnel was operating, 
the resulting airflow forced the oil to flow along the surface. Within one 
minute, the separation line along the leading edge was clearly visible. Behind 
this line were three different sections of oil: (1) oil flowing towards the 

leading edge, presumably part of the reverse flow inside the separation bubble; 
(2) a dark region where the turbulent flow Is reattaching to the surface, 
forcing the oil to ooze In both directions; (3) oil flowing towards the 
trailing edge via the turbulent attached boundary layer (fig. 8a). In all 
these experiments, the point chosen as the reattachment point was the solid 
line separating the first light region and the dark region, because this line 
was clear, definite and always visible. Actually, it was reasoned that the 
real reattachraent line is a few percent chord further downstream but the line 
chosen in this work is easily measured and is visibly affected by changes in 
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Rnynoldfi number nnU nnftlo of attack. The error in extremely nmall in canea of 
pofiitive and ?.ero anglen of attack, an gravity will pull the oil an far towardn 
the trailing edge (hence roattachment point) an ponnible. The error becomen a 
bit more aignificant, on the order of .3-4% chord, in the canea of negative 
anp.len of attack, an gravity forcen the oil towardn the leading edge, hence 
driving it nllghtly further away from the. actual reattachment point. 


Plotn of the chord Reynolda number veraua percent chord wore obtained for 
each given angle of att.ack and flap angle, e.g, figuren 9 and 10. Separation 
wan virtually unaffected by inereanlng the Reynolda number, while the tranal- 
tlon and roattachment pol'.ito moved toward the leading edge in an exponential 
manner, aa dlscuaaed by Huang and Hannon (ref. 2). As angle of attack 
increased, the region of turbulent mixing grow larger. At the higheat Reynolda 
number used in this won*. (470,000), the position of reattachment was always 
within 14 to 18% chord, rogurdlesa of angles of attack and flap angle. At 6° 
angle of attack and 30** flap, reattachment was shown at 40% chord. Figure 9 
shows at Rcc" 140,000 this point Is moved to 25% chord, and figure 10 placed 
reattachment Just beyond the 20% mark. Consequently, the slope of these curves 
is highly dependent on angle of attack. 

While in the process of experimenting to obtain these reattachment points, 
an unexpected phenomenon occurred. With the angle of attack set at -4°, 0® 
flap angle and a chord Reynolds number equal to 469,000, the separation line no 
longer appeared to be merely straight, rather a scalloped wave-like pattern was 
formed. Figure 8 demonstrates the formation process of this pattern. Three 
distinct regions were observed in the chord-wise direction of this oil flow. 
From 10 to 12% chord a thick build-up of fluid occurred in an almost regular 
sinusoidal pattern along the upper surface of the model. Following this con- 
figuration was a region extending to the 14% chord-line containing a thinner, 
more dispersed layer of oil, which accounts for most of the pattern seen in 
figure 8d. The third zone exhibits a very small amount of fluid, again 
stretching about 2% of the chord, which assumes a shape similar to trailing 
vortices. This basic shape formed at the higher Reynolds numbers (above 

400.000) for every angle of attack and flap angle tested. The scalloped line 
is a manifestation of the three-dimensionality of the transitional separated 
flow. The formation process was interesting in Itself. The oil film was 
placed on the model and the tunnel was turned on to idling speed (Pec = 

150.000) and a very steady separation line formed at angle of attack -4°. An 
Increase in Rec to about 300,000 saw the separation line become somewhat 
unsteady, as it oscillated back and forth. Further increases in speed (Rec ® 

370.000) resulted in the ceasing of the oscillatory motion in conjunction with 
the first appearance of the scalloped pattern. In this Reynolds number range, 
the line was very sensitive to velocity changes. An increase in tunnel speed 
resulted in the rebirth of the rapid oscillations for approximately 30 seconds 
before damping out; once again the oil assumed the wavy-like form. At maximum 
speed (Rec “ 450,000) the oil, in the scalloped pattern, was Impervious to 
incremental velocity changes. 

With the tunnel at peak velocity, it was desired to observe any changes 
evident in the surface oil pattern with respect to angle of attack Increments. 
The angle was raised from -4° to -2° to 0° without any real alteration encount- 
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ered. However, once positive angles were Introduced, a definite change in the 
wave pattern was noted, causing the frequency of the nodal points to Increase 
in the sinusoidal wave. As before, just prior to the alteration of shape, an 
oscillatory motion was observed, eventually ceasing within 30 seconds, A simi- 
lar physical situation was noted by Mattingly (ref, 3) in 1962, His experiment 
Involved the flow about a circular cylinder placed in a water tank. What he 
observed was an undulating separation line. He noted that the source of the 
separation line oscillation is related to the three-dimensional motion of the 
wake fluid directly behind the separation line, Mattingly calculated the fre- 
quency of these oscillations to be that of the prevailing shedding frequency. 
The frequency of the oscillatory motions observed in the oil flow line on the 
CLE-CT model was not determined. By visual observation, however, as in 
Mattingly's work, an increase in Reynolds number caused a corresponding 
increase in this frequency. These oscillations eventually damped out, however, 
and the oil f:\ow line formed a steady, wavy pattern. Further experiments using 
hot-wire anemometry are planned to study this phenomenon in greater detail. 


CONCLUDING REMARKS 


Results from this study indicate that visual techniques can be used to 
Increase our understanding of separated flow phenomenon. As adequate analytical 
methods are presently not available to predict the location and detailed struc- 
ture of transition, the insight gained from these visual experiments is of 
significant importance , 

Smoke flow photographs clearly visualized the free shear layer following 
separation from the leading edge. Also, due to the favorable thickness to 
length ratio of the bubble, the formation of Tollmein-Schlichting type waves in 
the laminar shear layer is quite distinct, as is the breakdown into a turbulent 
motion. These photographs will aid in the placement of the hot-wire probes in 
future investlgaticins. The effect, if any, of the probe on the flow in and 
around the separation bubble can be easily determined. The unsteadiness of the 
transition process was clearly observed in the high speed movies* 

An oil mixture placed along the upper surface of the model revealed the 
reattachment locations of the turbulent shear layer to the surface. These 
points, along with those taken from the smoke flow photographs, were plotted 
with the static pressure distributions. These plots indicate that the chosen 
parameters of chord Reynolds number, angle of attack and flap angle definitely 
affect the size and characteristics of leading edge separation bubbles. The 
surface oil flow produced a spanwlje scalloped pattern at the higher Reynolds 
numbers. This oil flow pattern is a manifestation of the three-dimensional 
nature of the transition and separated flows. Hot-wire anemometer studies in 
progress will provide additional insight into this important but extremely 
complex flow problem. 
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Figure r..** Low turbulence subsonic smoke wind tunnel. All dimensions 

are In millimeters. 
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4 30* FLAP deflection 

Figure 2.- The two-dlmenclonal cylindrical- leading-edge constant-thickness 

(CLE-CT) model. 
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Figure 3.- Turbulence Intensity at - 37 cm plotted against tunnel velocity 
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Figure 4.- Turbulence intensity at « 76 cm plotted against 

tunnel velocity. 
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Figure 5.- Static-pressure coefficient along the chord of the CLE-CT model 

for o • 4°, 6 » -30®, Re « 1.59 x 10^ and for a • 6®, 6 ■ 0®, 

Re - 1.54 X 105. ' 

c 



Figure 6.- Smoke-flow photograph at angle of attack of 4°, flap deflection 

of -30®, and Re^ - 1.59 x 105. 













Figure 9.- Measured separation, beginning transition, a^ reattach^nt 
locations against chord Reynolds number for the CLE-CT model for 
a ■ 4® and 6 • -30®. 



Figure 10.- Measured separation, beginning transition, rad reattachment 
locations against chord Reynolds number for the CLE-CT model for 


a -40 5 • 0®. 
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INSTRUMENTATION, TECHNIQUES AND DATA REDUCTION 

ASSOCIATED WITH AIRFOIL TESTING PROGRAMS 

AT WICHITA'STATE UNIVERSITY 
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Wichita State University 

H.C. Seetharam 

Boeing Commercial Airplane Company 


SUMMARY 


Two-dimensional airfoil testing has been conducted at the Wichita State 
University Beech Wind Tunnel for a number of years. The Instrumentation de- 
veloped and adapted during this period of testing for determination of flow 
fields along with the traversing mechanisms for these probes are discussed. 

In addition, some of the techniques used to account for Interference effects 
associated with the apparatus used for this two-dimensional testing are pre- 
sented. The application of a mini-computer to the data reduction and presen- 
tation Is discussed. 


INTRODUCTION 


During the wind tunnel testing program of two-dimensional airfoils at the 
Walter H. Beech Wind Tunnel at the Wichita State University during the last 
three years, certain Instrumentation has been developed or adapted. These 
Instruments Include pressure probes for obtaining measurements near the sur- 
face of an airfoil and pressure probes for the determination of the velocity 
fields In reverse flow regions. Traversing mechanisms were also developed for 
positioning these probes. These probes and traversing mechanisms are described 
and samples of the data obtained are shown In this paper. In addition to the 
pressure probes, commercially available hot-film probes were adapted for mea- 
surement of velocity In the reverse flow regions. The hot-film probe came In 
two versions, a single-element probe and a split-film probe. These probes and 
their use are discussed and data obtained with the hot-film probes are shown. 

Use of the probes and traversing mechanism In some regions and under some 
conditions Introduced Interference effects on the flow field about the model 
specimen. These interference effects arc discussed along with modifications 
and techniques used to evaluate, minimize and correct for these Interference ' 
effects. One of the effects discussed Is that on the separation of the bound- 
ary layer and consequent changes to the pressure distribution. 
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Large amounts of data generated during wind tunnel programs can become un- 
manageable. A mini-computer has been Implemented In the data acquisition and 
processing system at the Wichita State University wind tunnel. Use of the 
mini -computer system for data acquisition, reduction, analysis and display Is 
described. 

Symbol definitions are given In an appendix. 


INSTRUMENTATION 


measurements are those required to determine and describe the 
flow field about the airfoil. For an Incompressible flow, if one determines 
pe velocity field both In magnitude and direction then one has all of the In- 
formation necessary to understand the flow mechanism. The velocity field 
about the airfoil at the pre-stall through post-stall range of angles of at- 
tack of Interest varies quite widely In character. The flow can be steady or 
unsteady with varying degrees of turbulence In different regions about the alr- 
foll and also normal or reversed In direction. The regions of Interest about 

field, the far-fleld flow, and the 
wake difficult to get a single instrument that can operate 

under all of the conditions found In these different regions. A number of dlf- 
ferent Instruments are therefore used. These Instruments are discussed in the 
following sections, under two categories: 1) pressure probes, and 2) hot-film 
probes. In addition, the traversing mechanisms used to position these probes 
are also discussed. 


Pressure Probes 

The basic pressure probe used Is the five-tube probe shown In figure 1. 
This probe consists of a total head tube and four static taps spaced about the 
circumference of the probe. The description, calibration and use of this 
probe have been documented In references 1 and 2. The probe Is used to deter- 
mine the magnitude and direction of the velocity In the far-field region over 

< airfoil. In regions where the angularity 

of the local velocity with respect to the probe axis exceeds ±40®, the probe 
cannot sense the direction or magnitude. Because of the physical dimensions 
* cannot be used for near-wall (<0.2 mm) measurements. In order 

to obtain velocity measurements closer to the surface of the airfoil the flat- 
tube probe shown In figure 2 was developed (reference 1). This Is a total 
head tube whose prime component Is a flattened hypodermic tubing. This probe 
is capable of measurement as near as 0.15 mm from the local surface. The 
velocity magnitude Is determined using the static pressure tap on the surface 
of the body for static pressure quantity. The flat-tube probe cannot sense 
the velocity angularity or direction, and like the other pressure probes. Is 
Insensitive to velocity fluctuations. 
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In an attempt to determine the separation and reattachment points over the 
airfoil In the near-wall region, the cylindrical -tube probe of figure 3 was 
developed (reference 1). This probe consisted of two hypodermic tubes, s1de» 
by-slfl'e, sealed at the extreme end with a small total pressure head opening 
drilled Into each tube near Its sealed end. During use, the tube was brought 
close to the surface of the airfoil and the flow direction and total pressure 
were obtained from the higher two pressures. The cylindrical -tube probe could 
not sense the flow angularity. The velocity magnitude was determined using 
the airfoil surface static pressure reading obtained simultaneously. This cyl- 
indrical-tube probe was found to be rather Inaccurate for determining the sep- 
aration and reattachment points, because of sensitivity to si dewash velocity 
near separation. 

Another configuration developed for near-wall velocity measurements and 
also capable of determining flow reversal was the four- tube probe shown In 
figure 4. This probe has forward- and aft-facing flattened total tubes for 
near-wall flow direction and total pressure measurement. The pi tot- type tubes 
are less sensitive to sidewash than the cylindrical -tube probe, and are cap- 
able of iTsasurement nearer to the surface (down to 0.34 mm) than the cylin- 
drical probe. The two static pressure tubes added to this probe give It the 
capability of direct static and total pressure measurement, and therefore ve- 
locity magnitude. The probe Is not capable of determining flow angularity. 
Since flow direction near a wall Is usually nearly parallel to the wall, the 
lack of flow angularity measurement capability Is not a severe limitation. 

This probe was adapted as the secondary pressure device used In the flow mea- 
surement work supplementing the data from the five-tube probe. The four-tube 
probe Is described In more detail In reference 3. 

A comparison of typical velocity profiles measured over the airfoil using 
the cylindrical-tube, the five-tube, and the flat-tube probes are shown In fig- 
ure 5. The differences In velocity as measured by the different probes are 
rather small. The five- tube probe data Is probably the most accurate since It 
has been calibrated extensively as mentioned previously. 

Another example of the data obtained for determining the flow field about 
the airfoil Is the velocity field data shown In figure 6. The vectors are to 
scale and Indicate both the magnitude and the direction of the flow over the 
airfoil and In the wake. As can be seen on the figure, the five-tube data does 
not supply any Information In the reverse flow regions. The four- tube data 
shown by the dashed lines Is used to supplement the five- tube data although the 
velocity angularity In this reverse flow regli-n Is not determined. 

Another pressure probe developed and used was the boundary layer mouse 
shown In figure 7 and discussed more fully In reference 4. This probe or rake 
consists of a number of total head tubes stacked closely and attached tangent 
to the surface of the airfoil The total pressures determined were used with 
the surface static pressures to obtain the velocity magnitudes near the wall. 
Velocity magnitudes In the reversed flow regions cannot be determined with 
this rake. A typical total pressure measurement survey over the flap and the 
flap cove of an airfoil Is shown In figure 8. 
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Hot"Film Probes 

ciated^SltrSl“sianM research program, ^ bot-i’ilm sensor and asso- 

JJJo* « conditioning equipment were ofctaineo (reference 4) Th<» 

^’9t-f11m sensor shown in f1^ 9 Ih^ sen- 
"Ire. .060 n>n In dialler! JSunind tet^S; 

^ pidtinuTj. The film Is heated by proper circuitry and the rate 
of heat transfer depends on the velocity of the flow past the wire The hit 
film probe Is sensitive to unsteady flow velocities and Is caiih»*At 4 aH tA maa” 
sure the .mgnitude of the veloclty^^ HweleJ! ? «nSS? «n“ W^dlJeairof 
the f ow if reversed. One can speculate from the datrobta1%d whether 

bwome^zero®^ velocity would Intermittently 

Recently, split-film sensors have become available which eliminate the 
f ow direction ambiguity associated wnh singled m sin?orr( 5 ? 

spllt-ffln, sinsor dhln flliTStaMnum 
prob?2hlJh^ir */* portion of the cyHndrlcel sensor of this 

jiiowphrdieri^:?i;roVThVVioi''^?;cifs^ 

«JettertrnoiNrin1he'7i?fir*!:?t'^^ 

film sensors are described more thoioSghJrin referel!^^ 

Figure 11 shows how the output of the split-film [a voltaael and thp fmir 
tube pressure probe [a pressure differential] vary with the vefocity!^ ^ 

oscnwiL-fn?3?w^°'^ 3^? split-film eaulpment appear: on an 

versa! if flow turbulence and percentage of time of re- 

fiin, ^ sample of a calibrated velocity trace obtained bv the snlit- 

13. This temporal data c^be anal y lid ti 

a? Ssu fi? {hlllJaijsir “ «•*««" 

a ?r?i®"^ presenting the data Is that shown In figure 14. 

thi^S distinct regions about 

the airfoil and In Its wake where certain characteristics of the flow exist. 


TRAVERSING MECHANISMS 


ina !? P''®''1ous sectlon Were mounted onto a travers- 

schematically In figure 15. The probe positioning could be 
^^® wind tunnel both In the vertical and axial direc- 
tion. While trying to repeat previous data of a high angle of attack ca«tP if 

the prob^d1ffiiid^^iai^f^®‘^^l^®^ pressures obtained on the model without ^ 
the probe differed significantly from that when the probe was present. This 
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e. ;s.T«^ E E i 

n??,‘2;;sJ5%Sru ^s^^<M;;l!."'At ?S.*W »f .ttac. the pressure 

data shows insignificant interference. 

„ .,:: rr^!.=sr«:;T»'y&;:.K.3 

tional problems that f ‘’‘‘““S *''! oldInS the S?lbe rotates within 

retSsln,!^chrn^nS' thetn;^^^ .rp^rUEl lEl' 

S:rs5'r?:EE! ’■TSEnS fro’’EEiE S Jth the s^^ and probe Is ^u^ted ro the 


INTERFERENCE ON INSTRUMENTATION 


,n the course Of the ^ V."ZlTon ‘t”he1SE™5»nEE'onal 

t:^?;oE'^w^ir in ^rS nElo rEEe.rdlu ??Em previous tests It was fo.^d that 
section walls, in try g p deflected the data was somewhat in 

s.'rixw.'S'i.K s t:=; :!•.» 

velocity field about the airfoil. 

statiE“;EEno^riEr?EEJr„raErtEeEEtEroT^^^^^^^^ 

on figure 19 which depends on the square of the lift coeffic 
recti on is discussed more fully in reference 6. 
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DATA ACQUISITION, REDUCTION AND DISPLAY SYSTEM 


The components of the force data reduction and display system are shown 
In figures 20 and 21. The system Is described fully In reference 7. A mini- 
computer serves as the central element of this system. 

Figure 21 shows the complete elements of the configuration with the mini- 
computer at the center. Various input and output devices are associated with 
the mini -computer. This has allowed more refined treatment of the data and 
Improved the efficiency of the data acquisition and reduction. The engineer 
now has results available almost Instantaneously and can determine whether the 
data point Is satisfactory or whether It should be re-run. 

Data can also be plotted on line using the digital plotter. Typical re- 
sults are shown In figures 22 and 23. The first of these shows a sample of 
velocity field plot. The airfoil and velocity vectors are plotted about 6 
seconds after the data points had been obtained. A sample pressure distribu- 
tion plot of an airfoil with flaps Is shown In figure 23. 

Computer programs are being formulated and written to obtain statistical 
Information from the histograms obtained by the split-film probes. This should 
allow more Information to be obtained from the data presently generated during 
the wind tunnel tests. 
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APPENDIX 


SYMBOLS 


Values are given In both SI and U.S. Customary Units. The measurements 
and calculations were made in the U.S. Customary Units. 


C 

C, 



h 


Ps 

Pt 

Pco 

Q 

RM 


I 


X 

z 

r 


Airfoil chord 

Airfoil section lift coefficient, section lift 

PS ^ Pc» 

static pressure coefficient, — q — 

Pf “ P« 


Total pressure coefficient, — q- 


Wind tunnel test section height 

Location of wall static pressure tap from airfoil quarter chord 

Local static pressure 

Local total pressure 

Free stream static pressure 

Free stream dynamic pressure 

Reynolds number based on wing chord and free stream conditions 
, «;tream velocity 

'nsionalized component of local velocity in the free stream 

1 

Axial component of velocity at wall static pressure tap 
Resultant velocity at wall static pressure tap 
Streamwise coordinate 
Vertical coordinate 
Airfoil circulation 
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Figure 1.- Five- tube probe. 
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Figure 2.- Flat-tube probe. 
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Figure 3.- Cylindrical- tube probe. 
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Figure 4.- Four-tube probe. 
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Figure 5.- Typical velocity profiles over an airfoil. 



Figure 6.- Velocity field data 
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Figure 7.- Botmdary- layer mouse. 



Figure 8.- Mouse total-pressvire measurements. 
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Figure 9.“ Single hot-film sensor. 
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Figure 10.- Split-film sensor. 
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Tunnel Floor 


Figure 15.- Probe traversing mechanism before modification. 
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Figure 16.- Probe traversing mechanism after modification 
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Figure 17.- Surface static-pressure interference. 



Figure 18.- Probe traversing mechanism. 
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Figure 19*- Dynamic-pressure correction for wall static-pressure error 
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Figure 20*- Force data-reductlon system* 








rnmMiw 


Figure 21 ♦- Elements of the configure ! 



Figure 22,- Sample velocity field plot 
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Danpy R. Hoad* 

NASA Langley Research Center 
SUMMARY 

A laser velocimeter (LV) was installed in the Langley V/STOL tunnel to 
measure the velocity field about a wing with a NACA 0012 airfoil section. These 
i^dsurements were con^jared at low angle of attack with a two-dimensional viscous, 
flow prediction program. The velocity field over the wing in a fully stalled 
condition was also measured by the LV. The unique ability of the LV to measxnre 
absolute flovr magnitude and direction without prior knowledge of general flow 
direction was demonstrated in the complex separated reverse flows near the 
wing \Q>per surface at the high angle of attack. 

The general characteristics of the flow field over the stalled wing were 
substantiated by a vapor screen flow visualization technique. 

INTRODUCTION 

The laser velocimeter (LV) is a relatively new device to measxire fluid 
velocities. It is unique with its capabilities of obtaining these measurements 
in flow conditions where conventional devices either cannot obtain measurements 
or woxild seriotisly influence the measurement dcae to the presence of the device, 
^is fact indicates that it is a viable tool to determine velocity field con- 
ditions about an airfoil, and is particularly important near the surface where 
a conventional probe’s presence wotild induce velocity components, thereby 
biasing the measured data. 

The performance of an airfoil is typically judged by body force measure- 
ments or surface pressure measurements. In fact, two-dimensional prediction 
techniques are designed to estimate the local surface pressures as a prime 
objective. One such prediccion technique is described in reference 1. It has 
been shown to accurately predict the surface pressures on a single element 
airfoil at low angles of attack. 

The performance characteristics of the NACA 0012 airfoil are well-defined. 

It is not the intent of this paper to present anything new about the airfoil. 
Because of its application in the aerospace industries, particularly on heli- 
copters, it was chosen as the baseline model to demonstrate the application oi 
the LV to airfoil flow— field research. 


« 
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SYMBOLS 


Hie axes used for this investigation a.re presented in figure 1. The 
velocity measurement position was referenced to the airfoil chord line, and 
toe velocity measurement magnitude was referenced to the free-stream direction, 
^e for the physical quantities defined in this paper are the Interna- 

tional System of Units. Most quantities were measured in this system; however, 

^me were measured in the U. S. Customary Units and converted by using factors 
given in reference 2. 

c wing chord, O. 30 U 8 m 

U^ local velocity conq>onent, in direction of fi'ee-stream velocity, m/sec 

local total velocity , ^ , m/sec 

U^ free-stream velocity determined from pitot-static probe, m/sec 

Vj. local velocity conq)onent, perpendicular to direction of free-stream 

velocity, m/sec 

^c*^c coordinate axis relative to wing chord 

distance downstream from airfoil leading edge along chord, m 
y^ distance above and perpendicular to wing chord, m 

a wing angle of attadc, deg 

APPARATUS 

^ optics system operating in the backscatter mode was used 

for these tests. A sketch of the optics system is presented in figure 2 and a 
photograph is presented in figure 3. A high-speed burst counter was used to 
measure the period of the hi^-frequency signal contained in toe burst f^^jm a 
particle traversing the sample volume. LV system control, data acquisivion, 
and data reduction were handled by a minicomputer. A complete description of 
the LV optical system, electronics system, and data acquisition and reduction 
is available in reference 3. 

The model used in this investigation was a simple straight wing. It had a 
span of 2.438 m, a chord of 0.3048 m, and a NACA 0012 airfoil section. Velocity 
measurements were made at center span to obtain two-dimensional characteristics. 
The wing was supported by struts from the floor near the tunnel center line 
with no balance measurements taken. The location of the strut mount to the 
wing was chosen as far outboard as structurally feasible to minimize flow dis- 
t’lrbanoe at the wing center line. A photograph of the model with crossing laser 
bf oms is presented in figure 4. 
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This investigation was conducted in the Langley V/STOL tunnel at a nominal 
free-stream Mach nuniber of 0.I5. The Reynolds number based on the wing chord 
was approximately 1 x 10°. Local flow velocities were measured abotzb the 
center line at two geometric angles of attack: (l) at a « 0® in order to 

compare with a two-dimensional theoretical prediction technique at a condition 
for whi<di the flow is two-dimensional and the theoiy is appropriate » and (2) 
at o s 19.40 in order to investigate the flow-field characteristics over the 
separated airfoil at a constant angle of attack. A pitot-static probe was 
mounted 2.5 m below and 1 m ahead of the wing center line to provide accurate 
reference of the free-stream tunnel dynamic pressure.' A hygrometer was used 
to obtain wet-bulb tesperatxires , and total temperature was measured in the 
settling chaniber. Thus , the tunnel air density could be calculated and, with 
dynamic pressure measurements , the tunnel velocity could be accurately calcu- 
lated. 


DISCUSSION 

The velocily .neasurements at eacdhi measxirement location were first reduced 
to histogram form. Ihese data for the wing at a = 0® (relative to txmnel 
geometric center line) along with a description of data reduction techniq\ie, 
hxstogram interpretation, and conplete error analysis can be found in 
reference 3. 

Free-stream velocity measuiements were obtained with the LV with no wing 
in place at the location of the wing center line. These data indicated an 
average upwash angle of 0.6® (relative to txaimel geometric center line). The 
wing was installed with '•■’'e dbord line pareO-lel to the test-section center line; 
therefore, the effective angle of attack was assumed to be 0.6®. 


Prediction Technique 

The extexnal forces generated on a body in a fluid are manifested in the 
velocity distribution of the fliiid about the body. In developing a prediction 
tedmiqw, the calculations at the surface of the bo<ty are verified with con- 
ventional pressure and force measurements • Reference 1 presents an excellent 
comparison with measured surface pressures for this viscous-flow prediction. 
Since the local surface presstires are conputed frcxn predicted local surface 
velocities, it is justifiable to question the validity of the predicted velo- 
cities away from the surface. Th<=’ need to verily these predicted velocities is 
evident. The \ise of conventional probes near the surface raises questicms cibout 
the accuracy of the measurement with interferences caused by the presence of 
the probe. It was determined that the LV was a device capable of measuring 
this flow field without inducing any interference since nothing was present in 
the field but the wing and light beams. 

The theory for this prediction technique (ref. 1) involves an iterative 
procedure which first obtains an inviscid-flow solution for the basic air- 
foil. It computes a b our. d ary-layer solution based on the inviscid-flow solu- 
tion and constructs a modified airfoil by adding the boundary-layer displacement 
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thickness to the original airfoil* It ottains the Invlscid solution for the 
modified airfoil and repeats these steps until appropriate convergence cri- 
teria are satisfied* Hie field point velocities avay from the surface are 
then computed from the basic vorticity distribution along the modified airfoil* 


Experiment-Theory Comparison at a « 0® 

The velocity vectors as measured by the LV for the »fing at a geometric 
angle of attack of 0® are presented in figure 5* Eadi velocity vector (arrow) 
is an average of an ensenible of measurements taken over a short period of time 
at the desired location* This arrow plot indicates the relative location of 
the velocity measurements , magnitude , and angle of the velocity vector* The 
velocity magnitude and direction are indicated by the length and orientation 
of the arrow* The tunnel free-stream magnitude and direction reference are 
provided in the lower left comer of the . figure* The position of the velocity 
measurement is mariced by the tail end of the eurrow* 

The velocity measurements were obtained by positioning the sasple volume 
at a desired chordwise station (x /c) and incrementing the entire optics 
padtage downward along this dtiordwise station. This was acconplished remotely 
and was completely controlled by the miniconputer* Four of these series of 
measxu*ements (scans A, B, C, and D from fig* 5) are presented in more detail in 
figure 6* This figtire presents a compsulson between LV-measured velocities 
and the two-dimensional visco\Js-flow prediction* The compazdson is presented 
with the resultant velocity nondimensionalized by tunnel free— stream velocity 
as a function of the verticeQ. position of the measurement nondimensionalized 
by the wing chord* The scans near the leading edge have velocity gradients 
which are the iaost difficxilt to predict* The free-stream rpwash angle without 
wing or supports was measured at 0*6® at this Math number* Topically, flow 
angularity is affected by a model's presence* It is normally determined by 
model \pright and invi-rbed angle-of-attaCk ranges* Cmpeurison of balance 
data from these two conditions provides the total flow angularity* It is very 
difficult to obtain this type of measurement with discrete velocity measurements 
in the presence of the model* Since no balance meeeurements were obtained on 
this xnvestigation, there is some uncertainty in the effective angle of attack 
of the wing* Predicted velocities were calculated first using the measured 
0*6® tunnel flow angle without the wing* These are presented as dashed lines 
with maximum discrepancies on the order of 6 percent* Calculations were 
repeated with a 1-degree shift in angle of attack to provide an eissessment of 
the effect of uncertainty in this measurement* These calculations are pre- 
sented as solid lines with a = -0*U® and indicated better agreement with theory. 
It is obvious in these comparisons that the precise measurement of these velo- 
cities depended on the precise determination of the effective angle of attack 
of the wing* It is Justified to say, however, that these data provide a 
quantitative and qualitative measvire of the accuracy and applicability of LV 
measurements about a surface submerged in a fluid flow situation. 
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Airfoil at a ® 19.4® 


An arrow plot of the mean velocity held about the airfoil at o « 19.4® 

8 presented in figxire 7. Each velocity vector is presented with the length 
of the ^row indicating magnitude and direction of vector relative to the air- 

from this figure that the airfoil is in a fully separated 
condition from leading edge to trailing edge. The shear layer region between 
the free-stream and the separated turbulent area over the wing is broad, but 
ewily discernible. The velocity field in the separated region indicates 
e existen^ of a large recirculating eddy with reverse flow near the airfoil 
s^face. The velocity fluctuations within the shear layer were large; however, 

velocity fluctuations were smaller. A dashed 
line is provided indicating the approximate location of zero velocity in the 

trailing edge (see fig. 8) , a veiy shaip shear layer 
is evident with low, reverse flow velocities generated near the airfoil imper 
surface md with nearly free-stream velocity from the lower surface. The 
spatial distance across this shear layer is on the order of 0.005 y /c. The 
reverse flow in the wake region above the airfoil is also evident. ° 


^ese relocity measurements were obtained without prior knowledge of the 
dl region of the flow at each measurement point. The LV is unique in this 
capability unlike conventional probes, which require this information to 
reduce ambigtiity primarily caused by support structxire interference. 


Flow Visualization of Airfoil at a ® 19.4® 

To visualize the flow patterns over the airfoil, a thin plane of light 
(approximately 2 cm thick), perpendicular to the wing trailing edge and parallel 
to the free-stream direction , was projected from behind the wing along the 
TOnterspan. The only flow patterns visable were those depicted by smoke 
trarersing this plane. The patterns were recorded on video tape with still 
p ctwes taken later from a television monitor. The reversed recirculating flow 
in the separated region was observed, thus confirming the directions and general 
flow patterns as measured by the laser velocimeter* 


toe of the pictiires taken from the television monitor of the flow is pre- 
sented in figure 9. The television camera which recorded this frame was not 
positioned orthogonal to the wing chord due to dii'ficulties in locating the 

rear* “^he airfoil section evident is 
the tip and not the airfoil section at the location of the va;ior screen. The 
lea^ng edge of the airfoil at the location of the vapor screen is near the 
eading edge of the separated region. The trsdling edge of the eiirfoil at this 
location is at the trailing edge of the bright streak. This streak is the 
reflection of the sheet of light from the wing upper surface. 

The separated region is evident in this figure, although the reverse flow 
just above the airfoil is not. The sharp shear layer at the trailing edge can 
be seen. With better photographic techniques, the vapor screen technique with 
the laser as a powerful light source can be a very useful research tool in 
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assessing the character of the flow field as measured hy the laser velociraeter. 

CONCLUSIONS 

A laser velociraeter was installed In the Langley V/STOL tunnel to measure 
the velocity field about a straight wing with a NACA 0012 airfoil section. The 
wing was installed at zero angle of attack to provide data to con^are with a 
well-accepted two-dimensional viscous-flow prediction program. These results 
provide a qualitative confidence level in the accuracy of the LV measurements. 
The wing was also installed at a fully stalled condition (a = 19. U®) to charac- 
tei*ize the flow field in the separated region. 

The results of the investigation indicated that 

1. !I5ie Ifiser velociraeter is an effective and accurate instrument for 
measiiring the velocity field about a surface. 

2. The precision of the laser velociraeter for the low angle-of-attack 
data in this case depended on the precise determination of the effective angle 
of attack of the wing in the tunnel. 

3. The separated region over the wing in the fully stalled condition was 
well-defined. Hie shear layer between this region and the free stream was 
hroad euid highly turbiilent; however, in the reverse flow region the measure- 
ments indicated relatively lower turbvilent characteristics. 

U. The measvirements in the reversed flow region demonstrated the mique 
capability of the laser velociraeter for measuring velocity magnitude and 
direction without prior knowledge of the flow direction. 

5. The trailing-edge measurements demonstrated the capability of the laser 
velociraeter to raeasure the velocity characteristics across a very sharp velocity 
gradient. 


6. The vapor-screen technique with the laser as a powerful light source 
was deraonst rated to be effective in assessing the character of the laser 
velociraeter measured flow characteristics. 
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Figure 1.- Sketch of axis system used Including dlrectlcn of 
velocity components computed. 




Figure 2.- Schematic of laser veloclmeter optics. 
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SCANS 



Figure 5.- Velocity vectors computed from measurements over wing. 
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Figure 6.- Comparison of laser velocimeter flow-field velocity 
measurements with a two-dimensional viscous-flow prediction 
program. 
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Figure 7.- Velocity measurements over the stalled wing (a 
free-stream Mach numibeT • 0.13. 



Figure 8.- Velocity measurements at trailing edge 
of wing at a « 19. 



4 



Figure 9.- Flow patterns over wing depicted by vapor-screen 
flow-visualization technique. 
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APPLICATIONS OF DYNAMIC MEASUREMENT TECHNIQUES 

TO WIND-TUNNEL TESTING 

P. K. Owen* 

Consultant 


SUMMARY 


The uses of l,aser doppler velocimeter» hot wire and surface hot film tech- 
niques in the study of turbulent flows are described and data obtained in com- 
pressible flows are discussed. Applications are Illustrated with measurements 
of wind-tunnel freestream turbulence characteristics and with data obtained in 
transitional, turbulent and separated shear flows. A new method which has been 
developed for the study of time dependent and unsteady turbulent flows is also 
presented. 


INTRODUCTION 


Although the title of my paper suggests a broad subject area, in the 
interest of conciseness, I have decided to restrict my conments to the areas 
in which I have had direct experience and to stress measurements which I feel 
are particularly relevant to laminar flow and unsteady aerodynamic research. 
Thus, the paper will be restricted to applications of hot wire and hot film 
anemometer and dynamic pressure measurements and to laser velocimeter measure- 
ments of unsteady flowfields. 

At the present time little is known about the Influence of freestream 
flow turbulence on steady and dynamic measurements on models in wind tunnels 
at transonic speeds. Indeed, few measurements have been made of the charac- 
teristics of freestream unsteadiness in transonic wind tunnels. The result is 
that information on velocity, pressure and temperature fluctuations, their amp- 
litude, frequency, phase relation and space-time correlation, is lacking. This 
information is needed if we are to accurately assess the relationship between 
wind-tunnel and flight behaviour. 

Perhaps the major open question is the Influence of freestream disturban- 
ces on model boundary- layer transition. Recent developments in boundary-layer 
transition research, particularly those of the NASA Transition Study Group, 
have stressed the dominant role freestream fluctuations have on model boundary- 
layer stability at transonic and supersonic speeds. Not only do the external 
fluctuation amplitudes dominate transition but their spectral characteristics 
are particularly significant. This importance of the spectrum of the external 
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disturbances was first demonstrated by Spangler and Wells at subsonic speeds. 
They found a much higher transition Reynolds number on a flat plate than did 
Schubauer and Skramstad even though the freestream disturbance amplitudes were 
similar in the two experiments. However, the spectra of the freestream dis- 
turbances were radically different. 

The Influence of flow unsteadiness on shock-wave— boundary- layer Inter- 
action Is also determined by the connection between boundary— layer development 
and transition location. For example, freestream turbulence can determine 
whether the shock causes laminar boundary- layer separation with or without 
transition or whether the interaction Is with a turbulent boundary layer. 
Freestream flow unsteadiness may also Induce shock-wave oscillations. 

Experiments on flat plates with turbulent boundary layers at zero 
pressure gradient have shown the Influence of freestream turbulence on boundary- 
layer development. An increase in freestream turbulence level leads to an in- 
crease in skin friction coefficient, a fuller velocity profile and a thicker 
boundary layer. Measurements show that a small Increase in free-stream turbu- 
lence has the same effect on the shape of the boundary- layer velocity profile 
as a fractional Increase in Reynolds number roughly 60 times as great. Calcu- 
lations suggest that similar effects occur in flows with moderate pressure 
gradients, with an increase in turbulence level delaying separation onset. 

However, the flow around the wing of a wind-tunnel model goes through 
regions of high acceleration and deceleration which distort the turbulence, 
the further from freestream conditions the more anisotropic it becomes. It 
is possible that the interaction between the distorted turbulence and the 
shear layer could well be quite different from flat plate observations. In 
the case of dynamic measurements of wing buffeting, turbulent pressure fluc- 
tuations can mask the data and can, therefore, be difficult to detect. There 
are also indications that buffet boundaries change with changing turbulence 
level, presumably due to a movement of the mean shock position. 

It should be borne in mind that the object of wind-tunnel testing is to 
simulate atmospheric flight so that the key question remains whether the 
levels and scales of turbulence in the atmosphere have any significant effect 
on boundary- layer structure. Since it is to be expected that the most ener- 
getic scales will be many times the boundary- layer thickness even on large 
transport aircraft, the. primary influence on boundary layers in flight will 
probably be on the large-scale unsteadiness of the flow. 

Thus, it is essential that the freestream characteristics of wind tunnels 
used in advanced aerodynamic testing be thoroughly documented. In this way, 
disturbance levels and scales can be assessed in relationship to those 
thought likely to be encountered in flight. 

Symbols are defined in an appendix. 
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DISCUSSION 


Facility Disturbance Characteristics 

Since there is very little information available on the structure of the 
freestream turbulence, the scales involved, their magnitude and how they vary 
from tunnel to tunnel, an investigation was undertaken (ref. 1) in which 
transition data and hot wire turbulence and pressure fluctuation levels were 
measured in two wind-tunnel facilities namely the NASA Ames 3.5 ft. and the 
Langley 18 in. Variable Density Wind Tunnels. 

Fig. 1 shows rms values of the mass flow and total temperature fluctua- 
tions in the Ames 3.5 ft. wind tiinuel calculated assvttuing a correlation co- 
of -1.0. Due to the scatter in the data, no trends in turbulence 
level with operating pressure could be established. However, th« mean values 
indicated in Pig. 1, i.e., mass-flow and total-temperature fluctuations of 
2.65 and 0.83 per cent, respectively, should be representative of the free- 
stream turbulence levels over the unit Reynolds range. 

These data show the types of measurements that can be o tained in super- 
sonic flow. Although data interpretation is a little more c*> olex, similar 
techniques can be applied to transonic flows. However, it must be borne in 
mind that rms intensities give little if any indication of the turbulent 
scales involved. This information can be obtained from power spectra gn d 
space-time correlation measurements of the turbulent fluctuations. The 
pressure disturbance spectra, presented in Fig. 2, show that, although most 
of the energy is concentrated at low frequencies, the spectra levels are 
quite different at high frequencies, reflecting expected differences in fluc- 
tuation scale due to wind-tunnel size. 

Some interesting features of the freestream disturbances have been deter- 
mined from two-wire, space-time coarrelatlon measurements. Streamwlse distur- 
bance convection velocities were measured in both facilities and found to be 
independent of scale and equal to 70% of the freestream velocity. This 
result is in good agreement with an extrapolation of Laufer's lower Mach 
number data. Pig. 3 shows the variation of the optimum spatial correlation 
functions in the streamwlse and lateral directions as measured on the tunnel 
centerline in the Ames facility. These results, which indicate that the 
disturbance length scales (calculated for Rjjjj « 1/e) are several times their 
width, are consistent with the concept of radiated sound from the sidewall 
turbulent boundary lay^<r; i.e., the indicated ratio of streamwlse scale to 
lateral length scale of approximately 3.0 agrees with that predicted assuming 
wall boundary-layer source propagation angles originating upstream of the test 
section where the source Mach number was approximately 3.0. These streamwlse 
fluctuation scales correspond to overall wind-tunnel wall boundary-layer 
source lifetimes of several boundary- layer thicknesses. 

Recent cross-correlation measurements have also been made of the test 
section and diffuser static pressure fluctuations in the Ames 12 ft. wind 
tur ,el. Since the probe separation was sufficient (approximately 8 metres) 
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to make the correlation of the background turbulence negligibly small, corre- 
lations of the acoustic modes could be directly measured. With the diffuser 
probe output delayed, it is apparent that there is a coherent acoustic mode 
which propaptes upstream. The propagation speed is determined to be approxi- 
mately equal to the speed of sound minus the wind-tunnel freestream velocity. 


Boundary-Layer Transition 

Apart from the complex coupling of the usually unknown freestream turbu- 
lence spectra and model boundary-layer transition a major source of scatter 
in transition data can be attributed to inconsistent choice of transition 
point IndlMted by different techniques, mostly locating positions near the 
end of transition, which generally have a strong Mach number and unit Reynolds 
number dependence. A more complete picture of transition dependence on these 
par^eters can be obtained from experiments in which the positions of the 
beginning and end of transition are accurately determined. It is of interest 
to note that transition data reported for supersonic and hypersonic flows are 
not generally based on observations of turbulent spots but rather some macro- 
scopic quantity such as skin friction, heat transfer, or surface pitot pressure, 
whose departure from laminar values can be detected only when the turbulent 
Intermlttency is appreciably greater than zero. 

BeMuse of the Importance of correctly determining the onset and extent 
of boxmdary-layer transition, techniques which do not disturb the flowfleld 
and which respond to microscopic changes must be used. Such a technique, 
i.e., one that detects turbulent bursts at the model surface, is the surface 
hot film gauge. An example of the variation of the rms thin film gauge output 
over a range of unit Reynolds numbers on a cone model in the Ames 3.5 ft. wind 
tunnel is shown in Fig. 4’. The curve clearly shows a rise from the laminar 
to the turbulent level, with an intermediate peak. These curves enable three 
distinct points in the transition region to be accurately and consistently 
determined: namely, the onset of transition, defined as the point where the 

rms signal begins to Increase from its laminar value (this onset of Intermit- 
t^cy can be clearly seen on the oscilloscope traces) ; the peak rms signal, ' 
which coincides with the point of maximum turbulent burst frequency (ref. 2) ♦ 
and the end of transition. Examples of the characteristics of the film vol-** 
tage fluctuations through the transition region are also shown on Pig. 4. 

Measurements of this type enable the effects of Mach number and unit 
Reynolds number on the beginning and length of transition to be established 
more precisely than with previously used methods. Examples are given in 
refs. 1 and 2. 


Turbulence Structure Measurements 

As in the freestream, turbulence Intensity measurements of the mode 
fluctuations can also be obtained across compressible turbulent shear layers. 
Once again, however, this time from the turbulence modelling viewpoint 
information on the turbulence scales and lifetimes are of crucial importance. 
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since turbulent flows vary not only In tine but also In space, their 
Investigation nust involve an examination of both the spatial and temporal 
statistical structure. Space-time correlations can make a contribution to 
this study since they give evidence of the heredity and structure of turbulence, 
as well as values of the convection velocities of the vortlclty and entropy 
modes compared with the average mass transport velocities. 

Examples of both auto and space-time correlations In a compressible tur- 
bulent boundary layer (ref. 3) are given in Figs. 5 and 6. These data were 
obtained on a cone-oglve-cy Under model in the Ames 3.5 ft. wind tunnel. 

Fig. 5 shows the autocorrelation of the fluctuating signals on the 
cylindrical portion of the model 176 cm from the cone apex, at two positions 
in the boundary layer- and in the far field. It can be seen that there is a 
marked variation of energy distribution with frequency across the boundary 
layer and that, as expected, the far field contains proportionately much less 
energy' In the high wave number range than the wall region. 

The results of a series of filtered (4 kHz) cross-correlation measurements 
at several separations in the boundary layer are shown in Fig. 6. It can be 
seen that each cross-correlation curve reaches a maximum at a non-zero value 
of the time delay, clearly indicating the presence of convection. The ampli- 
tude of this maximum Is a function of the wire separation distance. A convec- 
tion velocity of these disturbances may be determined from the time delay at 
which the maximum of a particular cross-correlation occurs. 

The peaks of the cross-correlations obtained for various values of wire 
separation distance represent the autocorrelation in a reference frame moving 
with the disturbances. They are, therefore, a measure of the lifetime of the 
disturbance pattern as it is swept along with the mean flow. The long turbu- 
lence lifetimes which can be inferred from these space-time correlation measure- 
ments (refs. 3 and 4) illustrate a major objection to turbulence models based 
on local flow conditions. It cannot be assumed that turbulence is uniquely 
related to local conditions, and flow history must be considered, especially 
when attempting to calculate non-equilibrium flows. 


Turbulent Boundary-Layer Separation 

Turbulent separated flows occur in many types of engineering configura- 
tions. They may be unintentional features of some classes of equipment or 
they may be deliberately introduced. But in all cases such flows can have a 
significant effect on engineering performance. Furthermore, additional com- 
plications are added by the unsteady aspects of the turbulent boundary-layer 
separation and reattachment processes. Despite the fact that these flows 
have been extensively studied, detailed Information regarding the unsteady 
nature of turbulent separation is practically nonexistent for high-speed 
compressible flows. Conventional "time averaged" measurements such as surface 
pressure, skin friction, heat, transfer, and pitot pressure surveys cannot 
supply this Information. However, once again thin platinum films mounted 
flush with the model surface provide basic information on the significant 


575 


wltSw separation. The fluctuating 

abwJ Jhe^flLi!*^ ^ ^ provide neasurements related to the flow character 

at. variations of the ms thin film voltage fluctuations throush 

a shock-wave—boTOdary-layer Interaction region (ref. 5) are shown In Pig. 7 
Mso Indicated a^e the measured pressure distributions for the two cases. * 

flow (shock-wave generator wedge angle of 7 5®) 

detailed pltL ^ 

These ma Surveys, surface skin friction and surface oil flow data were obtained 
TOese m^n measurements Indicated attached flow for the 7.5® wedge angle aST* 
substantial region of separated flow for the 15 ® wedge angle. (The realon of 
^on**hsi shear, as determined from a floating element sH^ frlc- 

indicated on the figure.) The thin film results show a marked 
"f® ^®^een the attached and separated flows. Normalized power spectra 
*? th. TMlon «d .*t« 

r»t«l raiion ““ *’* ‘I** «»rgy lacrM.. in th. .epa- 

? is confined to a narrow band around 15 kHz while the Increased 

*V4^**® pressure rise after reattachment is broad band. Power spectra 
194 cm^frL^Jh* enerp peak were obtained at measuring stations between 183 and 
194 cm from the model tip. Similar measurements were also obtained for the 

an«r« puk F**** « aaootli ponar apectra irtth no 

1.^5" ® energy peak was not evident In these cases, it Is 

of^JhJ^^® associated with turbulent separation unsteadiness. The scale 

peak*^fre^eJci T convection velocities and the measured 

p k frequency, la of the order of the length of the separated region. 

J^e decrease and subsequent Increase in rms voltage after the first peak 
SL S! !®P®'®*^®^ 7) can also be explained by thiruisteadlnL^ 

region, which is where the measured skin friction was most 

the ''^®’^® *"^® remains separated most of 

llna^a^!u ^ affected by the increased voltage fluctuations due to the 

unsteadiness of the separation onset and reattachment regions. 

been ®f s^ck induced separation unsteadiness have also 

made during transonic airfoil testing. An example, provided by 

Ames, Pig. 9, shows a shadowgraph of a 64A010 airfoil 
In Jh« A ^ ®*^^® attack of 6 degrees and freestream Mach number of 0.6 

in the Ames 2 x 2 ft. wind tunnel. Included in Pig. 9 Is a power spectrum of 
ii® ® transverse laser schlieren system In the region of the shock 

e. The peak around 1400 Hz Is thought to be associated with similar 
larger scale motions of the separation zone. * 

Thus, from the turbulence modelling viewpoint. It Is important to 
turt^lent fleld^® ** l®t8e-scale contributions to the total 
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Measurement of Unsteady Plowflelds 

dlllOTal saitpltag 

experlments In fihtchihe tlowflnld sensor hnfl'l I”* «“‘'lcted to 

osn be used to generate the^Slwf « ^tlnuous output which Itself 

date. nessurleSn SI S nSSSS 1 “* ‘“'““““I eversgesi so that, to 

flows In which forlSSi- ^r i ^ ;*” to unidirectional she^r 

used. Shole clsseLlffJ^ l>ot wire snenonetry techniques c«, he 

have. thSefoSi bSnlegJS^;!"^^^^:^??^""? 1 «»»» 

to generate reliable analoa or diait-ai «. extremely dlfflctilt 

mentation, since they are ext^e^el^ outputs with conventional flow Instru- 
llnearlzed data interoretatlonr^^il sensitive to probe Interference and since 

flows. Thu^ reamrjSJJiJIt?v! TJ *^°w**^® ^ turbulent 

nonlntrusive linear tech^qJ^J? information can only be obtained using 

velocJSt«-t«SJ,°m SSSjJSlhlcTLTh'"*" s“^ **“*• 

tl^on the ph«.e .,.rsg«l and turhul«.t etructurSof t1S-SSndS.1Jr“' 

circuit c'Sl'SirS nSc™*2i1 1* *”%? * ' '“”'1 »” « 

numbsrs In the euheonlc and ttsSoSc1.glS?*’Thrcrtlos! **!™°lds 

with constant temperature surface hot fll* oIi... Inder was Instrumented 

The dynmslc gauge resSSe 1rSLr1hS1n*JS* “« ‘TP» >““1 1" «f. 2. 
tlon) was sufficient to determine the Mba with negligible phase dlator- 
These gauge outputs wSe uSrJrtSL!^! 1*“'! ••■U'idlng. 

which generated conditionally oamolaS*.»tfi*°'s*'^ scatter laser yelocimeter 
tlons in the unsteSy yoSS fC ShlS 

obtained by this new samnl-lnff "i®*^ the cylinder. Detailed information 

behind a clrcuUr cylinder Sa1f1S“? '^■^dpuudent mean flowfleld 

tore of Its waS wSe pSlSld turbulent struc- 

dlameSr^lSSSla^H tS1y!l.SS11'ri’,?‘“ '1“ ^.5 

conventional avetagJS ^outd ttJe f f® probable. Thus, 

large rms velocity fluctuation level which velocity close to zero and a 

changes in Induced mean flow velocity causefi hv° course, due to Instantaneous 
Any time-dependent information wouirbri^eJ ^ shedding, 
the probability density ^ L centerline (y/d - 0 . 5 ), 

velocities predominated as here the loJ^i ®^^hough negative vertical 

shedding frtm thrippe; surflll In ttf. determined more by vortex 

are more 
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v/d‘!!‘ flowfleld. However, outside the woke f . 

evereg^d’dlt; ""trlbutlone octur, ti 

dletrlLM^^r^et'Se P"® ‘’‘“"‘al velocity 

Since the flow repeats Itself oerlLlrrii^** “sing the thin film gauge output. 

the shedding vortices are at sJme^glJen p^ktLrin^th^^fl*"**® 

the time-dependent thin film as determined by 

sample again, and thus, over many cycles*\uild“iran 

constant phase. These velocli-i«a ,.L* * build up an ensemble average at 

random behavior of the fl^TiJ snd small-scale 

vortices fynZll aI ^ ® flowfleld with the 

Shedding cycle are sh^^ln^wP obtained throughout the 

sanq>led vertical velocity varlatlons'^sr^^ conditionally 

parlod corr..poudlS '-‘‘p 

«l/d“i ITT ”• 

the apparent vertical velr^n 4 ^v «. a..* 2. 5. As mentioned previously, 

of th^wtaf H~^“i, “*•’ «"«' 

fall below the axial small-scale turbulence data 

r "o* e 

X'r ^Lr^o'Sy' -r ^ 

couuuy ‘ turbulence. In general. It is 

t^tch rms velocities “-“ j T *?,, ‘°, ““ turbulence nodels scaled to 

unetesdy phenomena are llkelv to h. P°“ventioMl manner In any flow where 
prime eiemple. ^ ° *"~“tered, separated flows being a 


CONCLUDING REMARKS 


Hot wire, surface hot film, dynamic pressure and ia«„ i 
ment techniques have been described and i-Lff if ° 7 doppler measure- 
freestream turbulence tran!?M«! \ f applications in the areas of 

measurement w\"“’dS:::‘:j?"’ “"““'‘P ««» 

aLlhSr'“relatl1ns^r^r^^^^^ 

flight envelopes be aLessed. ^ encountered in particular 

complete^plcture°of^tMna<^7^^^a'^ transition it has been found that a more 
Piete picture of transition dependence on Mach and anlt Reynolds numbers 
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can be obtained by measuring the change in the root mean square of the voltage 
fluctuation across surface thin film gauges operated at constant temperature. 
This nonlntruslve technique not only enables the beginning and end of the 
transition region to be located, but also enables the turbulent intermlttency 
distribution through the transition region to be determined in both low- and 
high-speed flows. 

Compressible turbulent boundary-layer measurements have also shown the 
importance of spectra and length scale measurements of the mode fluctuations. 

In compressible flows large differences were detected between the mass— flow 
and total-temperature fluctuations. Not only were the scales different, but 
their probability densities and skewness were significantly different across 
most of the boundary layer. Long turbulence lifetimes suggest that future 
turbulence models must account for flow history especially when calculating 
non-equlllbrlxim flowflelds. 

In turbulent separated flows surface thin film gauges have been shown 
to provide basic Information on the significant unsteady character of these 
flows. Conventional "time averaged" techniques cannot supply this Information. 
The thin film gauges also Indicate a greater extent of the separated region 
since the onset and reattachment locations of separation are Intermittent. 
Time-averaged techniques can only locate regions where the flow is reversed 
at least 50 per cent of the time while the instantaneous thin film measure- 
ments are sensitive to regions which are separated for only a small fraction 
of time. 

A combined surface hot film and laser veloclmeter measurement technique 
which can be used to obtain new Information on the structure of time-dependent 
flowflelds has also been described. The data obtained in a cylinder wake show 
that mean and constant phase-averaged velocities can be determined. In 
addition, turbulence data can be obtained by conventional and conditional 
averaging of the velocity flucti .ttions. These data provide initial details 
of the small- and large-scale contributions to the total turbulent field. 


APPENDIX 


SYMBOLS 


D 

d 

e (f),e (o) 


f 

P 


w 



Re /in 
rms 


zz 



U 

c 

x.y 

bx^tiZ 

6 

P 


T 

< *> 


diameter of wind tunnel 
diameter uf cylinder 
mean square fluctuation levels 
frequency 

wall static pressure 

stagnation pressure 

correlation coefficient 

Reynolds number per meter 

root mean square 

stagnation tenqierature 

streamwlse velocity 

convection velocity 

streaimrlse and vertical coordinates 

streamwlse and lateral separation distances 

boundary-layer thickness 

density 

time delay 

rms value 
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Figure 5.- Autocorrelations in a turbulent boundary layer. 


4kHZ pQ»2iatm 
y/8*0.25 



6,“ Longitudinal cross-correlations in a turbulent boundary layer 




i i i iil i iM i yiw 









o SER^RATEO REGION 
— AFTER REATTAOIMENT 




IMSTAMTANtOUt VtLOCrrV. mftm 


Figure 10 .- Velocity probability deneities in a cylinder wake. 



Figure 11 .- Conditionally sampled velocities in a cylinder wake 
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Figure 12.- Root mean square velocity profiles across a cylinder wake. 
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HOLOGRAPHY AND LDV TECHNIQUES, THEIR STATUS 

AND USE IN AIRFOIL RESEARCH 

D. A. Johnson 
NASA Ames Research Center 
and 

W. D. Bachalo 

Spectron Development Laboratories, Inc. 


SUMMARY 


In recent experiments conducted In the NASA— Ames 2— by 2-Foot Transonic 
Wind Tunnel, the measurement capabilities of laser veloclmetry and holographic 
Interferometry In transonic airfoil testing were demonstrated. Presented In 
this paper are representative results obtained with these two nonlntruslve 
techniques on a 15.24-cm (6-ln.) chord 64A010 airfoil section. These results 
Include the density field about the airfoil, flow angles In the Invlscld flow 
and viscous flow properties Including the turbulent Reynolds stresses. The 
accuracies of the density fields obtained by interferometry were verified from 
comparisons with surface pressure and laser velocimeter measurements. 


INTRODUCTION 


Considerable advances have been made toward the numerical solution of the 
transonic flow past two-dimensional airfoil sections (e.g., refs. 1—5). How- 
ever, these numerical methods currently have a limited range of applicability 
either because they ignore viscous effects altogether or because they inade- 
quately predict the viscous effects as the shock wave strengthens on the air- 
foil’s upper surface. Accurate predictions are especially difficult to obtain 
once shock- induced separation occurs. The poor agreement with experiment in 
this case is believed due to deficiencies in the models employed for the 
turbulent Reynolds stresses (ref . 4) . 

Although, It Is the airfoil surface pressures, for which experimental 
data are plentiful, that the computational methods attempt to predict (the 
prediction of skin friction drag is much more elusive), quantitative measure- 
ments of the external flow are needed to truly assess the ability of the 
methods to describe the flow behavior and to provide insight into how to 
improve upon these methods. In view of this, there has been an effort within 
the Aerodynamics Research Branch at NASA Ames Research Center to establish 
measurement capabilities for the study of the flow past airfoils at transonic 
speeds. Two nonlntruslve techniques applicable to airfoil research that have 
been under development are laser veloclmetry and holographic interferometry. 

In this paper, the measurement capabilities of these two techniques in two- 
dimensional airfoil testing are described. The data presented are from experi- 
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of dftta obtainable with thaae two optical nathoda* a few representative 
results are presented. A more complete presentation of results with emphasis 
on fluid mechanical Interpretation can be found Inaa paper presented at the 
11th AIM Fluid and Plasma Dynamics Conference, Seattle, Washington. July 10- 


The symbols used herein are defined In an appendix. 


MEASUREMENT TECHNIQUES 


Laser Veloclmetry 

Figure 1 Is a schematic of the laser veloclmeter system for the NASA 
Mes 2- by 2-Foot Transonic Wind Tunnel. This fringe-mode veloclmeter Is a 
dual-color system utilizing the 4880 X and 5145 A lines of an argon-ion laser. 
One color Is used to measure the streamwlse velocity component; the other to 
measure the vertical velocity component. Bragg cell frequency shifting neces- 
sary for probing highly turbulent and separated flow regions Is Incorporated 
In both colors. It also facilitates the measurement of the vertical velocity 
component (l.e., +45® beam orientations to resolve the vertical velocity are 
unnecessary) . 


As seen In figure 1, most of the optical components are located outside 
the tunnel plenum chamber. There, color separation, Bragg cell frequency 
shifting, and the establishment of the four beam matrix are accomplished. 

transmitting lens, collecting lens, and photo detectors are mounted 
inside the plenum chamber. Two traversing systems are shown inside the plenum 
chamber. The one on the opposite side of the test section from the laser holds 
the collecting lens and photo detectors when forward-scatter light collection 
is used. This is the case in airfoil testing since there Is no problem with 
the model blocking the field of view. The traversing system on the laser side 
of the test section supports the transmitting lens and the light collection 
optics when backscatter light collection is used. Mirrors affixed to this 
traversing system permit three-dimensional scanning of the veloclmeter *s 
sensing volume; the optics outside the plenum chamber remain stationary. Both 
traversing systems are driven with computer controlled stepper motors. 


Signal processing Is accomplished with single-particle, burst counters 
and the Individual realizations from the two channels are simultaneously 
recorded with a digital con^uter. Hence, the velocity correlation u'v' can 
be obtained straightforwardly by multiplying and averaging rather than by the 
less accurate method of subtracting the variances of the two signals obtained 
from +45 beam orientations. 


Naturally occurring particles are sufficiently abundant in this facility 
to obtain data rates approaching several thousand per second. These particles 
arise from vaporized lubricating oil in the drive system which recondenses in 
the nozzle section of the tunnel. Flow response measurements across a normal 
shock wave have shown these particles to be 1 ym in diameter and smaller. 
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However, rates as high as 50,000 sec“^ can be realised by Injecting a DOP 
aerosol downstream of the test section wherein the entire tunnel Is seeded. 
This aerosol produced by an ultrasonic generator has a mean diameter of 
0,7 


Holographic Interferometry 

.Figure 2 Is a schematic of the holographic Interferometer setup used In 
the Amiss 2** by 2-Foot Transonic Wind Tunnel, As seen In figure 2, the holo- 
graphy tyystem Is designed to utilise the existing tunnel Schllereu mirrors. 

This ali( 3 W 8 a field of view equal to the dimensions of the tunnel test section 
windows (<i\l cm In diameter). The light source for this 8yst«*m Is a Q-swltched 
ruby laser v<d.th a pulse duration of a few nanoseconds which Is more than 
sufflclentl/y short to freese any mechanical motions. A series of mirrors 
shown on t^te opposite side of the test section from the laser are used to 
match path lengths between the reference and object beams. The dual-plate 
method Is employed wherein only one reference shot Is needed, thereby eliminat- 
ing the need to shut the tunnel down for each new test condition. Due to the 
relatively large span of the tunnel, the sensitivity of the system 
to density changes Is quite good In two-dimensional testing. For * 0.8, 
for exanq>le, one fringe shift corresponds to only nbout a 0.5Z change In 
density and a corresponding change In Mach uuinber of about IZ. This and the 
high aspect ratio of the wings tested in this facility (typically four), which 
Improves the two-dlmenslonallty of the flow, both contribute to the accuracies 
obtainable with this optical method. 


DISCUSSION 


A representative infinite-fringe incerferogram taken on a 15.24-cm (6-ln.) 
chord 6AA010 airfoil section at - 0.8 and an angle of attack, a » 0®, Is 
shown in figure 3, The fringes of the Interferogram represent lines of constant 
density. As noted earlier, the change in density between adjacent fringes 
is approximately 0.5Z. 


Holographic Interferometry, unfortunately, does have the disadvantage of 
being sensitive to flow disturbances along the total optical path. For most 
practical applications, this limits quantitative measurements to two- 
dimensional and axlsymmetric flows. Its range of applicability becomes even 
more questionable when one considers that most two-dimensional experiments 
exhibit some three-dimensional effects. Where applicable though, one hologram 
as Illustrated lu figure 3 can provide detailed density information which 
would be very laborious to realise with a point density measurement device If 
It existed. A means for determining whether minor three dimensionalities in 
the flow severely degrade the accuracy of the Interferograms Is available 
through comparison with laser veloclmeter and surface pressure results. Where 
the flow can be considered to be Isentroplc (a good assumption due to the 
small total pressure losses across the shock wave at transonic conditions) , 
the density rt^sults obtained by holographic interferometry can be compared to 
velocity and pressure results. 
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Figure 4 showa a comparison between direct surface pressure measurements 
and surface pressures Inferred from an Interferogram assuming Isentroplc flow 
and constant pressure across the airfoil boundary layer for a • 3.5 . At 
this angle of attack, the boundary layer Is separated at the foot of the shock 
wave, but the airfoil Is not near a stall condition. As seen from figure 4, 
the agreement with on-centerllne surface pressure measurements Is very good. 
What three dimensionalities are present In the flow appear to have a negllgi'ble 
effect on the accuracy of the Interferogram. There Is a difference, however. 

In the peak suction ahead of the shock wave recorded by the pressure trans- 
ducers and t>*« interferogram. This may be due to differences in the averaging 
times of the two techniques. Further Investigation vliereln a number of holo- 
grams are taken for the same test condition are needed to resolve this point. 

An even more thorough evaluation of the accuracy of the Interfeirograms 
can be made by making comparisons with localized laser velocimeter mensure- 
ments. Figure 5 shows such a comparison for a « 0*, Mg, ■ 0.8. In this 
figure, Mach contours were obtained by linearly interpolating between atations 
where point velocity measurements were obtained with the velocimeter. Away 
from the airfoil, the agreement is seen to be excellent. Close to the airfoil 
at the midchord position where the shock wave is located, the agreement is 
not as good* The explanation for these differences may be the same as with 
the pressure measurements since the velocimeter results are obtained over a 
much longer time period. Overall, the comparisons of figures 4 and 5 demon- 
strate that the interferometer can provide quantitative information of suf- 
ficient accuracy to give new insight into the character of the flow past 
airfoil sections at transonic conditions. 

Interferometry, although a very powerful technique, cannot provide 
information, for example, on flow direction or local turbulence pri'jpertles. 

To obtain these quantities, the laser velocimeter technique must be utilized. 

Flow direction measurements in the Invlscid flow regions become extremely 
important when tunnel wall effects on the flow field need to be considered. 

Flow angle measurements obtainable with the laser velocimeter are Illustrated 
in figure 6 for a scan at a fixed height (y/c * 0.167) above the chord line 
at a - 0”. At this condition, both tunnel wall and viscous effects are 
minimal and agreement with theory should be expected as is seen in figure 6. 
However, as the angle of attack is increased, this has not been the case, 
primarily due to tunnel wall effects (the Ames 2- by 2-Foot Transonic Wind 
Tunnel has slotted upper and lower walls). To totally account for any wall 
effects in comparisons with theory when slotted upper and lower walls are 
used, it appears that far-fleld flow angle measurements with the velocimeter 
will be needed to establish valid boundary conditions. 

The pacing item in advancing our ability to predict the transonic flow 
past airfoil sections is the development of Improved turbulence models. The 
realization of this does not appear forthcoming by numerical experimentation 
on large-scale computers devoid of any new physical Insights about the flows. 

If a solution to the turbulence modelling problem for transonic airfoils is 
to be realized, it seems that this will come about from measurements of the 
quantities that need to be modeled. Until recently, the fluid dynamlcist ^d 
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not have the means to measure these flow quantities (l.e#, the turbulent 
Reynolds stresses) on a transonic airfoil even when the flow was attached. 

Now» In principle, with the laser veloclmeter, this measurement capability 
has become a reality even In regions of separated flow. An Illustration of 
the laser veloclmeter's capabilities In measuring the turbulent flow proper- 
ties Is given In figure 7. The profile data shown were obtained In the 
separated flow region for an angle of attack close to Ct . Measurements 
were realized within 0.5 mm of the wing's surface. Note the smoothness In the 
Reynolds shear stress (pu'v') distribution (the mean de nsity has not been 
added) A check on the correlation coefficient, R^y - u'v'/u'v', at the 
point o f max imum shear shows It to be nearly -0.5. Also, the mixing length, 
i/6 ^ (u'v')I'2/[(3u/3y)5]^ at this point Is approximately 0.09. The self- 
consistency of these two results, which Indicate the flow to be In near- 
equlllbrlum In this part of the layer, supports the validity of the shear- 
stress measurements. From data like that shown In figure 7, obtained along 
the alrfo**.! and In the wake. Improved turbulence models can be formulated for 
airfoil flow field predictions. 


CONCLUDING REMARKS 


Laser veloclmetry and holographic Interferometry have been known to offer 
great promise In the study of the transonic flow past two-dimensional airfoils. 
In recent experiments in the NASA Ames 2- by 2-Foot Transonic Wind Tunnel, 
the measurement capabilities of these two techniques In transonic airfoil 
testing were demonstrated. A sample of the results obtained for a 15.24-cm 
(6- In.) chord 64A010 airfoil section has been presented in this paper. The 
detailed flow field Information that can be realized by these measurement 
techniques should provide the und'irstandlng needed to formulate improved turbu- 
lence models for airfoil flow field prediction methods. 


la 


APPENDIX 


SYMBOLS 


Measurementa and calculatlona were made In the U.S. Cuatomary Unite. They 
are preaented herein in the International Syatem of Unite (SI) with the equiva- 
lent values given parenthetically in the U.S, Cuatontary Unite. 
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airfoil chord 
preaaure coefficient, 


P - P« 
<1 


local static pressure 
free-stream static pressure 
free-stream Mach number 
velocity correlation coefficient 
streamwise coordinate 
normal coordinate 

free-stream streaoiwlse velocity component 

mean streamwise velocity component 

mean normal velocity component 

r.m.s. value of streamwise velocity component 


r.m. s. value of normal velocity component 

velocity correlation 

angle of attack 

boundary- layer thickness 

mixing length 

mean local density 
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Figure 2.- Schematic diagram of holographic interferometer installation in 
the Ames 2- by 2-foot transonic wind tunnel. 
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Comparison of surface pressures Inferred from an interferogram 
with surface pressures measured directly. 
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Figure 5.*> Comparison of interferometer and laser velocimeter results. 



Figure 6.- Flow>angle determinations with laser velocimeter. 
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SECTION DRAG COEFFICIENTS FROM PRESSURE PROBE 
TRAVERSES OF A WING WAKE AT LOW SPEEDS 

Lawrence C . Montoya , Paul F . Bikle* , and Richard D . Banner 
NASA Dryden Flight Research Center 


SUMMARY 


An in-flight wing wake section dre^ investigation was conducted using trav- 
ersing pitot and static probes. The primary objective of this investigation was to 
develop measurement techniques and impro/e the accuracy of in-flight wing proDle 
drag measurmnents for low values of dynamic pressure and Reynolds number. Data 
were obtained on a sailplane for speeds hrom about 40 knots to 125 knots at chord 

Reynolds numbers between 1 X 10 and 3 X 10 . Tests were conducted with zero flap 
deflection, deflected flaps, and various degprees of surface roughness, and for smooth 
and rough atmospheric conditions. 

Several techniques were used to increase data reliability and to minimize certain 
bias errors. A discussion of the effects of a total pressure probe in a pressure 
gpradient , and the effects of discrete turbulence levels , on the data presented herein 
and other experimental results is also included. 


INTRODUCTION 


In the fall of 1973 a joint NASA-SSA (Soaring Society of America) flight experi- 
ment was initiated to define the wake characteristics of a low speed airfoil section. 

The primary objective of the investigation was to develop measurement techniques 
and improve the accuracy of in-flight wing profile drag measurements for low values 
of dynamic pressime and Reynolds number . This improvement in accuracy is neces- 
sary if small airfoil performance differmices due to such things as airfoil surface , 
ambient air turbulence , or Reynolds number are to be measured . The accuracy effort 
was prompted by the desire for aerodynamic efficiency , which improves fuel consump- 
tion , and the desire of the designer to choose the optimum airfoil for a given mission . 

This paper reviews tlie techniques used to increase data reliability and to mini- 
mize certain bias errors during a series of wing profile drag measurements perfomed 
in flight on a sailplane airfoil. The pitot-traverse method developed by Jones (ref. 1) 
was used in this study. Jones' equation was rearranged so that incremental pressure 
measurements rather than absolute pressures in the wake region could be used to 
determine section drag. This was done so that (1) only one transducer could be used 
to measure both the total and static pressures in the wake , thereby eliminating errors 
from additional transducers; (2) differential pressure measurements could be used 
instead of the difference in absolute pressure measurements (which made it easier to 
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scale the transducer to the pressure measurement range); (3) in*flight transducer 
tare readings could be routinely obtained; and (4) the wake edges were well defined. 

Unresolved questions concerning errors In the use of total pressure probes 
In this and other studies are discussed. 


SYMBOLS 


Physical quantities In this report are given In the International System of Units 
(SI) and parenthetically In U.S. Customary Units. The measur^mts were taken 
in Customary Units . Factors relating the two systems are presented in reference 2 . 

c chord at test station, cm (in.) 

C(i section drag coefficient 

D external diameter of circular pitot probe or external he^ht of rectangular 

pitot probe 

Pq f)*ee-stream static pressure, kN/m^ Ob/ft^) 

p^ free-stream total pressure, kN/m^ (lb /ft*) 

0 

p. wake total pressure, kN/m* (Ib/ft*) 

w 

p^ wake static pressure, kN/m* (Ib/ft*) 

qQ free-stream dynamic pressure, p. - p_, kN/m* (Ib/ft*) 

q^ wake dynamic pressure, p. - p , kN/m* (Ib/ft*) 

w ^ 

X distance along chord , cm (In . ) 

y vertical wake width, cm (in.) 

P pitot probe displacement 

6 boundary layer thickness 

AP =P^-Pq» Wm* Ob/ft*) 

Ap. = p. - p. . kN/m* (Ib/ft*) 

*0 V 
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Ay 

A 


vertical wake width increment , cm (in . ) 
turbulence scale 


AIRPLANE AND TEST CONDITIONS 


The test vehicle was a T-6 sailplane with a modified Wortmann FX61-163 airfoil 
(fiCT. 1) . The airfoil modification consisted of straightening the aft lower surface cusp 
re^on and incorporating a 17**percent flap hinged on the lower surface (insert, 
fig. 1) . The coordinates of the airfoil are presented in table I . A more complete 
description of the test vehicle is given in reference 3 . 

The wing surface finish was smooth, with a maximum waviness for the test 
region near the midsemispan of about 0.008 centimeter (0.003 inch) in a 5.08-centi- 
meter (2-inch) section of surface. 

Data were obtained for calibrated airspeeds from about 40 knots to 125 knots, 

g 

which provided in-fiight test section chord Reynolds numbers between 1 X 10 and 

g 

3 X 10 , respectively . Tests were performed on a physically clean wing (for flap 
deflections measured at the test sections from 3^ to -10®) , for smooth and rough 
atmospheric conditions, and for various wing surface conditions (which included 
600 and 220 grit sandpaper finish, dusty wing, simulated insect impacts near the 
leading edge, and boundary-layer trips at the 5-percent chord) . 


INSTRUMENTATION AND TECHNIQUE 


The pitot-traverse method developed by Jones (ref. 1) was used to convert the 
wake data into section drag coefficients. Some rearrangement of Jones' equation 
was made to take advantage of the differential pressure measurement scheme, which 
is described below . Details of the way in which Jones' equation was used are given 
in the appendix . 

The wake measurements were made 24.4 centimeters (9.6 inches) behind the 
wing trailing edge, a distance corresponding to about 32 percent of the 75.9-cemi- 
meter (29.9-inch) local chord. 

The wake probe had both total and static pressure heads , as shown in figure 2 . 
Alsc shown are the sources for the reference total pressure (a Kiel tube) and reference 
static pressure (a trailing boom) . The Kiel tube and trailing boom were mounted 
on the wing near the probe to remove lag effects. The probe traversed about 20.32 
centimeters (8 inches) above and below the wing's trailing edge at a rate of about 
7.62 centimeters (3 inches) per second. 

An important part of the wake probe unit was the switching function provided 
by the valve (fig. 3) . Pneumatic lag was minimized by incorporating a design that 
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resulted in a small internal volume. The design also minimized the possibility of 
pneumado leaks. In the incremental total pressure (App mode, the different 

between the reference total pressure sensed by the Kiel tube and the wake total 
pressure was m^siared. Thus, the probe transducer sensed the total pressure 
defect (Ap^) in the wake. When the probe moved outside the wake, both sides of 

the trjmsducer's senstog element were exposed to free-stream total pressure, 

tSfwZ fo** the transducer . This feature minimized 

the bias error for the transducer and provided well-defined wake edges. 

incrementel static pressure (Ap) mode, the difference between the wake 
static pressure and the reference static pressure sensed by the trailing boom was 
measured. Therefore, the difference between wake static pressure and the free- 
st^am static pressure was obtained from a direct measurement and from the appli- 
cation of the position error correction . The position error calibration for the trail- 
ing boom source was obtained from the airplane static pressure source. The air- 
plane airspeed system is discussed in reference 4. 

switching valve feature just described, the same transducer pro- 
vided both Ap^ and Ap . As a result , the wake static pressure bias errors were also 

minimized through toe in-flight tare measurements, which were made before and after 
toe incremental static pressure measurements when the switching mechanism was 

in toe Ap^ mode. A 11.72 kN/m^ (10.25 psid) low range pressure transducer, 

which had an infinite resolution and natural frequency of 5000 hertz, was used. 

Six Ap^ traverses and several Ap traverses were made in succession , while 

indicated airspeed, and consequently q^, was held constant. Individual wake probe 

traverses were made alternately in the upward and downward directions and were 
averaged for determining the section drag coefficients for each test point (see 
appendix) . The incremental static pressures through toe wakes are not presented 

herein, but they were obtained with toe same accuracy as the incremental total 
pressiir es • 

Ifrobe position and the pressures were recorded on tape by utilizing a system 
mounted on a shelf behind the pilot's headrest. Aircraft dynamic pressure, q^, 

was d^ermined from a calibrated airspeed system and also measured with a trans- 

^ the recorder package. Pressure altitude, air temperature, indicated 
airspeed, general atmospheric conditions (smooth or rough air) , and pilot comments 

werehandrecordedby the pilot for each test point. 

coefficient values at the test section were not measured directly but were 
determined from known aircraft lift coefficients (±1 percent) adjusted for measured 
tail loads (from pressure distributions) and nonu;.iform span lift distributions. The 
nonuniform span lift distribution was determined from measured flap and aileron 
deflections at a number of spanwise stations. 
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RESULTS AND DISCUSSION 


The data presented herein illustrate some of the results of the experiment , The 
section drag data presented (flight, computed, and wind tunnel) are for chord 
Reynolds numbers corresponding to those obtained in flight, which varied from 

1 X 10® at the lowest speed to 3 X 10® at the highest speed. 

Examples of the wake profiles obtained by using this method are shown in 
fi^Pe 4. Data beyond the wake edges are not included, since the total pressure 
differences were zero using the measurement method described previously. These 
profiles were obtained at a calibrated velocity of 44 knots. Note the low pressures 

(0.10 kN/m (2.1 Ib/ft^)) that the instrumentation must measure at these meeds. 
The ve^ low random scatter illustrates the good repeatebility and low magnitude of 
the random errors . The fact that the shapes and displacements are the same for 
alternate upward and downward traverses Indicates that lag effects were essenttaUy 


The ability of the wake probe measurements to reflect the effects of small changes 
to dynamic pressure is clearly demonstrated to figure 5 by the magnitude of the 
mcrease to separation as speed decreases from 42.0 knots to 41.5 knots. All of the 
increased separation occurs on the upper surface; the wake from the lower surface 
remains essentially unchanged. In-flight tuft photos (fig. 6) confirm these results. 
Drag polars were defined for each of five incremental flap settings from 3® to -10® 
as measured at the test section . The flap results were cross plotted and adjusted 
to zero flap deflection and are summarized to figure 7. These data represent 
approximately 360 wake traverses and were obtained from seven flights over a period 
of 6 months. The line represents a fairing of the data. Most data fall within about 
3 percent of the fairing, which gives another indication of the repeatability of the 
results. A few data points fall outside this 3-percent band. The four isolated points 
between lift coefficients of 0 . 4 and 0 . 7 represent the largest scatter (approximately 
14 percent) . The reason for the large scatter of these foxir points is not known. 

Figure 8 shows section drag coefficients firom the test airfoil with boundary layer 
strips, which were 0.63 centimeter (0.25 inch) wide, had 0.089-centimeter 
(0.035-toch) grit, and were located at the 5-percent chord. The figure also shows the 
fairmg of the clean wing data from figure 7 • The artificially fixed transition increased 
section drag to approximately double the level measured for natural transition , indi- 
cattog that there was some laminar flow on the clean wing. This finding is consistent 
with wind-tunnel results for this class of airfoil . 

With the flaps deflected 6®, the mean camber line and the maximum lift coefficient 
were essentially the same as for the unmodified airfoU (airfoil without flaps) . Flight 
data for both 0® and 6® flap deflections are shown in figure 9. This comparison 
shows the effects of the ^nodifications on measured airfoil performance . 


The flight data for the 6® flap setting is compared with wind-tunnel and 

(Without nap,) in figure XO. One wt of 

wind tunnel data la reported In reference 5; the other wind-tunnel date are reported 
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in reference 6 . The average difference in drag between these sets of data is about 
12 percent, although both sets of data were derived from the same type of tests and 
airfoil model. The only changes were in instrumentation. The 1972 (ref. 6) data are 
believed to be the more accurate . This assumption appears to be conhrmed by the 
comparison in figure 10, which shows that the flight and computed data tend to 
support the 1972 wind~tunnel data. The computed data are based on S()uire~Young 
computational methods. As flgure 10 shows, the agn*eement between the flight data, 
computed data, and 1972 wind-tunnel data is generally good. 

Section drag coefficients obtained during flight in rough air which was typical 
of turbulent convection are presented in figure 11. Even though the data are some- 
what limited (37 wake traverses from seven flights) , they indicate that the rough air 
did not increase the drag until stall speeds were approached, and then only because 
the lift coefficient extended into the higher section drag coefficient regions. 

Figure 12 shows the effects of several s\arface conditions on the airfoil section 
lift-drag polar . The clean wing results (solid line) which were shown previously 
(fig. 7) are for a standard 600 grit sandpaper finish. The artifically tripped bound- 
ary layer results (dasl.ed line, shown previously in fig. 8) are also included. The 
circle symbols are for a dusty wing (the wing was not cleaned after the sailplane 
was tied down for 3 we^s) , and ttie triangular symbols represent a standard 
220 grit sandpaper finish. The square symbols represent the drag coefficients that 
resulted -when roughness particles simulating insect impacts were placed near the 
leading edge. As expected, the insect impact simulation appreciably increased the 
section drag (the increase was roughly 40 percent of the increase caused by the 
transition strips) , while the other surface conditions had no detectable effects. The 
effects of insect impacts, as is well known, depend on the number and size of the 
insects and their distribution near the leading edge of the wing. 


DATA UNCEBTAINTY 


As shovm in data presented previously in this paper , the low level of scatter 
indicates that random errors were quite small. Systematic errors such as bias 
errors and lag errors were also essentially eliminated. 

Error in the section lift coefficients resulting from the adjustments mentioned 
in the Instrumentation and Technique section are systematic in nature and may 
approach 4 percent at high lift coefficients and 10 percent at low lift coefficients. 
These errors are not felt to be of concern because of the insensitivity of the drag 
levels to lift coefficient through the lift coefficient range of greatest interest (section 
lift coefficients less than or equal to 1.0) . 

The systematio error which remains of concern is that associated with the 
total pressure obtained from a total pressure probe in the environment of a wake 
or boundary layer where (1) there is a pressure gradient across the face of the 
probe, (2) there is a discrete turbulence scale factor, (3) there are viscosity effects, 
or (4) the streamlines are deflected by the presence of the probe . Literature on the 
subject (refs. 7 to 20) acknowledges these problems; a summary of the findings is 
listed in table II, which was adapted from reference 7. The findings shown in 
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^ ^ equivalent dlaplacement of the probe center- 

line required to obtain correct totol pressure values. 

condttloM ®*'® Inconsistent, even where test 

P differences can probably be attributed to the varylntr 

for te2 l>revlou8ly . The reference 8 results are 

for test conditions essentially identical to those used in acqulrintr the fliffht data 

K® ®/f section drag coefficients (for the data from this experi- 

ment) by 0.0004, or about 6 percent, over the entire range of section lift coefficients. 

av<«*1h °J f^*!*'* “I? ground facility results indicates Uiat such errors also 
nS data from those experiments. However, these errors have generally 

not been accounted for or even acknowledged. It is probable that the probe effects 

inSfferMit faelHHM under identical conditions or from Identical airfoils 

m deferent facilities. Ihe designer may have difficulty in selecting an airfoil from 

such performance experiments , which are subject to significant errors. The current 
energy situation is such that smaU performance differewes c^ 
when extrapolated to long distances or large numbers of aircraft . It is therefore 
the^ur ®*P®**^enters to Investigate and attempt to reduce the impact of 

identified. As a first step, measurement techniques could 
hoii^^v io!f probe configuration and size relative to the wake or 

bounda^ layer being measured. An attempt should also be made to determine some 
form of turbulence criterion at the probe measuring station; Ste^iSively^^^ 

eliminate this problem with new remote measurement 
techniques or equipment such as lasers . 

Uatp? authors should acknowledge the four factors 

listed, indicate which are unaccounted for, and provide some estimate of the result- 

■?ate" • “ 

It is commonly assumed that published airfoil data are accurate at least to with- 
?erformanfp°H?ff‘ sometimes make airfoil selections on the basis of 

® percent. Published data are not necessarily 

to ,o 


CONCLUDING REMARKS 


versfn^cr^ITt?! '^®^® ®®®^®" investigation was conducted using tra- 

I«Hon li ♦ 1 ^®^® pressure probes. The primary objective of the in^sti- 

fn measurement techniques and to improve the accuracy of 

in flight wing profile drag measurements for low values of dynamic pressure and 
Reynolds number so that small differences in airfoil performance could be 
determined . 
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could be obtained*u^g ^averSngtotal^n^ measurements 

nlques described. Small dralr P^ssure probes with the teoh- 

condltiona and ambient air tufbulenoe wwe weU d”efi?^^^^ 

drag. Data obtained in rouah fturbulen?» -I * cnpacts affected 

affected at conditions near the wina^Sall t* ®"*y 

in lift coefficient extended into the hi^drorr^ns ” variation 

eiw. . which arc not adnuJtah?“nde?rtSd k*’”"”® ^heae' 

acknowledged In similar merimenta ’ have not been accotmted for or even 

This has made metmingftil comparisons between°vartou?** fllffht or in wind tunnels. 

«««d uaer. u, ha m 
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appendix METHOD USED TO OBTAIN 
SE TION DRAG COEFFICIENTS FROM MEASURED INCREMENTAL PRESSURES 

Jones' equation (ref. 1) is as follows: 




dy 


Since 


Jones' equation can be rewritten as 


Cd 




Step 1: Determine (a direct measurement from the aircraft system. Includes 
a position error correction) . 

Step 2. Determine Ap^ (a direct measurement (see Instrumentation and Tech- 
nique section) ) . 

niqursLuon^r.®***”^® ireasurement (see Instrumentation and Tech- 

Step 4: 


+ Ap = (Step 1) - (Step 2) 


Step 5: 


Step 6: 


q^ = (Step 4) - (Step 3) 



V(Step 5) /(Step 1) 


609 



step 7: 


Step 8: 
Step 9: 
Step 10: 
utep 11: 



" ^(Stop 4)/ (Stop 1) 


1 “ (Step 7) 
(Step 6) X (Step 8)Ay 
S(Step 9) 


Cd = (Step 9) X I 
c 
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Coordinates for ~ 


x/c 

Upper surface t 
cm (in.) 

Lower surface, 
cm (in.) 

0 

0 (0) 

0 (0) 

0.00102 

0.3556 (0.14) 

-0.3302 (-0.13) 

0.00422 

0.8128 (0.32) 

-0.6350 (-0.25) 

0.00960 

1.2954 (0.51) 

-0.8890 (-0.35) 

0.01702 

1.8288 (0.72) 

-1.1684 (-0.46) 

0.02650 

2.3622 (0.93) 

-1.4732 (-0.58) 

0.03802 

2.9464 (1.16) 

-1.7526 (-0.69) 

0.05158 

3.5052 (1.38) 

-2.0574 (-0.81) 

0.06694 

4.0894 (1.61) 

-2.3628 (-0.93) 

0.08422 

4.6482 (1.83) 

-2.6416 (-1.04) 

0.10330 

5.1816 (2.04) 

-2.9210 (-1.15) 

0.12403 

5.6642 (2.23) 

-3.2004 (-1.26) 

0.44643 

6.1468 (2.42) 

-3.4544 (-1.36) 

0.17037 

6.5532 (2.58) 

-3.7084 (-1.46) 

0.19558 

6.9342 (2.73) 

-3.9370 (-1.55) 

0.22221 

7.2390 (2.85) 

-4.1402 (-1.63) 

0.24998 

7.4930 (2.95) 

-4.3180 (-1.70) 

0.27891 

7.6708 (3.02) 

-4.4704 (-1.76) 

0.30861 

7.8232 (3.08) 

-4.5720 (-1.80) 

0.33923 

7.8740 (3.10) 

-4.6228 (-1.82) 

0.37056 

7.8740 (3.10) 

-4.6482 (1,83) 

0.40243 

7.7724 (3.06) 

-4.6482 (-1.83) 

0.43469 

7.5946 (2.99) 

-4.5720 (-1.80) 

0.46733 

7.3152 (2.88) 

-4.4450 (-1,75) 

0.49997 

0.9850 (2.75) 

-4,2672 (-1,68) 

0.53274 

0.6040 (2.60) 

-4.0132 (-1,58) 

0.56525 

6.1722 (2.43) 

-3.7592 ( 1 ,48) 

0.50750 

5.7404 (2.20) 

-3.5052 (-1,38) 

0.62938 

5.2832 (2.08) 

3.27GG ( 1,29) 

0.66074 

4.8000 (1.89) 

2.9718 ( 1.17) 

0.69133 

4. .3088 (1.72) 

2.G924 ( l.CG) 

0.72115 

3. 9110 (1.54) 

2.4384 { l),9G) 

0.74995 

3.5052 (1.38) 

2,1844 ( 0.8G) 

0.77773 

:i.07.)-l (1 21) 

1.9304 ( 0.7G) 

0.8043.5 

2.0924 (l.rO) 

1.7272 ( 0.G8) 

0.82970 

2.3(i22 (0.93) 

1 .4478 ( 0.57) 

0. 8.53.50 

1 .9.'..*>H (0,77) 

1 .2192 ( 0 4H) 

O.S7.5'.:o 

1 .(i2r»t; ((Mil) 

O.Oitori ( f).39) 

0. S96i 4 

1.3208 (0.52) 

o.7i»2o ( o.:in> 

0.91571 

1.0922 (043) 

0..5:i34 ( 0.2n 

0 . 9329!* 

II.H890 (II 3.‘>) 

0.3302 ( 0.13) 

0. 'I4k4m 

(1 (iSTiK (0,27) 

0.2032 ( O.OH) 

tl. !0ilMj 

0.'»(I8(. (0.2(1) 

0 . 1524 ( O.or.) 

0 . 1 

o..to4‘( ((1.12) 

0 lolo ( 0.04) 

1 . imiloo 

o.o'iOK (0.02) 

0 O'lOK ( 0 0/) 






















OlffOtlNIIAL rOTAl racSSURE MODE 


PROBE tOtAl PRESSURE^ *TATIC PRESSURE 

REFERENCE tOIAl PRB$URe;;> \ ^ "SrV”"® 



DiFFERfNTIAL STATIC PRESSURE MODE 






Flgup 4.“ Typical total>preasure wake profiles, 
flap deflection, 0°; velocity, 44.0 knots; q^, 
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Figure 5 


TOTAL PRESSURE, IWft* 

0 12 3 



- Total-pressure wake profiles showing the development of flow 
separation. Flap deflection, 0°. 
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Figure 9.-^ Comparison of section drag coefficients between 0^ and 6^ 

flap-deflection data* 


1.4 ; 



Figure 10,- Comparisons of section drag coefficients between flight, 
wind tunnel, and computed results for the FX61-163 airfoil. 
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SECTION DRAncoFmar.T 

Figure 11,- Comparison of smooth and rout;li air .section drag 
coefficients. 0° flap defleci lun. 




Figure 12.- Section drag coefficient loi ^ u i i ; :uii.ice 
conditions. 0*^ flap dellectUu. 
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FLIGHT TEST TECHNIQUES 

FOR LOW SPEED AIRFOIL EVALUATION*^ 

M. J. Hoffmaiui, a. M. Oregorek, and G. S. Welslogel 
The Aeronautical and Astronautloal Research Laboratox'y 
The Ohio State University 


SUMMARY 


Techniques for In-fll^t evaluation of new airfoils by modifying a single 
engine general aviation aircraft and measuring and recording airfoil surface 
pressures^ airfoil wake pressures^ and aircraft angle of attach and airspeed 
are presented. Included are descriptions of the a5rcraft modifications, Instru- 
mentation, data reduction techniques. Illustrations of typical results and 
comments on new equipment for flight test applications. 


INTRODUCTION 


The NASA LS(l)-0413 airfoil section characteristics have been evalxiated In 
a flight test program (ref. 1). A single engine aircraft was modified. Instru- 
mented and flown, and pressure data vrais acquired, reduced and summarized. The 
program was successfully accomplished by Implementing certain effective techni- 
ques. The existing wii\g of a "Beech Sundowner" testbed was "gloved" over the 
existing full-span to the contour of the LS(l)-0413 airfoil by the Beech Air- 
craft Corporation, an active participant In the flight program, therelQr saving 
the expense of construction of an entirely new wing. The aircraft was instini- 
mented with existing equipment supplemented by specially developed pressure 
measxirlng systems. A sophisticated and efficient data processing scheme was 
developed to handle the large quantities of data. 

The purpose of this paper is to summarize the details of the techniques 
used to effectively complete the flight test program and to comment on some 
new Instrumentatioh systems that could enhance future flight test efforts. 


The symbols used herein are defined in an appendix. 


xSoine of the techniques described hei?eln trere developed and implemented while 
under contract to NASA Langley Research Center. 
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AIRCRAFT MODIFICATION 


BaaJtgroimd 


To DQtlofy the objective of dotormlnlnR tho ooctlon chnractoplotlcn of the 
LS(l)-0413 airfoil in flight, it wan flrot nocossary to change the wing of a 
teotbod aircraft, ohown in figwe 1, to tho appropriate contour. Two options 
wore available; build on entirely new wing, or modify the oxiotlng wing. De- 
signing, building and proving a new \«dng is, unfortunately, time conoumlng and 
expensive. Modifying the existing wing by bonding on a new surface seemed 
attractive for this particular program for various reasons: 

1) budget constraints favored the relatively inexpensive nature of this 
modification approach, 

2) structural integrity of the existing wing structure could be 
utilized, 

3) cable and tubing routing could be easily facilitated, 

4) total time to modify the aircraft wing was a fraction of the time 
of the former alternative, 

5) by use of bonding, a smooth finished surface could be obtained 
without extensive structural proving. 

Based on these observations the ‘'gloving" approach to wing modification 
was adopted. This simply meant that the new contour would be obtained by 
bonding formers to the old wing surface and bondirjg a new skin to those formers. 


Design and Fabrication 


The LS(l)-0413 was a 13% thickness ratio airfoil while the existing 
632415 airfoil of the Sundowner was 1555. Figure 2 shows the modified wing with 
the larger chord. In order to accommodate the modification in the easiest 
manner, the leading edge was extended by 17.5 cm (7 in.) and the trailing edge 
by 7.5 cm (3 in. ) to obtain a "good" range of center of pressure relative to 
the aircraft center of gravity (CG). The new gloved wing had an incidence 
angle 1.4 degrees larger than the existing wing, and the original linear wing 
twist of 2 degrees washout from root to tip was preserved. Also, an internal 
channel was incorporated into the modification for running cables and tubing 
from the wing without significantly disturbing the airflow. 

The modification was begun by stripping the painted surfaces of the exist- 
ing wing and ailerons. Balsa former^ c\ 2.5 cm (1 in.) thickness were then 
bonded to the "old" skin on 20.3 cm (S .’n.) centers with an epoxy type adhesive. 
Spanwlse stringers were used for contour uniformity with wing results as shown 
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In figure 3 and aileron reeulte in figure 4. These formers were wrapped 
with 0.5 im (0.020 in.) sheet aluminum which was also bonded, thereby leaving a 
smooth, unlfomn rivetleso surface. Also, external mass balances were used to 
statically balance the now modified ailerons. 


INSTRUMENTATION SYSTQJ 


Overview 


Once modified, the aircraft was instrumented as shown schematically in 
figure 5. Aircraft angle of attack was monitored as was dynamic pressure (not 
pictured). A scanl valve/cut-off valve system, operated by a remote controller, 
was used to acquire surface pressure data. A wake survey probe, sensing total 
and static pressui s. ■?nd a rotary drive mechanism were used to obtain momentum 
deficit information, ultimately resulting in drag coefficients. All these 
systems were powered from an Instrument rack equipped with power supplies and 
signal conditioning equipment. 


Sensors 

The angle of attack sensor (vane) pictured in figure 6 consisted of a 
15.2 cm (6 in.) stem fitted with a 3.6 cm (1.5 in.) fin, driving a one turn 
potentiometer. A collar was used to limit the sensor travel to 5 degrees nose 
up and 40 degrees nose down relative to the mounting boom. The vane was 
located 0.75 chord ahead of the quarter chord point. To correlate the vane 
reading to local section, a deck angle inclinometer (bubble level, figure 7) 
was used. This added piece of instrumentation allowed for the determination of 
local geometric angle of attack. 

To sense surface pressures along the chord a strip-a-tube belt was used. 
The belt was formed from 5.1 mm (0.2 in.) plastic tube arranged in a group of 
twenty. Each tube was plugged appropriately to obtain forty active lines from 
the twenty tubes. Orifices were located once the belt was fastened to tjxe 
aircraft wing surface with double sided tape. Lead in lines were run frun the 
belt sensor, through the wing channel, into the cabin smd connected to the 
scanivalve/cut-off valve system. Figure 8 shows a typical belt installation. 

The wake survey probe was actually two sensors - one total pressure and 
one static pressure sensor - which was rotated through the wake of the wing at 
one of two spanwise stations during a given flight. The sensors were separated 
frc»n each other by 5.1 cm (2 in.) in the wing spanwise direction, thus allowing 
total and static pressures to be measured at the same chordwise station (0.14 
chord aft of trailing edge). By properly locating the wake probe, either a 
baseline (basic airfoil) or aileron station could be surveyed. Figure 9 
illustrates this wake survey probe. 
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Supporting Equipment 



excitation voltagen and olgnal conditioning an Inairumont ra<-)- 
(flgi^r. 10) was constructed. The rack contained all neceoofiry power -junnlle- 

denodulators. A^rthe scan^ 

valve system was mounted to the rack, and the 7 channel ra arialor 
tape recorder used for data recording was mounted atop the rock. TIk* <otai 

aW fully equipped was 351 N (79 lbs. ) with api^oxlmatelltn. 

slono of 48 cm W X v3 cm H X '1 cm D (19 In. W x 21 In. H x 1 ^/ In. D). 

was instrumentation systems a remote controller 

^roi o? ?hf r? controller allowed full manual or automatic con- 

trol of the syst^s from the flight test engineer position. The controller 
clock used to sequence the scanivalv e/cut-off system and the 
pitch-pause motion (described later) of the wake survey probe. Once the 

intervene until all data 

bHeset forLothei recorded. The systems would then 

con^JSU condition. Flgtare 11 is a photograph of this 


DATA PROCESSING 


flight would produce twelve test points resulting In larre 

effieleitfv^L^nT a processing scheme was davel^d to " 

f 1 5 ^ completion of a flight, the Ff.f recorder 

l?l^nH ground based digital computer system (figure 

vtdeH^hv ® sipals played back by the recorder. Timing pulses pro- 

vided by the instrumentation controller greatly assisted In the diglti-^ing 
.equenee. FORT^ coded progreice were used to maMpulete the Lfdlgktod 
data ^d also allowed the operator to select options as to how the data ^o-ii,i 
presented. For 1, stance, plot/end printouts lould Sr^d!" 
tvnl.fl for each test condition end/or suniroary plots could be made. It 

of airir^?t 1*HI« 1 houf 


TYPICAL RESULTS 


Surface Pressures 

acc 4 1^1 partially reduced surface pressure data were pressure '• >- 

efflclent-chordwls® location plots shown In figure 14. fvo flight te.st 

shown by symbols and the corresponding analytic computation;;’ 

oJ auLk Lt comparisons are at matched a.igle 

hv matched lift. The two-dimensional angle of attack was obtained 

of Beer.f'^At^"^ C4 induced ^le calculated by a three dimcnslorml analvUe cod.- 

te-t data^irrAAn ^^® scatter in the flight 

•teot data is seen to be low and compares well with analytic calculations At 
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the hJehcr nnglc (lower dynamic preanure for the cone at hand) the acattcr la 
;il li'htl.y more due to the preonure belt Itcelf and the decreonlni; accuracy of 
th<* pressure meuoureraenta at low dynamic pressure. Overall, however, these 
results are highly acceptable. 


Wake Surveys 

Partial reduction of the raw data from wake surveys lead to pressiiTe- 
posltion plots similar to those shown In figure 15. The first plot shows a 
continuously scanned survey while the second plot shows the stepped or "pitch- 
pause" method of surveying. The latter method was generally used to eliminate 
an.f potential response problems. The result of the method was a physical 
averaging of the data due to a finite number of points through the wake (shown 
by the relative smoothness of the plot). Both total and static pressures were 
measured and presented and the static pressure variation Is seen as signifi- 
cantly different from free stream static (the reference pressure). Based on 
these kinds of plots, limits of Integration were chosen and drag coefficients 
produced . 


Baseline Lift and Moment 

Carrying the surface pressure data reduction to completion by Integration 
of pressure distribution resulted In lift and moment coefficients as functions 
of angle of attack (figure 16). Lift data from three spanwise stations are 
shown to coincide very well by applying the three-dimensional analytic Induced 
angle correction and con^jares well with the faired wind tunnel data (solid 
line) of McGhee, et al. (ref. 3). The lift coefficient data becomes somewhat 
scattered at low dynamic pressures again due to the lower accuracy of the 
transducers in that regime. A small error in the dynamic pressure measurement 
unfortunately comes through strongly In the final reduction. The moment co- 
efficients are scattered and are generally more positive than wind tunnel 
measurements. The significant deviation could be due to slight trailing edge 
differences between the LS(l)-0413 modified aircraft wing and the wind tuxmel 
model used for comparison. 


Drag Polars 

Full reduction of the wake pressures lead to the baseline drag polar shown 
in figure 17. The symbols represent two test flights taken almost one year 
apart. The wind tunnel data is again that of McGhee. The cases shown are for 
smooth wing surface and smooth model. Due to the varying Reynolds number in 
the flight data two bracketing wind tunnel cases are presented and the proper 
trend of the flight data can be seen. A similar drag polar is shown in figure 
18. Here, however, the wind tunnel model boundary layer was tripped at 7.5/? 
chord as was the flight test airfoil. The baseline drag polar is also shown 
for reference as a solid line. Very good agreement can be seen as the flight 
test data trends from almost exact agreement at lower lift/higher Reynolds 
number to good agreement at higher llft/lower Reynolds number. 
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' f ,"i/i l.Ii .n Drit.H 


l‘ro.;.:'jTV' |mmi. . .. ■ . v.’ wore altjo taken at a!lorf-’>n otatlonr!, 

l’lr'ur*o r,' t Ml;; i.f ibutlon tit an aileron otatlc.n with tho 

no] M ]ln- ].< it.i- •• • ■ ■ ■ ' ■ I'l’/i for f.-larlty. The "pinched in" (reduced 

lin.) na'uf>' ■ ‘ •< M- i.o /eip flow Ic evident. No attempt wan 

made t.. 'If.! ';, ‘ • ’ I •ititi lyMo induced angle computatlona at 

t]ie;a- id -i ! i < ■! i.: in’ ' ' ' ii,’’ of flow phenomenon. Figure 20 

nhov.'.: the !V| 1 !,. i- jf.yr Incurred by having gap flow. The 

c .'mpar innti i.' ,vl "i ; . dirl’oil (nolid line) and reprenenta a 45!? 

inerean'/ In 'U-ir . 


The ndd?td.-'F>M ; ■■■.,. ’-.ippri by developing and using the afore- 

nienUrned it , i -a';,' • • ' • ■!. led to the design and construction of at 

leant two ru’W i.-o' ',' ' "Irure 21 shows a rotating type probe 

nc"*h;inj:'m ■.v'll 1. ■! • • i nojuntable in the chordwise direction. The 

rotallnr na; ■ ;/.vo npanwlse stations to be surveyed. The 
drive ii. :vir I. •. ■ .. ) wide and high and 30.5 cm (12 in.) long. 

Overall l•■ni’•• t- platform is 55.9 cm (22 in.) i«^ile the 

t'jt.-il w. pe... : . . : . . forjjn formers (streamlining) is 26.7 

tJ (0 lb;;. ). 

T'. ■•rii 'I! • ‘ ■ • i fi-ducUon capability while in flight, a 

Digital Oai.:i • r: ;’yn‘*.em ( DDARS ) has been developed in- 

h'juno (r'.‘f. ■ ■ n. t:;.| 111 components of this system, the main- 
frame (oJi i"'' D*, and the ground based disc drives (on 

rirha). ■■ : 1 1 microprocessor, 32K words of memory, all 

nece;a:ary > ! ape drives and a 32 channel analog-to- 

dhtltaJ '• nvi ‘ ‘ o.ivo throughput using all channels). 

Also In.'Jiid' ; iapif. relays for on/off type controls. 

Conlal’i'-'d I'! •' . ■ • ■ nh'jm) are 6 bridge balance units a:id 

6 dlfrer<-nt!';j a-,.; .uill Honing. The ccanponents used in a 

flight ;•! ; • I.* cignal conditioner rack and the small 

console l.< f’-il r 1 1 „. airborne package is 600 N (135 lbs.) 

and Uie rnx!—i • • t ; o watts. The total volume required to 

loieiti/ ,lii‘ ■ .. ' Mibl'’ imH.ers (8.8 cubic ft.). 


■'.V oUn DBriERVATIONS 


iv " ' 

airfoil in a !':i !■ 

iiavc bo« 'p it i si-'i ■ >•; 

Tli<.> dal. a alum' .v'!' 
most of I il' ' .1 ' C.' i‘: 
dynamic pros.sni'o . • ■ 
problim. ■ 


• •: ' 1 . 0 1 I (.1 evaluate the performance of an 
M'orafi. modification and Instrumentation 
'■•iia ;'roL*ussing scheme and typical results. 

•o':; lo indicate the system was acceptable over 
. ' ■ ' .’iilng questionable only at the very low 
l■•PlS'i'!ccr replacement could cure that 
!• '.V . ystoms which can be effectively applied 
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APPENDIX 


SYMBOLS 

Measurements and calculations v/ere made in the U.S. Customary Units. 
They are presented herein in the International System of Units (SI) with the 
equivalent values given parenthetically in the U.S. Customary Units. 

c chord of an airfoil 

C, section drag coefficient, 


m 


P 

P» 

«e 

X, x/c 
®2D 


P • pQ 
<3co 


section lift coefficient, Section lift 

section pitching moment coefficient with respect to 0.25 chord. 

Section moment 

— 

static pressure coefficient, 

measured local pressure 

free stream static pressure 

free stream dynamic pressure 
Reynolds number based on chord 

distance ilong chord, non-dimensional distance along chord 
aiigle of attack in two-dimensional flow 
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Figure 1.- Testbed aircraft in flight. 
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igure 2.- Air foil -modification cross section. 
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Figure 3.- Balsa-wood wing formers. 



Figure A.- Aileron modification. 
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Finure 6.- Angle-of-attack senaor. 
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Figure 11.- Remote-control panel. 



Figure 12.- FM analog recorder with patch panel. 
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Figure 13.- Ground-based digital computer system 



Figure 1A«- Comparison of pressure distributions 



position 


Figure 15.- Typical wake-survey profiles 











VVWV- 

1.6 0.8 ] 


Figure 19.- Aileron-station pressure distribution. 
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Figure 21.- New wake survey mechanism. 



Figure 22.- Digital— data acquisition and reduction system. 
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IN-* (GHT THREE-DIMENSIONAL BOUNDARY LAYER AND WAKE 
MEASUREMENTS FROM A SWEPT SUPERCRITICAL WING 

David P . Lux 

NASA Dryden Flight Research Cmter 
SUMMARY 


Three-dimensional boundary layer and wake velocity profiles were meastired in 
flight on the supercritical wing of the F-111 transonic aircraft technology (TACT) 
aircraft. These data, along with pressure distributions, were obtained to establish 
a data base with which data obtained by three-dimensional analytical techniques 
could be correlated. Only a brief summary of the total data base is given in this 
paper. 

The data presented represent one chord station at a wing leading-edge sweep 
angle of 26®. They cover an angle of attack range from 6® to 9® at free-stream Mach 
numbers from 0. 85 to 0.90. A brief discussion of the techniques used to obtain the 
boudary layer and wake profiles is included. 


INTRODUCTION 


A recurring problem during the initial design of transonic aircraft is the lack 
of proven analytical methods to reliably predict three-dimensional, viscous flow 
effects. Many analytical methods have been or are being developed to pzedict these 
effects. A few of these methods are described in references 1 to 3. However, there 
is a lack of adequate full-scale Reynolds number , three-dimensional data with which 
to compare and verify these analytical methods. 

Therefore, an investigation of the three-dimensional boundary layer and wake 
flow characteristics was conducted at one semispan station on the F-111 transonic 
aircraft technology (TACT) aircraft. The purpose was to establish a data base with 
which three-dimensional analytical techniques could be correlated . 

The three-dimensional upper-surface boundary layer measurements were obtained 
at 96 percent chord and the wake measurements at 108 percent chord of the TACT wing 
for a leading-edge sweep angle of 26® . The data presented are for a range of airplane 
angles of attack from 6® to 9® at Mach numbers from 0.85 to 0.90. The chord Reynolds 
number of these data is approximately 20 million. These data provide examples of 
boundary layer and wake flow characteristics at, above, and below the optimum 
airfoil design conditions . A brief discussion of the techniques used in obtaining the 
three-dimensional data is included . 


SYMBOLS 


The angle of attack used in this paper is the aircraft angle of attack referenced 
to the wing reference plane. However, because of wing twist and aeroelastic effects, 
the angle of attack for the section characteristics is 4° lower than that for the airplane . 

B aircraft span at wing sweep angle of 26° , cm 

Cl Uft coefficient , Lift/qS 


h 

M 


P 

Poo 

q 

s 

V /V 


pressure coefHcient, (p - p^)/q 

distance above or below wing trailing edge, cm 
Mach number 

local static pressure, kN/m^ 
free-stream static pressure, kN/m^ 
free-stream dynamic pressure , kN/m^ 
wing area, m^ 

ratio of chordwise velocity to free-stream velocity 


Vy/V^ ratio of spanwise velocity to free-stream velocity 

c ratio of distance from leading edge to local chord length 

Y distance from aircraft centerline , cm 

z distance above and perpendicular to wing upper surface , cm 

a aircraft angle of attack referenced to wing reference plane 

q semispan station, 2Y/B 


EXPERIMENT DESCRIPTION AND TECHNIQUE 


The TACT test-bed airplane is a modified F-lllA airplane (fig. 1) . The modi- 
fication consisted of the installation of a new wing incorporating an earlier NASA 
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supercritical airfoil. Although the aircraft's maximum speed is in excess c* Mach 2.0 
the supercritical wing was designed for transonic speeds and represents a compromise 
between cruise conditions (M = 0.85, Cj^ = 0.47) and maneuver conditions (M ~ 0.90, 

Cj^ = 0.67) Si a wing sweep angle of 26® (ref. 4) . 


For the boundary layer and wake data presented in this paper, the wing leading- 
edge sweep angle was fixed at 26® where chordwise pressures were previously 
obtained . The boundary layer , wake , and surface pressure measurements were 
made at semispan station 0.705, as shown in figure 2. The boundary layer measure- 
ments were ir.ade on the upper surface using the fixed, trailing-edge rake, which is 
shown in fig^ire 3. This rake was approximately 11 centimeters high with 12 multi- 
oriHce probes to sense local flow velocity and direction (fig. 4) . The calibration 
technique used to determine total pressure and flow direction angles for this probe 
is described in reference 5 . For the present experiment , the rake calibration was 
checked with tests conducted in Ames Research Center's 2- by 2-Foot Wind Tunnel. 
These results agreed with the results of reference 5 . Airfoil wake measurements were 
made using the rotating probe shown in figfure 5. The drive motor, position trans- 
ducer, and pressure manifold were mounted under the trailing edge of the wing, and 
the probe head rotated behind the wing at 108 percent chord. Details of the probe 
head are shown in figure 6 . The calibration used to determine total pressure and flow 
angles for this probe head is described in reference 5 . Velocity profiles of the wake 
relative to the trailing edge of the airfoil were obtained by applying an axis transfor- 
mation to the local flow angles at the probe head as it rotated behind the wing. 

The rake and rotating probes described above and the associated sensors were 
found to perform adequately , exhibiting good repeatability for all flight conditions 
where attached flow existed. However, the 10-second rotation time limited the amount 
of wake data that could be obtained with the rotating probe . For some applications , 
a fixed wake rake may be more desirable . 

All pressures from the boundary layer rake and wake probe were measured by 
differential pressure transducers located in temperature-controlled transducer boxes 
within the wing . The transducer boxes were interconnected by a common refereme 
line with the reference pressure being measured by a precision , variable- capaci- 
tance , absolute transducer . The pressures and air data parameters were recorded on 
an onboard aircraft data tape . 

The boundary layer data were gathered by flying the aircraft at a quasi-stabilized 
Mach-altitude condition . Then , to obtain data for a range of angles of attack , the 
aircraft was flown into a pushover-pullup maneuver followed by a slow-rate windup 
turn . This series of maneuvers generated an angle of attack range from approximately 
2® to 10® during a time interval of about 20 seconds. The wake data were obtained 
at the same nominal conditions with one exception . During the windup turn maneuver, 
the pilot would incrementally hold each angle of attack for 30 seconds. This allowed 
rotating probe data to be obtained for a minimum of three complete probe cycles for 
each angle of attack . 
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RESULTS AND DISCUSSION 


Representative three-dimensional boundary layer and wake data from the experi- 
ment are shown in figures 7 to 10. Each boundary lay^r and wake velocity profile 
is presented with the corresponding wing upper-surface pressure distribution. All 
velocity distributions are nondimensionalized with respect to the free-stream velocity 
obtained from the aircraft's onboard air data system. A complete description of this 
system is gfiven in reference 6 . 

The pressure distribution and boundary layer velocity profiles for a Mach number 

of 0.85 and an angle of attack of 8^ at a dynamic pressure of 14.36 kN/m* are shown 
in figure 7 (a) . Wake velocity profiles ara presented in figure 7(b) for the same 
flight condition . The chordwise pressure distribution data suggest that this flight 
condition is nearly optimum for upper-suiface flow . There is a long chordwise region 
of supercritical flow that terminates in an aft shock wave. The aft shock may cause 
a separation bubble at its base » but the pressure makes a complete recovery at the 
trailing edge and there is no second velocity peak . This flight condition is the basis 
for all other data comparisons in this paper. 

Note that in figure 7(b) the chordwise wake edge velocity' ratio exceeds 1.0. 

This is attributed to a static pressure error induced by the location of the static 
pressure orifices only two diameters behind the shoulder of the probe head (fig. 5) . 
While this error does not affect the total pressure measurements, it causes the abso- 
lute values of the calculated velocities to be incorrect . This static pressure position 
error is evident in all the wake data presented . 

It is of interest to compare the semispanwise flow characteristics for the wing 
upper surface at an x/c of 0.96 with the corresponding resnlt« at an x/c of 1.08 
(figs. 7 (a) and 7 (b)) . Note that at an x/c of 0.96 , the spanwise tiow component is 
inboard (V /V^ = -0.1) throughout most of the boundary layer. In the wake, at a 

distance of more than 8 centimeters above the surface of the wing, the flow turns 
outboard . This general trend can also be seen in the data for the other flight condi- 
tions considered in this paper . 

Figures 8(a) and 8(b) show a comparison of the data for the optimum upper- 
surface flow condition (figs. 7(a) and 7(b)) with those obtained at an angle of attack 
of 6® . At this angle of attack , the shock wave is farther forward and a strong second 
velocity peak is followed by recouipression to the trailing edge . Relative to the data 
for an angle of attack of 8® , the boundary layer velocity distribution for an angle of 
attack of 6® shows a thinning of the boundary layer and the wake velocity profiles show 
a similar thinning of the wake . This is , of course , indicative of lower section profile 
drag. While this lower section profile drag may seem more desirable, it should be 
remembered that less section lift is generated at an angle of attack of 6® , and airfoil 
efficiency is determined by the ratio of maximum lift to drag. 

It is of interest to relate these resulis to corresponding wing boundary layer and 
wake data from the F-8 supercritical wing airplane (ref. 7) . The F-8 data showed 
that the boundary layer and wake were thinnest at the airfoil design condition . This 
difference between the results of the present experiment and those obtained with 
the F-8 supercritical wing airplane is attributed to the fact that the F-8 supercritical 
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winpr was designed for a specific flight condition (M = 0.99, C, « 0.40), whereas 

in the chordwise velocity distribution In addition^'if nun ^®*®**^*^® edge condition 
veloctft, dUlribudon hJ "sVSmSVioS ' “ *’® 

indicate *'** *» ‘"‘“O"- 

not surprisin* whon one “nSdSS thifel, outboard) How. These results are 
as indicated by the pressure disMbution ***** *^* Oow has separated 

sLrt^distance forww^ fhe^ ase^oTth’ region extends only a 

«i?s“s»Si“ra“a5— 


CONCLUDING REMARKS 


number of 0.85 and m^ter“at*S fwm s" tfATthi 

nesses increase with increasing angle of attack. “"°ary layer and wake thick- 

In ^ner^, the data also indicate that there is little change in the throA-dimon 


647 


REFERENCES 


1. Bushnell, Dennis M.; Cary, Aubrey M., Jr.; and Harris, Jiiiius E.: Calculation 

Methods for Compressible Turbulent Boundai>y Layers—1976. NASA SP-422, 1977. 

2. Wortman, A.: Three-Dimensional Turbulent Boundary Layer Calculations - Exact 

and Simplified Solutions. AIAA Paper 75-854, June 1975. 

3 . Baker , A . J . : Finite Element Computational Theory for Three-Dimensional Boundary 

Layer Flow. AIAA Paper 72-108, Jan. 1972. 

4. Hallissy, James B.; and Ayers, Theodore G.: Transonic Wind-Tunnel Investi- 

gation of the Maneuver Potential of the NASA Supercritical Wing Concept. Phase I . 
NASA TM X-3534, 1977. 

5. Dudzinski , Thomas J . ; and Krause, Lloyd N.: Flow-Direction Measurement 

With Fixed-Position Probes. NASA TM X-1904, 1969. 

6. Cooper, James M. , Jr.; Hughes, Donald L.; and Rawlings, Kenneth III: Transonic 

Aircraft Technology - Flight-Derived Lift and Drag Characteristics. Volume I. 
AFFTC-TR-77-12, Air Force Flight Test Center, Edwards, Calif. , July 1977. 

7. Montoya, Lawrence C.; and Banner, Richard D.: F-8 Supercritical Wing Flight 

Pressure, Boundary-Layer, and Wake Measurements and Comparisons With Wind 
Tunnel Data. NASA TM X- 3544, 1977. 


648 



Figure 1.- F-111 TACT aircraft. Wing sweep angle, 26®. 
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Figure 2.- Location of boundary layer, wake, and pressure-distribution 
measurements on the F-111 TACT wing. 
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Figure 3.- Boundary- layer rake installed on wing of F-111 TACl aircraft 



Figure A.- Directional probe head for defining flow angle within the 

boundary layer. 
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(a) Pressure distribution and boundary-layer velocity profiles 



(b) Wake-velocity profiles. 

2 

Figure 7.- Section flow characteristics. M = 0.85: q “ 14.36 kN/m : 
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(a) Pressure distribution and boundary-layer velocity profiles. 




(b) Wake-velocity profiles. 

Figure 8.- Section flow characteristics. M - 0.85; q • 14.36 kN/m^; o 


0 20 40 60 80 100 

x/c, PERCENT 


''.4 .6 .8 1.0 -.2 0 


.2 


V /V„ 

X «o y «o 

(a) Pressure distribution and boundary- layer velocity profiles. 



(b) Wake-velocity profiles. 
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Figure 9.- Section flow characteristics. M ■ 0.85; q • 14.36 kN/m ; o 
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A PROCEDURE FOR ANALYZING TRANSONIC FLOW OVER 43 

HARMONICALLY OSCILLATING AIRFOILS* 

Warren H. Weatherill, F. Edward Ehlers, and 
James D. Sebastian 
The Boeing Company 

SUMMARY 

Finii difference procedures have been successfully used to solve the steady transonic flow about 
airfoils and appear to provide a practical means for calculating the corresponding unsteady flow. The 
purpose of the current paper is to describe a finite difference procedure derived from the equations for 
the potential flow by assuming small perturbations and harmonic motion. The velocity potential is 
divided into steady and unsteady parts, and the resulting unsteady equation is linearized on the basis 
of small amplitudes of oscillation. The steady velocity potential, which must be calculated first, is 
described by the classical nonlinear transonic differential equation. The initial research on this 
procedure was presented by Ehlers in reference 1, and is a direct extension of the steady state 
difference procedure of Murman, Cole, and Krupp (refs. 2, 3, and 4). The current authors have 
described further research pertaining to the procedure in references 6 and 6. 

INTRODUCTION 

The intent of the research described in this paper has been to develop a means for calculating 
air forces for use in the analysis of flutter and other aeroelastic phenomena. Thus, there is a need 
for developing a relatively efficient computational procedure as well as, of course, providing 
adequate accuracy in the representation of the physical phenomenon. 

The procedure of this paper is intended to be intermediate in terms of computer machine 
resource usage and is based on a finite difference method. The assumption of small perturbations 
from a uniform stream near the speed of sound retains the necessary complexity for describing 
flows with local supersonic regions. The application of the perturbation velocity potential restricts 
the solution to weak shocks, which, for thin wings of reasonably good design, is not too limiting 
an assumption. When the flow is steady, the resulting nonlinear differential equation reduces to 
the welLknown transonic small perturbation equation. The unsteady differential equation is 
simplified by considering the flow as consisting of the sum of two separate potentials representing 
the steady and unsteady effects. The assumption of small amplitudes of harmonic oscillation leads 
to a linear differential equation for the unsteady potential with variable coefficients depending on 
the steady flow. The resulting air forces are thus superposable and may be directly used in 
conventional flutter analysis formulations. 

The effect of thickness is included in the steady flow analysis. The unsteady analysis is carried 
out for a wing of vanishing thickness but submerged in a velocity potential distribution resulting 
from the steady analysis. As formulated, the shock is fixed by the steady flow, and it does not 
move with the wing motion. It is noted that shock motion could be included in a linear fashion; 
see, for example, reference 7. 
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FORMULATION AND SOLUTION 


The complete nonlinear differential equation for the velocity potential waa simplified by 
assuming the flow to be a small perturbation from a uniform stream near the speed of sound. The 
resulting equation for unsteady flow is 

[K ♦ (y • D^t • (y + IVxl^xx + ^yy • ^2^xt + “ 0 (1) 

where K « ( 1 • M is the freestream Mach number of velocity Uo in the x-direction, x and y are 

made dimensionless to the semichord b of the airfoil and the time t to the ratio b/Uo> y is the ratio of 
specific heats and e is defined in terms of thickness ratio 8 as « = (8/M)*'^. 

The potential is separated into steady and harmonic parts and written as: 

Ifi - + <^i(x.y)e^"‘ (2) 

For the steady potential, we obtain the usual nonlinear differential equation: 

[K - (y + + <f>0 yy “ 0 

While the unsteady potential, <p \ , is given by a linear equation of the form: 

|[K - (y - + V’lyy * (2i«/€)^pi^ + q^t = 0 (4) 


where: 


q = o»^/€ - i«(y - 1)^„ 


( 6 ) 


XX 


Linear boundary conditions are applied on the slit in the plane of the wing for both equations. 


A computer program for solving the steady>state transonic flow about lifting airfoils based on 
equation (3) was developed by Krupp and Murman (refs. 3 and 4). The output of this or a similar 
program can be used in computing the coefficients for the differential equation of the unsteady 
potential. The similarity of the unsteady differential equation to the steady-state equation suggests 
that the method of column relaxation used by Krupp for the nonlinear steady-state problem should be 
an effective way to solve equation (4) for . Note that equation (4) is of a mixed type, being elliptic or 
hyperbolic whenever equation (3) is elliptic or hyperbolic. Central differencing was used at all points 
for the y derivative and at all subsonic or elliptic points for the x derivatives. Backward (or upstream) 
differences were used for the x derivatives at all hyperbolic points. 

For the set of difference equations to be determinate, the boundary conditions on the outer 
edges of the mesh must be specified. In the original unsteady formulation, these boundary 
conditions were derived from asymptotic integral relations in a manner parallel to that used by 
Klunker (ref. 8) for steady flow. A later formulation in reference 6 applies an outgoing wave 
boundary condition to the outer edges of the mesh. This boundary condition is numerically 
simpler to apply and, on the basis of limited experience, appears to provide better correlation (sec 
the following section on numerical accuracy). 




The preferred numerical solution approach to the resulting large order set of difference 
equations is a relaxation procedure, which permit, the calculation to be made a. a ocgucncc of 
relatively small problems. The initial solutions were obtamed using a line 

Snve«Lr to rterminca b, monitoring ERROR, the maximum change n the *e octy potcntml 
between iteration steps. ERROR is defined as the maximum value over all i and j of 


•ft 


(n> 


U 


Vl 


(n*l) 


U 


( 6 ) 


where is the unsteady velocity potential for the n’th iteration. 

Itentitl for the preceding iteration, and r is the relaxation factor. The solutio .1 was cons dered 
ronverged when ERROR « 10'\ In some cases, particularly for finer meshes and for the pitch 
mode, convergence was considered complete when ERROR ^ 10 . 

RESULTS 

The above procedure has been used for problems with mixed flow and low reduced frequences 
with generally satisfactory results. In this section, results, as calculated with this procedure, ar 
compared (1), with linear theory appropriate to a flat plate of vanishing thickness and (2). with 
experimental results for a NACA 64A006 airfoil with an oscillating quarter-chord flap. 

The flat plate results are shown in figures 1 and 2. Figure 1 shows the 
coefficient across a flat plate oscillating in pitch. Figure 2 shows t^ 

flat plate with a harmonically oscillating quarter-chord control surface. In each fi^re, ^esulte 
obtained using the finite difference procedure are compared with ®®"®®P®"f 
using a kernel function program. For these calculations, the program is the subsonic routine 
generated for NASA-Langley by Rowe, Winther, and Redman (refs. 9 and f 

fhe results from the two theories should match exactly. The differences may be attributed to the 
finite difference representation together with the limited solution region. 

Corresponding results for a wing with thickness are shown in ^^8^®® f * f Sm^n and 
from the finite difference theory are compared directly with experimental data from Ty^^^ 
Schippers (ref. 11). Figure 3 shows the jump in pressure distribution for the airfoil at M - 0.80 
with\ reduced frequency of 0.25. At this Mach number, the pressure distribution does not have 
shock and the singularity in the pressure distribution is solely due to the presence ®^ ^® ®®®^ 
surface hingeline. Generally, the pattern of the calculated pressure distnbution matches that of 
the measured values very well. The calculated amplitude of the real part exceeds the measured 
value near the leading edge for the real part and over the front 76% of the chord 
imaginarv part. Parallel results at M = 0.875 and w = 0.06 are shown in figure 4. Here, 
is mixed with the shock appearing just aft of midchord. The characteristic P^®®^® "®® ^® ® '‘'® 
; shock and the pressure singularity due to the presence of the hingeline are 

both measured and calculated distributions. Figure 6 presents corresponding results at M 0.. . 
The shock has moved aft to a position just in front of the hingeline, but the two pressure rises 
Ifill appear as Separate peaks Also, the amplitudes are in better agreement than in the preceding 

two examples. 

These examples are typical of results obtained for low reduced frequencies (e.g.. refs. 1. 5. 6. 

12. 13. and 14). The reason for the discrepency between theory and experiment is not known, but 
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Jh‘orv »■• ‘■“‘h. or t. unknown problema with tho 

RELAXATION SOLUTION STABILITY 

.toWltv‘o.Sl"”““^ ‘" • “«>or8 (refo. 6 and 6), aignificant 

r '"Jf 0"00""*ored with tho relaxation procedure, ured to rel,e tl^iZT 


The characteristics of the solution instability are as follows: 


L T’”” *’"'’ ‘"u o-xl «>«» » not inw.l,od with 

the presence or absence of transonic shock flow* 

of X, and the dimensions of the finite difference remon An 
of X, above which the^rei^JltionTl^^^^^ °f 

1/2 


'(’RITKAL 


(7) 


3. 


‘’'[l**Kh‘‘] 

Where a is the streamwise dimension of the mesh region, b is the height 


Xxx + Xyy + = 0 


( 8 ) 


For7he”recto*^^*l'^*”* waves for which homogeneous boundary conditions occur on the boundary* 


lA]{^,} = {R} 


(9) 
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with certain conditicns on the structure of A, converge only when [A] in positive de(initt\ and 
this holds for the niodel unsteady problem of eq, (8l when h\ is less than |^l 


Several concepts have been studied in hopes of moderating or removing the relaxation solution 
stability problem. For example, integral equation formulations currently in use for the linearized 
subsonic unsteady problems employ only the outgoing wave soluti<m for the Green's function. 
Boundary conditions designed to do this in the finite difference formulation do not appear to be 
adequate. An investigation of a variety of different kinds of boundary conditions has not resulted 
in any significant improvement in solution stability. Other concepts that we hav(' tried include 
( 1 ) using a coordinate transformation so that the boundary conditions in the physical plane at 
infinity could be applied to the outer boundaries of a finite mesh region, (2) replacing the 
iterative relaxation solution with a full direct solution, (3) using an overlapping subregion 
concept, (4) artificial manipulation of the elements in [A] in order to provide a bett(*r conditioned 
matrix, (5) using a sequential mesh refinement system* and (6) applying a mathematical 
technique for making [A] positive definite for values of above premultiplication 

by the conjugate transpose of [A] (see ref. 7), 

The studies performed on these concepts are described in reference 6, Only the full direc t 
solution appears promising. Recently, we have been investigating the cfl'ect of adding a small 
amount of viscosity by including a v^i^xx ® one-dimensional formulation. Although this 

study is not complete, results to date show little reason for optimism. 


In the full direct solution, the problem is solved '*all rt once’' rather than ''sequentially." 
Mathematically, this involves inverting [A], a process that should be possible except at values of 
X| for which A is singular; that is for values of Xt that are eigenvalues. This procedure has been 
explored with very mixed results. 

In light of the original formulation, it was natural to try a semidirect solution of the form 

where contains an element for each interior mesh point and the right-hand side, which 

applies the boundary conditions, is a function of the yrj distribution of the preceding (iterative! 
solution. The vector of unknown t^j’s is found directly but iteration is required to update the 
boundary conditions. Although this procedure was very efficient for small meshc's f»r which it was 
used (i.e., permitted by an available in-core solution routine), it was subject to the same type ol 
solution instability as the relaxation solutions. However, it is possible to rewrite thi* equations so 
that all unknowns are on the left-hand side of the equation and the solution may then bi> obtaim>d 
without iteration, i.e.. a full direct solution procedure. 

Tin* full direct solution procedure has been tested with both one- and two-dimensional 
problems, and no difficulty was encountered in obtaining solutions for values of X| wtdi ahovt* tin* 
critical values. However, the accuracy of these solutions, as measured against analytic 
calculations, is good for values of \] smaller than deteriorates rapidly as A| is 

incrcasi'd. This point is discussed in the next section. However, the full direct .solution does 
provide a means of obtaining solutions where the relaxation procedure does not work. 
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NUMERICAL ACCURACY 


The dincuBHion of pccuracy for thifl paper will be very limited^ but there ore Bovoral pointR that 
should be mode. First, as noted above, for certain values of Mach number and reduced frequency, 
jjenerally satisfactory results have been obtained when correlated with linear theory and with 
available experimental data. Second, it appears that the use of outgoing wave boundary 
conditions, together with a coordinate transformation, provides better correlation with linear 
theory than Klunker-type boundary conditions. This is demonstrated in figure 6. Note that two 
different mesh patterns have been used. 

The problems of accuracy noted in connection with the direct solution have been studied with 
both one- and two-dimensional formulations. First, the one-dimensional analysis shows the 
difference between an exact analytic analysis and the finite difference solution may be a function 
of the form of the boundary conditions used. The nature of the eigenvalues (i.e., values of Xj) of 
[AJ are determined by the kind of boundary conditions used to formulate the problem. If the 
boundary conditions result in eigenvalues that are real values of Xj then the error distribution 
(the difference between the analytic and finite difference solution) is singular for the values of X| 
corresponding to these eigenvalues. If the eigenvalues of [A] are complex values X^ then this 
singular behavior in error curve is suppressed since Xi take on only real values. This is shown 
very clearly in figure 7 where Dirichlet boundary conditions were used in the first case and 
Dirichlet and Cauchy boundary conditions were combined in the second case. 

The error between analytic and finite difference solutions can be shown to be generally 
proportional to h*X|^. This is illustrated in figure 8 and it is noted that a conventional truncation 
analysis would indicate that the error was proportional to 

The distribution of the jump in the pressure coefficient across a flat plate oscillating in pitch is 
shown for a freestream Mach number of 0.4 and two reduced frequencies in figures 9 and 10. The 
point of these figures is to show the accuracy problem encountered when solutions are calculated 
at values of X| above using the full direct solution. Since pivoting was used in the 

solution routine, it is currently thought that this problem is related to the limited number of 
mesh points used in the calrulation, due to use of ^vcore matrix inversion routine. The effect of 
increasing the number of mesh points will be investigated with the completion of an out-of-core 
routine currently under development. 
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MNKARTHCOHV 

• i^v fiRCAiNT METHOD 



Figure 1.- Jump in pressure coefficient Cp across a flat plate 
oscillating in harmonic pitch* M » 0.35; o) « 0.18. 



Figure 2.- Jump in pressure coefficient across a flat plate with harmonically 
oscillating quarter-chord control surface. M * 0.9; tx) » 0.120. 


664 



10 .O.e ^.2 0 0.2 0.6 1.0 

X 


(a) Real pressure distribution. 
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(b) Imaginary pressure distribution. 

Figure 3.- Jump in pressure coefficient across an airfoil with a harmonically 
oscillating quarter-chord control surface. M = 0.804; a' - 0.253. 



(a) Real pressure distribution. 



(b) Imaginary pressure distribution. 


PiBure 4.- Jump in pressure coefficient across an airfoil with hambnically 
® Lcillating quLter-chord control surface. M - 0.875; o) - 0.06. 



(a) Real pressure distribution. 



(b) Imaginary pressure distribution. 

Figure 5.- Jump in pressure coefficient across an airfoil with harmonically 
oscillating quarter-chord control surface. M = 0.9; ui = 0.06. 


MACH NUMBER 0.B 
REDUCED FREQUENCY = 0.06 

mCH ABOUT LEADING EDGE 
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Pigure 6«— Jump in pressure coefficient across a flat plate oscillating in 

pitch. M = 0.9; u) = 0.06. 



Figure 7.- Comparison of error curves for boundary conditions yielding 

all real and complex eigenvalues. 
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pressure^coefflclent across a flat plate oscillating in 





Figure 10.- Jump in pressure coefficient across a flat plate oscillating in 

pitch. M = 0.4: u » 0.89. 
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A NEW TWO-DIMENSIONAL OSCILLATING WING APPARATUS FOR 

UNSTEADY AERODYNAMICS RESEARCH 

Sanford S. Davis and Gerald N. Malcolm 
NASA Ames Research Center 


SUMMARY 


A new apparatus for experimental research into unsteady transonic flows is 
described in this paper. The apparatus, as Installed in the NASA-Ames 11- by 
11-Foot Transonic Wind Tunnel, impart full two-degree-of-freedom motions at 
reduced frequencies to 0.3, oscillatory amplitudes to ±2®, mean angles to 12®, 
Mach numbers to 1.4 and Reynolds numbers to 12x10®. The test wing is fully 
instrumented for dynamic waveform measurements and the data can be acquired, 
processed, and displayed in real-time with a new computational data acquisi- 
tion system. Following a description of the apparatus, sample data from a 
recently completed test program will be presented. 


INTRODUCTION 


Rapid advances in computational techniques have made it possible to compute 
the unsteady aerodynamics of oscillating airfoils using a variety of governing 
equations. Until recently, analytical and numerical solutions were only avail- 
able for incompressible flows (analytical solutions based on the Theodorsen 
Function) and linearized compressible flows (numerical solutions based on the 
Possio Integral Equation) . Today, solutions are available for a relatively 
complete spectrum of equations including the Euler Equations, the full Potential 
Equation, the Small Disturbance Transonic Equation (low-frequency approximation) 
and the linearized Small Disturbance Transonic Equation. For two-dimensional 
flows, only the effects of viscosity arc lacking. It is Important to have 
experimental data to compare with these computations. 

The only detailed experimental data which are available are the results of 
the recent NLR (National Laboratory for Research, Netherlands) investigations 
(ref. 1). These tests proved to be Invaluable for their physical Insights into 
the complex flow fields surrounding oscillating airfoils. For subsequent 
investigations, however, it would be desirable to obtain data at higher Reynolds 
numbers, more favorable ratios of chord/wind-tunnel height and more general 
airfoil motions than were possible with the NLR Investigations. Due to tht wide 
operating range of the NASA-Ames 11- by 11-Foot Transonic Wind Tunnel, this 
facility was chosen for the Installation of a new two-dimensional pitch-plunge 
apparatus for the study of oscillating airfoils. The results of the first test 
series with the new apparatus will be used to validate transonic theories, to 
compare with recent NLR teats, and to develop a data base for unsteady transonic 
flows. 
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GENERAL DESCRIPTION OF APPARATUS HARDWARE 


Arbitrary pitch-plunge motions require a more complicated wing support 
mechanism than for a simple pitching wing. The final design uses two splitter 
plates to form a two-dimensional channel in the tunnel. Figure 1 shows a sketch 
of the wing-splitter-actuator system as Installed in the wind-tunnel test 
section. The normal 3.4 m x 3.4 m test section is augmented with two splitter 
plates which are 3.4 m high x 2.8 m long x 0.025 m thick. In order to preclude 
excessive deflections of the splitters, side struts provide lateral support at 
mid-span. The splitters extend into the tunnel’s plenum area at the top and 
bottom where they are bolted to I-beam anchors. Access panels for the instru- 
mentation cables and four slots for the push-pull drive rods are included in 
the splitter plate. 

Four 0.0412-m- (1.625 in.) dlam carbon-epoxy push-pull rods connect the four 
corners of the wing to independently controlled hydraulic actuators. The 
hydraulic actuators, located in the lower plenum area, are supported by flexures 
and bear directly onto a massive concrete foimdatlon through the four support 
columns. With this design, the tunnel pressure shell does not have to support 
the oscillatory reaction loads due to the actuator’s motion. 

The wing model itself is free to pitch and plunge in response to the 
actuator’s command signal. It is restrained in the fore-aft direction by a 
pair of carbon-epoxy drag rods, and in the lateral, roll, and yaw directions 
by sliding cover plates which move with the wing on the inner surface of the 
splitter plates. The 1.35-m-span x 0.5-m-chord graphite-epoxy wing is exten- 
sively Instrumented near the mid- span station. 

The capabilities of the test apparatus Include sinusoidal oscillations over 
a frequency range of 0 to 60 Hz. The maximum angle— of— attack oscillation varies 
from ±2® at low frequencies to ±0.8“ at 60 Hz around any point along the chord 
axis from — to +*, and a vertical displacement in heaving motion up to 
±5 cm (2 in. ) . The mean angle of attack is manually adjustable from -5“ to 
+15“. 


The various components that make up the system just discussed will now be 
described in more detail, since the basic performance requirements dictated 
state-of-the-art design in many cases. Many of the components that will be 
described can be seen in the photograph in figure 2 which shows the system 
installed in the tunnel, and in the photograph in figure 3 which shows the 
basic system set up on a test stand to be described later. In the following 
description it may be helpful to refer to these photographs to visualize tho 
various components and their interrelationship. 


P. 
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DETAILED DESCRIPTION OF APPARATUS COMPONENTS 


Motion Generators 

The servo-hydraulic actuator system Is driven by two 11 kW (150 hp) 
hydraulic pump units rated at 4. 1 x 10 “ 3 n3/s (65 gal/mln) at 21x10® N/m^ 

(3000 psl) . Each of the four actuators consists of two separate pistons on a 
single rod enclosed In a dual chamber cylinder. The upper piston Is used for 
generating dynamic forces, the lower piston for load biasing. The load-bias 
system Is necessary to support the mean aerodynamic lift load, thereby reducing 
the power required to drive the dynamic piston. As static bias requirements 
change, the servo-valve system responds accordingly to maintain the required 
force output. Velocity and position transducers are combined Into a single 

physical unit with colls and cores axially aligned for mounting In the center 
of the actuator. 


Wing and Push-Pull Rods 

The wing model Is connected to the push-pull rods through flexure bearings, 
and the rods are In turn screwed directly Into the actuator pistons. The wing, 
as well as the push-pull rods. Is a lightweight graphite-epoxy structure. It Is 
designed to withstand a 2300 m/s^ (230 g) acceleration end a 44,000 N 
(10,000 lb) aerodynamic load. The push-pull rods are each capable of withstand- 
ing a 22,000 N (5000 lb) tension load. The flexures, which are mounted between 
the push-pull rods and the wing, are also designed for a 22,000 N (5000 lb) load. 
The flexures are also strain-gauged to provide a direct measure of the lift load 
on the wing. A pair of graphite-epoxy rods, mounted to the wing with a flexure 
support and attached forward of the wing to the splitter plates, provides a 
means of holding the wing in place In the splitter plate to counteract drag 
loads. These are capable of 6700 N (1500 lb) each. 

The wing Is instrumented with both static pressure orifices and dynamic 
pressure transducers. As many as 41 channels of static and dynamic data can 
be accommodated. These are all located approximately at midspan. Static 
pressure tubes are routed from the end of the wing (see fig. 4, which also 
shows the drag restraints discussed earlier), down through a cavity in the 
plate to the tunnel plenum chamber, and out an access port to 
scanlvalve/transducer units located outside of the tunnel shell. Dynamic 
transducers are mounted in the wing by Inserting the transducer (2.36 mm in 
dlam) In the end of a long plastic sleeve, which is In turn Inserted Into a 
cylindrical channel molded into the interior of the wing. The channel termi- 
nates at the center of the wing to an orifice communicating to the wing sur- 
face. The transducer wires are then routed out the end of the wing (see 
fig* 5), through the splitter plates and out through the tunnel walls to 
the data acquisition equipment in the tunnel control room. A single reference 
pressure tube from each dynamic transducer Is also inserted into the plastic 
sleeve and routed through the splitter plate to the scanlvalve/transducer 
assembly outside the tunnel. The dynamic reference pressure can be selected to 
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be the static pressure of the adjacent static orifice on the wing or any other 
convenient pressure (such as the tunnel static pressure). Six accelerometers 
have also been mounted Inside the wing, one at each of the attachment points 
of the four push-pull rods near the comers of the wing, and two at the raid- 
span near the leading and trailing edges. The actual motion of the wing can 
be determined from the accelerometer output and compared to the output of the 
motion transducers located In the actuator piston rods. 

Two airfoil sections have been tested in this series. An NACA 64A010 
laminar-flow symmetrical airfoil (shown Installed in the splitter plates in 
^^8* b) was tested to obtain unsteady aerodynamic data to compare with the 
numerical computation of Magnus and Yoshlhara based on the Euler equations of 
motion, which is being performed under contract to the Air Force Flight Dynamics 
Laboratory. A NLR 7301 supercritical airfoil section was tested to obtain 
unsteady aerodynamic data to compare with experiments at NLR. 


Splitter Plates 

Vertical splitter plates with tralllng-edge flaps and horizontal side 
struts form the support structure within the test section for the wing and 
connected apparatus (see fig. 2) . They each have a sharp leading edge and a 
movable tralllng-edge flap which is manually adjustable between +2° from the 
plane of the splitter plate. All testing has been done with the flaps at 0°. 
Horizontal side struts (see fig. 4) attach to the outside of the splitter 
plates just below the vertical center, and protrude through the test-section 
wall to the exterior tunnel structure. These provide stabilization to the 
splitter plates and eliminate any excessive deflection in the lateral direction 
due to aerodynamic loads. The splitter plates are installed with a 0.1° 
diverging angle from tunnel centerline to account for boundary- layer growth. 
There are openings (figs. 4 and 5) in the splitter plate to attach the wing* to 
the top of the push-pull rods which are centered in the four channels cut in 
the splitter plates. In order to seal these openings when the wing is oscillat- 
ing, sliding covers (fig. 6) are attached to the wing end plates and slide with 
the wing on the inside surface of the splitter plate. These are also made from 
graphite-epoxy to reduce weight and are teflon lined in order to slide freely 
on the surface. 

The splitter plates and trailing-edge flaps contain a total of 130 
pressure orifices distributed over the inside and outside surfaces of both 
plates. The Inside orifices are utilized to select the proper channel Mach 
number, and the outer ones, in conjunction with the inner ones, are used to 
monitor the loading on the splitter plates. While testing, accelerometers are 
mounted on the trailing-edge flaps to monitor any large or potentially destruc- 
tive flutter motions on the flaps or main splitter plates, such as might be 
produced from the oscillating flow behind the wing or naturally induced from 
the channel air flow. Previous testing in the NASA-Ames 2- by 2-Foot Transonic 
Wind Tunnel (ref. 2) have demonstrated the viability of the spxltter— plate 
concept. 
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Preteat Setup 


Since every part of thlfl ayatem waa new and there waa no teat Information 
available to Judge the performance and reliability of the ayatem aa a whole, a 
apecial preteat facility waa built to permit a detailed checkout program. 

Many of the componenta, Including the wing, puah-pull roda, drag reatralnta 
and the hydraulic actuator motion generator ayatem are new designs and would 
not be a satisfactory risk in the wind tunnel without preteat experiments. 
Figure 3 is a photograph of the assembly in the test area. A support structure 
was constructed to which the various components were attached. The hydraulic 
actuators are mounted at the base with the push-pull rods screwed in to the 
top of the pistons. The wing is then mounted on the push-pull rods with flex- 
ures and angle-o£-attack blocks between the rod end and the wing end cap. The 
drag restralr t is fastened on top of the rear flexures and the other end tied 
to the support frame. Lift loads were simulated by an Inflatable bag between 
the lower surface of the wing and a support cradle fastened to the support 
stand. Drag loads were simulated by a pneumatic cylinder coupled to cables 
and straps looped over the wing. A nearly complete envelope of test conditions 
could be evaluated on the test stand. In the early stages of the test check- 
out, a dummy wing constructed of fiberglass (shown in fig. 3) was used before 
risking the graphite-epoxy test wing. This proved to be an extremely valuable 
and low-risk method of evaluating the performance of the entire system. The 
only real limitations were that tW fiberglass wing was not stiff enough to 
prevent large deflections at the mid-span, particularly i.i heaving, at the 
higher frequencies (above 30 Hz), and was not strong enough to accept ^he 
maximum lift loads. A limited amount of testing was done with the carbon-epoxy 
wing before installation in the wind tunnel. 


DATA ACQUISITION SYSTEM 


In the past, multichannel unsteady aerodynamic data were acquired with 
analog tape recorders. Raw data were recorded and stored for future analysis. 
On-line analysis was restricted to a few selected channels using special pur- 
pose "boxes” to extract limited usable data from the great mass of incoming 
data. These systems suffered from long-time lags between acquisition and 
analysis and the high probability of unknowingly recording spurious data. 

In the present test a new computational data acquisition and analysis 
system was developed for on-line display of steady and unsteady aerodynamic 
data. Figure 7 depicts the main elements of the new system. It has the 
capability of graphically displaying thu firat-harmonic, pressure distribution 
(both magnitude and phase) due to arbitrary pitch-plunge motions of the air- 
foil, along wxth the conventional static pressure dlscribution. At the 
user’s option, an overlay of selected theoretical or experimental pressure 
distributions from computer-resident codes or from a dedicated data bank can 
be accessed. 

The system is centered about a Data General Eclipse minicomputer, a 
high-speed (500 kHz) multichannel analog-to-digital converter, a large 
capacity (92 Mbyte) storage device, and a graphics terminal. The software 
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nyHtom conolRtfl of npproxlmntcly 50 Independent Fortrnn’'Coded progrnnifl which 
nre controlled by two executive programs? one for dynamic data, the other for 
fltatlc data. 


Dynamic Data Acquisition 

A sinusoidal signal generator drives the four-channel hydraulic actuator, 
which In turn drives the four push-pull rods attached to the four corners of the 
wing. The actuator's control system Is adjusted to Impart the desired pitching 
and/or heaving motion to the wing. The motion of the four push-pull rods Is 
continuously monitored and is acquired along with the unsteady pressure data. 

The signal generator Is also used to trigger a pulse to Initiate the unsteady 
data acquisition process. 

The dynamic signals from 41 miniature pressure transducers are amplified 
and filtered before entering the analog-to-dlgital converter. Since the signal 
is periodic. It is possible to obtain good waveform samples with minimum storage 
per data point by signal-averaging the data. Theoretically, a periodic signal 
Is completely defined by just one cycle of data (e.g., a 40 ms record is all 
that is necessary to characterize a 25 Hz periodic oscillation). However, the 
experimental signal is usually so contaminated by random pressure fluctuations 
due to wind-tunnel turbulence and model vibrations that one cycle of data is 
not very useful. 

The signal-averaging technique Is implemented as follows: we have a 

sample waveform and a pulse train which is triggered at the same phase position 
for each cycle of the airfoil's motion. These timing relations are sho\m in 
figure 8. At time t , the sample waveform is recorded for T seconds. At 
time tjj + nT, the waveform is recorded again for T seconds, where T denotes 
the period. The process is repeated M times. These M samples, each being 
initiated by the phase-locked pulse, are then ensemble averaged to obtain the 
averaged signal. In the current experiment T is chosen to be slightly greater 
than one period; n = 2, and M = 100 is sufficient for a good average. At the 
user's option, the signal averaged waveform and the Mth realization for any 
selected channel can be displayed on the graphics unit. 

For on-line analysis, the first harmonic of the response is most useful. 

A simple Fourier analysis algorithm is implemented to extract the magnitude 
and phase information at the fundamental frequency. This data is displayed 
in tabular form on the graphics unit within 30 sec of the termination of data 
acquisition. This data is sufficient to determine if the unsteady data 
acquisition process was successful. If more on-line analysis is required, the 
first harmonic data may be displayed graphically in pressure coefficient form. 

The magnitude and phase of the chordwise pressure distributions on the upper 
and lower surfaces of the airfoil are displayed along with certain theoretical 
curves. The software package currently includes two theoretical options: 

(1) linear, incompressible small disturbance theory (Theodorsen function) and 

(2) linear, compressible small disturbance theory (Possio Integral Equation). 

For time-aff Iclent on-line analysis it does not seem feasible to include 
unsteady transonic codes on the current generation of minicomputers. 
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Alflo available for comparlaon are the reaulta of other InveHtigatlona 
(theoretical and/or experimental) which have been ntorod in the data bonk. 

For comparing with NACA 6AA010 data, the theoretical inveatlgatlonn of Magnun- 
Yoflhlhara are available. For the. NLR 7301 wing, experimental data obtained at 
NhR''Amate.rdam are available. It is posaiblo to obtain a comparlaon between the 
current data and the Bolected theoretical/experimental overlay in approximately 
45 HOC after the. termination of data acquisition. 


Static Data Acquisition 

The static pressures are sensed with a conventional system using pneumatic 
tubing connected to a pressure scanning valve. The electrical output of the 
pressure cell to which the unknown pressures are multiplexed are read with a 
digital voltmeter whose BCD output feeds directly into the minicomputer. 


Setting Test Conditions 

The splitter-plate arrangement used for the oscillatory airfoil test 
requires special attention with regard to the free-stream Mach number (M„) . As 
discussed in a previous report (ref. 2), the Mach number in the channel between 
the plates is not the same as computed from a static tap in the plenum chamber. 
In order to obtain the approach Mach number, the splitter plates are equipped 
with approximately 150 static pressure orifices which are distributed among 
10 rows above and below the plane of the wing on the inner and outer walls of 
the splitter plates. These pressures are also sensed by the scanning system. 

The computed Mach numbers on the splitters are displayed on the graphics unit, 
and the approach Mach number is selected interactively by fairing the horizontal 
cursor of the graphics unit to the data. Using this procedure, the velocity can 
be selected to ±0.002 in Mach number. Once the Mach number has been chosen, 
the static pressure distribution is displayed along with selected overlays. 

A static pressure distribution with overlays can be displayed in approximately 
30 sec after the raw data have been acquired. 


REPRESENTATIVE DATA FROM THE OSCILLATING AIRFOIL TEST 


Steady and unsteady pressure distributions were measured on two airfoils: 
(1) an NACA 64A010 and (2) a NLR 7301 supercritical. The purpose of the 
NACA 64A010 tests is to compare the measurements with numerical solutions to 
the invlscld Euler equations obtained by Magnus and Yoshihara. These calcula- 
tions were made under AFFDL (U.S. Air Force Flight Dynamics Laboratory, Wright- 
Patterson AFB, Ohio) sponsorship and were made available to NASA-Ames under 
the terms of a Joint NASA-AFFDL cooperative program on the measurement and 
analysis of unsteady transonic llows,^ The measurements on the supercritical 


.1. Olsen of the AFFDL coordinated these calculations with our 
testing program. 


677 


airfoil will be compared to similar measurements on an 18-cm-chord model that was 
recently tested at NLR-Amsterdam (ref. 1), The NLR data was supplied to NASA- 
Ames under a cooperative program between NASA-Ames and the NLR. - Both the 
NACA 64A010 and NLR 7301 comparison data were stored in the data bank for on-line 
comparisons with the experimental data. 

For the purpose of indicating the capability of the new testing technique, 
some results will be shown which indicate the versatility and economy of the 
on-line data acquisition scheme. A comparison between the steady and unsteady 
(first harmonic) calculations of Magnus and Yoshlhara with the experimental data 
will be shown for the NACA 64A010 airfoil. The usefulness of the on-line analy- 
sis technique will be demonstrated by comparing the NLR experiments with the 
current experiment for the condition of shockless flow on the supercritical 
airfoil. 

The data to be presented In figures 9 through 12 are copies made directly 
from the graphics terminal. These unedited results contain some spurious data 
points from plugged tubes, broken transducers, etc., but the value of these 
displays for on-line analysis will be evident. 

Figure 9 depicts the measured and computed static pressure distribution 
for the NACA 64A010 airfoil. The strength of the shock wave Is predicted 
quite well by the Invlscld theory. However, the measured shock position Is 
slightly upstream of the computed position. Further analysis of the viscous 
effects will be made using measured data at other Reynolds numbers. The 
calculations do not include the effect of wind-tunnel walls, and the measure- 
ments Indicate that the ratio of airfoil chord to test-section height for the 
current experiment was sufficiently small to preclude large interference effects 
due to blockage. 

Figure 10 shows a comparison between the experimental and calculated first 
harmonic unsteady pressure distribution. The mean conditions are those shown 
in figure 9. The model was oscillated at 33 Hz about the 0.25 chord locr tlon 
with an amplitude of ±1". The first harmonic response can be expressed as the 
first term of a Fourier series; 

p(x/c,t) = aj cos(wt) + bj sin(wt) 

where p is the unsteady pressure, x/c is the fraction of chord, u is the 
radian frequency, and a^ and bj are Fourier coefficients. The data in figure 
10 show the magnitude and phase of the pressure normalized by the dynamic 
pressure q^: 


phase Cp » tan“k“bi/a^) -it < phase < it 
Tlj deman supplied the NLR experimental data. 
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The phase reference for the Cp data Is arbitrary. The phase reference for the 
Cp ca..culatlons la the a motions The phase angle for the experimental data 
(x) has been shifted (®) so that both experimental and calculated results are 
keyed to the same reference* 

Upstream of the shock wave (indicated by a bump in the magnitude and a 
large phase shift), the magnitudes agree very well. This is probably due to the 
fact that the mean flows agree with one another and the unsteady response is 
tied so closely to the mean flow. At the shock, the measured peak lies somewhat 
below the calculated one due to viscous effects. Downstream agreement is again 
quite good. The phase agrees quite well upstream of the shock. As mentioned 
above, the Euler equation calculations were made for free-field conditions. The 
good agreement shown here indicates that wall interference is probably not a 
problem. These results are typical of other comparison data between the experi- 
ment and calculations. Data have also been obtained at conditions beyond the 
capabilities of current computer codes (e.g., where strong shock wave boundary- 
layer effects exist). Comparisons such as those in figure 10 serve as valuable 
baseline data for confirming the correctness of both the calculations and the 
experiments. 

Figures 11 and 12 depict the mean pressure distribution on the NLR 7301 
airfoil for two different angles of attack a. The NLR data were obtained at an 
angle of attack of 0.85*. This was the experimentally determined shockless 
condition in their facility. In figure 11, data from the NASA-Ames facility 
at this condition is shown in comparison with the NLR data. The presence 
or absence of a shock wave on the upper surface is hard to perceive, but the 
agreement Is not very good. After some trial and error, best agreement was 
found at a = 0.37* as shown in figure 12. As shown by the data key, the 
on-line analysis enabled the experiment to be successfully compared with the 
NLR data in approximately 4 min. Due to the time constraints, the unsteady data 
from the supercritical airfoil has not been examined very closely. A complete 
®^^“li®6 analysis of all the data Is currently underway. 


CONCLUDING REMARKS 


A new dynamic oscillation apparatus, capable of testing two-dimensional 
wing sections in transonic flow with motions ranging from a pure rotational 
oscillation to a pure heaving motion, has been developed for operation in the 
NASA- Ames 11— by 11— Foot Transonic Wind Tunnel. It provides unsteady pressure 
measurements at Reynolds numbers and reduced frequencies previously unobtain- 
able. Two airfoil sections have recently been tested, an NACA 64A010 and an 
NLR 7301 supercritical. Preliminary comparisons between theory and experiment 
show good agreement. A more complete analysis of the steady and unsteady data 
is underway. 
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Figure 1 



- Two-dimensional oscillating airfoil test apparatus Installed in 
NASA Ames 11- by 11- foot transonic wind tunnel. 



Figure 2.- Photograph of two-dimensional oscillating airfoil apparatus 
Installed in the NASA Ames 11- by 11-foot transonic wind tunnel. 






Figure A.- Photograph of wing end section, drag restraint, side strut, and 

splitter plate. 
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Figure 5.- Photograph of wing end section and access to instnimentation. 



Figure 6.— Photograph of NACA 64A010 airfoil model mounted between splitte 




Figure 7.- Block diagram of the computational data-acquisition scheme for 
the two-dimensional oscillating airfoil apparatus. 





Figure 8.- Timing diagram for dynamic-data acquisition. Upper trace: 
dynamic-data signal, T is slightly greater than 1 period. Lower 
trace: trigger for analog-to-digital conversion; T = Period; n = 2. 
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Figure 11.- Static pressure distribution on NLR 7301 supercritical airfoil. 
M “ 0.74, a = 0.8°. Record of 3/18/78: 19 h, 55 niin. 
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N79-20005 

SOKE CALClOJiTIONS OF TRANSONIC POTENTIAL FLOW FOR 

THE NACA 64A006 AIRFOIL WITH AN OSCILLATING FLAP 

Robert M. Bennett and Samuel R. Bland 
NASA Langley Roaearch Center 


SUMMARY 

A tnethod for calculating the transonic flow over steady and oscillating 
airfoils has recently been developed by Isogai. It solves the full potential 
equatiot; with a 8eml-inq»llcit , time-marching, finite difference technique. 
Steady ilow solutions are obtained from time asymptotic solutions for a 
steady airfoil. Corresponding oscillatory solutions are obtained by ini- 
tlatlng an oscillation and marching in time for several cycles until a con- 
verged periodic solution is achieved. In this paper the method is described 
in general terms, and results are compared with experimental data for both 
steady flow and for oscillations at several values of reduced frequency. 

Good agreement for static pressures is shown for subcritlcal speeds, with 
increasing deviation as Mach number is Increased into the supercritical 
speed range. Fair agreement with experiment was obtained at high reduced 
frequencies with larger deviations at low reduced frequencies. 


INTRODUCTION 

The flutter critical portion of the aircraft flight envelope gen- 
erally occurs at transonic speeds. This critical condition results from 
both the high dynamic pressures of operation and the dip in flutter speed 
or %ucket” that occurs at transonic Mach numbers. The dip in flutter 
speed is Influenced by airfoil thickness and shape and cannot be satisfac- 
torily treated by state-of-the-art aerodynamic analyses. Thus an important 
current topic in aerodynamic research is the development of methods for the 
calculation of unsteady aerodynamics for use in transonic flutter analysis. 
Many of the current efforts have built on the recent success of steady flow 
numerical flnive nlfierence solution procedures, and to date, have pri- 
marily been applied to two-dimensional airfoils as a means of evaluating and 
pjfining the analyses and algorithms Involved. In this paper, a method 
developed at NASA Langley Research Center by Kojl Isogai* (refs. 1 and 2) is 
described in general terms, and transonic results recently generated for the 
case of an oscillating flap arc presented and discussed. 


“NRC-NASA Resident Research Associate, on leave from National Aerospace 
Laboratoiy, Japan, 1975-1977. 
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SYMBOLS 


n 


c 

Cl 


c 

m 


C 

P 

C* 

P 

Ac 

P 

Im 

k 

M 

P 

*lco 

Re 

t 

U 

V 

x.y 


nondlmenslonfll speed of sound 
nondlmcnalonal freestream speed of sound 

airfoil chord 

flap hinge moment coefficient, taken about hinge line, 
in direction of 6^, (hinge moment) /q^^c^ 

lift coefficient, llft/q^^c 

pitching moment coefficient, taken about c/4, positive 
moment /q^c^ 

pressure coefficient, (p-p„)/q^ 

pressure coefficient for sonic flow 

difference in pressure coefficient, C - C 

^lower ^upper 

imaginary or out'-of— phase part 
reduced frequency, u)c/2V^ 

Mach number 
pressure 

freestream dynamic pressure, 

real or In-phase part 
nondimens ional time 

nondimensional velocity component in x-dlrection 
nondimene ional velocity component in y-dlrectlon 
free-stream total velocity 

coordinate distances 


positive 


nose up. 
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a angle of attack 

Y ratio of specific heats 

6^ flap deflection, positive trailing edge dovm 

Pjjg free-stream density 

^ perturbation potential 

0) frequency of oscillation, rad/s 


GENERAL DESCRIPTION OF THE METHOD 

The full potential equation for two-dimensional time dependent flow Is 
(a^ - - 2UV(j)^ + (a^ - 


- 2U4) - 2Vd) - 4 « 0 

^xt ^yt ^tt 


where 

4 « perturbation potential 
U ■ cos a + 4 

X 

V « sin ct + 4 

y 

a^ - a^ - 0.5 (y - 1) (24^ + -1) 


This equation is nonlinear because a, U, and V are functions of 4> and 
numerical finite difference techniques are generally used to obtain solutions. 
Use of the full potential equation yields a method that is Intermediate in 
completeness and computational effort between methods that use the Euler 
equations (ref. 3) and those that use the small disturbance equation (refs. 

4 at.d 5). One advantage of using the potential equation as compared with 
the Euler equations is that for the potential equation only the single vari- 
able 4 has to be stored, whereas for the EuJer equations p, p, U, and V 
must be stored. 
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A seml-lmpliclt finite difference method has been developed by Isogal to 
obtain numerical solutions of the full potential equation In a time asymptotic 
manner for both steady and oscillatory cases. In reference 1, the basic 
method and several results are presented with enq>hasls on verifying the method. 
Additional computed exanq>les are given In reference 2. 

The seml-lmpllclt finite difference technique alleviates some of the limits 
on time step, required for stability, that restrict explicit techniques. An 
unpublished addition to the algorithm described In reference 1 has also been 
made by Isogal which extends the permissible time step. It essentially Is an 
additional Implicit pass through the flow field. The time step Is still 
restricted, however, which lnq>lles that the con^uter time Increases as 
reduced frequency k decreases. The method also uses a rotated difference 
scheme to maintain numerical stability. 

A stretched rectangular Cartesian grid system similar to that of ref- 
erence 6 Is used to map the Infinite physical space to a finite compu- 
tational region. The airfoil motion boundary condition Is applied at the 
mean airfoil position. This assumption may restrict the valid range of 
amplitude of motion, but considerably 8lsq>llfles the coiiq>uter program as 
the airfoil and confutation grid remain fixed In time. In addition, the 
finite difference method uses a quasl-conservatlve shock-capturing differ- 
ence scheme to treat the moving shock waves. This treatment of shock waves 
ensures that the correct shock junf relations are maintained. It also 
slnfllfles the confuter program as the moving shocks are treated automatically, 
but they are smeared over a few mesh spaces. 

The current version of the confuter program Is dimensioned for a grid 
of up to 61 X 61 points. It requires 37 000 (114 Kg) locations of central 
computer memory. On a CYBER 175, operating under the NOS operating system with 
the FTN compiler, the program takes about 1/2 seconds of CPU time for each time 
step. About 3000 time steps are required to converge for a steady flow case 
and 1000 to 5000 steps for a typical oscillatory case. 


RESULTS AND DISCUSSION 

Several cases have been calculated for an NACA 64A006 airfoil with and 
without an oscillating quarter-chord flap. Experimental data for these 
cases have been obtained by Tij deman and Schippers (refs. 7 and 8). For 
all calculations a 57 x 57 point grid system was used. The upper half of 
the grid is shown in figure 1 In physical space. One finite y-grid line 

and the lines at Infinity are not shown. The mapping clusters points near 

the leading and trailing edges and near the airfoil surface, but no special 
consideration is given to shocks or hinge lines. For the grid system used 

(fig, 1), there are 29 points on each airfoil surface with the first and 

last points on the airfoil at O.Olc and 0.99c, respectively. 


The steady pressure distribution for the NACA 64A006 airfoil with a ® 

6 0 (considering airfoil thickness only) has been calculated for Mach numbers 


692 



ranging from 0.500 to 0.900. The results for Mach numbers 0.800 and 0.875 are 
shown in figures 2 and 3. Generally, good agreement exists for Mach number 
of 0.800 with some deviation from the experimental data at the trailing edge. 
The calculated results for a Mach number of 0.875 do not agree as well with 
experiment, primarily as a result of a calculated shock that is stronger 
and further aft than the experimental shock. These two cases Illustrate 
portions of a trend found in comparing the calculated and experimental 
results. Good agreement was found at subsonic speeds, with some deviation 
over the aft end of the airfoil, but a gradual and increasing deviation 
between calculated and experimental results is apparent as Mach nunfl)er 
increases into the supercritical speed range. The mismatch of experimental 
and computed shock locations has a strong effect on the unsteady results. 

The reasons for the mismatch of shock locations are not fully known but are 
possibly a result of boundary layer and wind ttinnel wall effects. The 
calculations of reference 3 show chat inclusion of porous wind tunnel wall 
boimdary conditions can change the shock location by the same order of 
magnitude as the difference between the free*air calculations and experiment 
shown here (tunnel height/chord • 3.06 for these experimental data). 

Results for the NACA 64A006 airfoil with an oscillating quarter-chord 
flap are shown in figures 4 and 5 for two values of reduced frequency k and 
for Mach number 0.875. The in-phase and out-of-phase parts of the first 
harmonic of the oscillatory pressure distributions are shown (fig. 4 and 5). 

For the higher reduced frequency, k •» 0.234, the real part has two peaks, 
one at the shock location and one at the hinge line, and the imaginary part 
has a single peak at the shock location. Calculated trends correspond to 
experimental trends, but differences in shock location are apparent. Isogai 
(ref. 2) has shown improved agreement for this case if the calculations are 
made at a Mach number of 0.860 rather than 0.875 so that steady shock is near 
the experimental shock location. For the lower frequency case, there is 
significantly more deviation, particularly for the Imaginary part. The 
lower frequency case Involves a larger shock motion, and the unsteady results 
are apparently more sensitive to the steady shock location as well. 

The calculated real and imaginary parts of the first harmonic of 
pitching moment, lift, and hinge moment coefficients (per unit flap deflection) 
are shown for several reduced frequencies in figures 6 to 8 and are compared 
with the experimental values of referencs 8. All calculations shown are 
for a 1.08-degree flap amplitude which is the value for the experimental point 
at k = 0.234. The calculated results fluctuate with reduced frequency, 
whereas such fluctuations are not apparent in the experimental data. There 
is fair agreement of the theory with experiment at the two higher frequency 
experimental points. However, at the lower reduced frequencies there is a 
sizeable deviation from the experimental data. The relative agreement with 
experiment might be anticipated from figures 3 and 4, as the coefficients 
shown are integrals of the load distribution. The imaginary parts of the 
coefficients are zero at k = 0, and thus there must be a rapid variation near 
k “ 0. Calculations for static deflection of the flap show that the coef- 
ficients are very nonlinear functions of the flap deflection. There may be 
similar nonlinearities in the low frequency unsteady results shown in 
figures 6 to 8. 




Although a large number of calculated points are shown in figures 6 to 8, 

® single grid slze» time step» and amplitude of oscillation have been 
used. It is desirable to evaluate the Influence of these parameters also, 
but such an Investigation would require a significant expenditure of computer 
resources. For the cases of figures 6 to 8, it took about 5 minutes of CPU time 
on a CYBER 175 to calculate 6 cycles of oscillations at k •» 0.500. This 
time Increased to about 30 t Inutes to execute 3 cycles of oscillation at 
k = 0.059. Although the method gives results that offer considerable Insight 
to transonic flows, It needs further development to reduce the computer time 
in order for It to be used as a production tool to study a wide variety of 
cases. 


CONCLUDING REMARKS 

The n^thod of Isogal for calculating the transonic flow over steady and 
oscillating airfoils has been used to calculate several cases for the NACA 
64A006 airfoil with an oscillating quarter-chord flap. The thickness pressure 
distributions were in good agreement with measured results at subcrltlcal 
speeds but gradually deviated further from the data as Mach number was 
increased. Calculated unsteady results for a Mach number of 0.875 showed 
large frequency effects not apparent in the experimental results. Fair 
agreement was obtained at the higher reduced frequencies with larger 
deviations at low reduced frequencies. The method gives results that offer 
considerable insight to transonic unsteady flows, but needs further develop- 
ment to reduce the computer resources required* 
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Figure 1.- Grid system used for calculations, 57 x 57 . 



2.- Steady pressure distribution for the NACA 64A006 airfoil 
at M = 0.800 and a » 0°. 




Figure 3.- Steady pressure distribution for the NACA 64A006 airfoil 

at M ■ 0.875 and a « 0°, 
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x/c 

(a) Real part. 

Figure A.- Load distribution for an NACA 64A006 airfoil with an oscillating 
quarter-chord flap; k = 0.234, M = 0.875, a = O.O”, and 16^1 = 1.08°. 




x/c 

(b) Imaginary part. 


Figure 4.- Concluded. 



x/c 


(a) Real part- 


Figure 5.- lx.ad distribution for an DACA 64A006 airfoil with an 

ouarter-thord flap; k = 0.059, M - 0.875. n - 0.0”, and |6j| - 1.08 . 


1.08‘ 








Figure 8*- Unsteady flap hinge moment due to flap motion 
at M = 0.875 and a = 0*^ with = 1.08*’. 
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OBSERVATIONS ON THE DYNAMIC STALL CHARACTERISTICS 
OP ADVANCED HELICOPTER ROTOR AIRFOILS 
L. Oadone 

Boeing Vertol Company 


SUMMARY 


A significant amount of research has been devoted to 
undersli^pdin9 ni6chdnlsiii of dynamic s^&ll dol^y as appli— 
cable to che flow environment of a helicopter rotor in forward 
flight. One aspect of such research deals with the unsteady 
characteristics of two-dimensional airfoil sections over a Mach 
number range from 0.3 to 0.6» since such characteristics 
can be meaningfully related to rotor performance and loads. 

This paper summarizes the results of several oscillatory 
tfsts carried out on conventional, transonic and BLC-equipped 

airfoils. ^ 


INTRODUCTION 


There are two reasons to conduct oscillating airfoil tests. 
First/ to identify the key elements in the unsteady flow environ- 
ment and/ second/ to determine the unsteady characteristics of 
specific airfoils. 


Until now an experimental approach has been necessary because 
d^amxc stall cannot yet be quantified by theoretical means. Al- 
though empirical dynamic stall representation techniques have been 
developed/ a fundamental understanding of the problem remains to 
be achieved. 

, considerable efforts have been directed to the 

depnition of transonic airfoils applicable to helicopter rotors. 
This research was encouraged by the availability of new methods of 
analysis which greatly reduce the difficulty in designing airfo’.ls 
with favorable transonic characteristics. 

Now we have a large body of information from the oscillatory 
tests of several airfoils, and we can at least identify the most 
significant trends in the data. A review of these trends pro- 
vides preliminary answers to questions concerning the comparison 
of transonic and conventional sections, further refinements in 
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airfoil design objectives, and generally, the purpose and useful- 
ness of future unsteady airfoil tasting. 


SYMBOLS 

Measuroments and calculations wore made in U.S. Customary 
Units and are presented in both the International System of Units 
(SI) and U.S. Customary Units. 

b airfoil semichord, m 

c airfoil chord, m 

Cj) normal force coefficient 
Cm pitching moment coefficient 

C^ blowing momentum coefficient, for sections employing 
active boundary layer control (BLC) 

f drive frequency of airfoil motion in pitch, Hz 

)c reduced frequency, irfc/v 

M Mach number 

Pij tunnel test-section total pressure, N/m^ (psia) 

t airfoil thickness, m 

V t\mnel velocity, m/sec 

a angle of attack, deg 

a first differential of a with respect to time, deg/sec 

Y stall delay function 

Ao amplitude of pitching motion, deg 

SOURCE OP DATA 


A thorough survey of current unsteady aerodynamics research 
is present in reference 1. However, until now the main source of 
oscillatory airfoil data applicable to helicopter rotors has been 
the two-dimensional subsonic insert of the Boeing Supersonic Wind 
Tunnel (BSWT) . 

702 



A sketch of the oscillatory rig in the two-dimensional test 

blow-down facility with provisions for variable density operation, 
so that tests can bo run at Reynolds numbers corresponding to 
full'-flealQ blade chords with 6'^ to 7-inch chord models# 

Solid floor and ceiling have been employed in all the os- 
cillatory tests. Quasi-steady baseline 
in the presence of both solid and porous (4.9%) 

The oscillatory tests have been conducted within the following 

ranges of conditions. 


Mach number, M 
Mean angle of attack, a 
Ampiltude of oscillation, A a 
Drive frequency, f 
Total Pressure, P-j 


0.2 to 0 . 7 
-20® to +20® 

2.5® to 10® 

12 Hz to 100 Hz 

2 

170,000 N/m^ to 520,000 N/m 
(25 psia to 75 psia) 


The models of the airfoils were equipped with differential 
pressure transducers mounted as close to the pressure ports as 
possible. Each transducer/orifice cavity was sized to place the 
acoustic resonance frequency well beyond tne tenth harmonic or 
the highest drive frequency of interest. 


T 2 d 3 le I lists the oscillating airfoil tests conducted in 
BSWT. Most of the tests have been documented and the data are 
available in references (2) through (5) . The boundary layer con- 
trol experiments on the V23010-1.58, and the VR-1 test have not 
been f'^rmally documented. Limited data from these tests are sho^ 
only to illustrate the relationship between static stall and ae 
dynamic damping. 

Figure 2 shows typical pressures and integrated 
data were acquired during the tests of 

ential pressure coefficients are displayed as time-histories, 
while the normal force and pitching moment coefficients are shown 
as a function of angle of attack. 
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KEY PARAMETERS 


Definition of Maximum Attainable Normal Force, Cnj^^x 


The highest normal force achieved while an airfoil undergoes 

sinusoidal pitching motions over conditions lelij So" 

tached flow to stall has been used as a measure of stall delay po 

tential. The maximum attainable ^it^f^SSt^beSS 

reduced freguency and Mach number. Until now, it has not oeen 

successfully predicted by either theoretical or empirical me^s. 
Because of telt limitations toe maximum no^l force cannot be al 
ways reached at model drive frequencies beyond 6Q Hz. 


Definition of Aerodynamic Damping 

The area enclosed by the Cni/ o trace, see figure 2, and the 
sense of mot'^on around the loop have an important physica ^ 
ni?“aS~ iSe net work done by toe airfoil on to. surrounding 
air is proportional to the integral 


W 


^ Cm da 


This integral is proportional to the area enclosed by the 
c a trace 2 mid it is positive for a counterclockwise circuit. If 
S!4 circSn Gileses a^ubstantial area in a 

contribution of that area is negative, i.e., it represents energy 
extracted from the airstream by the airfoil. 

For pitch oscillation about the quarter chord the theoretical 
damping is 

D^unping = ‘“k 

The theoretical damping is used to 
cycle damping value computed for each test condition. 
the derivation of the theoretical damping are shown in referenc 


Definition of Stall Delay Parameters 

Rftfarance (6; deb. *ibes the formulation of an empirical 
netho? to di.«mi= stall. This method was specif ioally 

dsflnsd for rotor blStUs calculations. 

for=eTnd1ItrhiJ^SSr,tL^^^^^^ 

srtSrsraiftSgreriis^So^Lsf^^^^ 
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has been. called "y- function.” This quantity is different for the 
normal force and the pitching moment, and it decreases with Mach 
number. Figure 3 illustrates the derivation of y-f unctions from 
oscillating airfoil data. 


REVIEW OP CHARACTERISTICS 


Maximum Normal Force 

.^’igures 4 and 5 summarize the largest normal force attainable 
by several airfoils at M = 0.4 and 0.6, respectively, over the 

range of reduced frequencies of interest for helicopter rotor op- 
eration. ^ 


The NACA 0006 and V13006-0.7 are thin airfoils (t/c » 0.06). 
The VOOll and V23010-1.58 are the airfoils employed on the CH-47A 
helicopters (t/c » 0.11 and 0.102) and, finally, the 

(t/c » 0.11) and NLR 7223-62 (t/c » 0.086) are transonic sec- 
txons . 


The M » 0.4 condition is representative of the flow environ- 
ment within which retreating blade stall takes place. The reduced 
frequency value of k * 0.06 is typical of 1/rev pitch changes. 

The maximum lift capability of all sections increases with re- 
duced frequency above the static Cnmax level. The rate of in- 
crease in stall margin with reduced frequency varies somewhat from 
axrfoil to airfoil but, on the 1/rev basis, airfoils with low 
^njnax not gain enough additional lifting capability 
through unsteady effects to overcome a poor static stall perform- 
ance. This fact was not obvious at the start of the stall delav 
research. 


Compressibility affects reduce the boundaries at M = 0.6 for 
the same airfoils of figure 4. At M = 0.6 the stall delay effects 
are not as critical as they are at lower Mach numbers. Above 

detrimentaJ effects could result from operation at high 
lift in the presence of negative rates of change in angle of attack. 
Thin and transonic airfoils display larger Cnmax changes at the 
higher subsonic Mach numbers than do the thick and conventional 
sections . 


Aerodynamic Damping 

The damping characteristics of six airfoils in 1/rev pitch 
oscillation at M » 0.4 are shown in figure 6. It is immediately 
evident that the thin and the transonic sectiou .3 are positively 
damped through dynamic stall, while the conventional airfoils. 
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VOOll and V2 3010-1. 58, are negatively damped, 
dynamic damping from oscillating airfoil tests 
related to the rotor environment, it is useful 
stall recovery characteristics of the airfoils 


Although the aero- 
cannot be directly 
as a measure of the 
employed. 


The aerodynamic damping is sensitive to Mach number and fre- 
quency changes, as illustrated in figure 7, but rotor test ex- 
perience has shown the M » 0.4 condition and the lower frequencies 
to have dominant effects over the retreating blade. 


4 . ^ review of all the data for the airfoils of table I shows 
that positive damping at the 1/rev to 2/rev frequencies (k up to 
0.12) is always associated with gradual quasi-steady stall, while 
negative doping during oscillation is typical of sections with 

conside?2^hf f f transonic sections 
ffafi limited maximum lift range and gradual 

Sotiirt*? while at the same condition the VOOll and 

leading edge stall (abrupt stall). This cor- 
ali^^t°*^ static stall and damping appears to hold gener- 


Figure 8 shows that as the stall character of the VR-1 air- 
foil chanps with Mach number so do the 1/rev damping character- 

illustrated in figure 8, one way to quantify 
llnll experimentally the static stall hysteresis, 

fects stall IS associated with significant hysteresis ef- 


•..u data mentioned up to this point were for airfoils 

with high maximum lift potential and abrupt stall characteristics, 
a low maximiOT lift range and gradual stall. The only way 
of achiev-ing high lift with gradual stall is through the attain- 
ment of trailing edge stall characteristics* Figure 9 is the only 
experunental evidence we have showing that trailing edge stall is 
beneficial in changing the aerodynamic 1/rev damping. Of course, 
since leading edge boundary layer control did not improve the 
maxim^ lift characteristics of the V2 3 01 0-1. 58 at M » 0.4, this 
method is not recommended unless positive damping is clearly es- 
sontrax • 


stall Delay Characteristics 


The stall delay functions of three conventional airfoils are 
shown in figure 10. Figure 11 shows the stall characteristics of 

the VR-1 and NLR 7223-62, both transonic sections. A few qeneral 
observations can be made; general 
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(a) Normal force alwa^ i 'Experiences a larger stall delay 
than the pitching moments. 


(b) Thin airfoils experience a larger lift-stall delay 
than thicker sections, but not enough to overcome 
the static maximum lift deficiency associated with 
thin airfoil stall. 

(d) Stall delay effects on the thin sections remain 
significant to higher Mach mnnbers. 


The two transonic sections of figure 11 display strikingly 
similar trends. These trends are not substantially different 
from those of the conventional airfoils of figure 10, but the 
stall delay fxinctions of the transonic sections do not decrease 
uniformly or vanish with increasing Mach number. The lift stall 
delay at Mach numbers above M » 0.5 is attributable to beneficial 
transonic effects not better quantiifed at this time, but it must 
be pointed out that stall delay at Mach nxambers above M = 0.5 is 
not useful on present-day rotors. • 


CONCLUDING REM2UOCS 


There are differences between the unsteady characteristics of 
conventional and transonic airfoils, but the transonic sections 
have no unusual stall delay characteristics at a Mach number M 
of 0.4, where such effect would be beneficial in compensating for 
a low static— lift capability. The transonic airfoils experience 
significant stall delay at Mach ntimbers above 0.5, a benefit 
which cannot be taken advantage of by present rotor technology. 

The only favorable unsteady behavior we can deliberately de- 
sign for is positive 1/rev and 2/rev aerodynamic damping at Mach 
numbers near M = O.4.- This can be accomplished by prescribing 
gradual static stall characteristics . Because of the high maximvim 
lift requirements at M =« 0.4/ this means prescribing trailing edge 
stall . 

Low maximum lift characteristics in the quasi-steady environ- 
ment cannot be offset by a large dynamic stall delay. 

The highest priority in unsteady aerodynamics research is 
the detailed understanding of the dynamic stall process. Until 
such understanding is reached, additional testing of airfoils is 
secondary to the detailed analysis of the data we already have, 
because approximate stall delay functions for most sections can 
be obtained by interpolating or extrapolating existing trends. 
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Figure 3.- Dynamic stall delay characteristics. 
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Figure 9.- Effect of leading-edge boundary-layer control on stall and damping. 
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TRANSONIC FTX>W OVKR THE NACA G4A0fjfl WITH AN OSCILLATING FLAP - 

CALCULATIONS BASED ON THE EULER EQUATIONS* 

R. J. Magnun 
Conornl DynamlcH Convalr 


SUMMARY 

ExiJloratoiy calculations of transonic flows over the NACA 04A006 airfoil with a 
quarter-chord oscillating flap have been made using a program which obtains approxi- 
mate solutions to the Euler equations with an explicit, shock-capturing, finite- 
difference scheme. The calculations, essentially invlscid and for the airfoil at zero 
anglo-of- attack in a free-stream, are at Mach numbers and X’educcd frequencies which 
wore tested in ejqjeriments by Tijdeman. The oscillatory lifts from analogous cal- 
culations by various investigators generally agree with one another better than they 
agree with Tijdeman* s data. Inclusion in the calculations of an approximate model- 
ling of boundary conditions expected at slotted wind tunnel walls tends to shift some of 
the results closer to the experimental values. 


INTRODUCTION 


Experimental work on transonic unsteady flows over the NACA 64A006 airfoil 
with a sinusoidally oscillating quarter- chord flap was described by Tijdeman and 
Bergh (reference 1) more than 10 years ago. Their results have Inspired a number 
of attempts to calculate the flows over the configuration (references 2 through 9). The 
physical arrangement seems ideal for methods based upon perturbation equations 
(references 2, 3, 5, 6 and 8) — a thin airfoil with small angle-of- attack and flap 
deflection. 

By and large, the calculations made to date have assumed invlscid flow and that 
tile airfoil is immersed In an unconstrained stream. In this paper, an abstraction of 
material in reference 9, it is demonstrated that including an approximate treatment of 
wind tunnel wall effects tends to bring some of the calculated features of the flow 
closer to the experimental data from reference 10, 


♦Work performed under contract with Air Force Flight Dynamics Laboratory, 
Contract F33615-76-C-3018. 


SYMBOLS 


C Airfoil chord 

Cl^ Airfoil lift coefflcient per radian of flap deflection amplitude 

The difference between lower and upper airfoil surface pressure coefBcients 
per radian of flap deflection 

|C' 1 Amplitude of first harmonic of osclllatoiy pressure coefficient per radian 
^8 of flap deflection 

k Reduced frequency of flap oscillation based on airfoil semi-chord 

p Airfoil surface pressure 

p^ Free- stream total pressure 

p* Free-Stream ciitical pressure 

t Time 

u Perturbation of chordwise velocity from free-stream value 

Free-stream velocity 

V Perturbation of normal (to airfoil chord) velocity from free-stream value 

X Distance aft of airfoil nose 

y Distance normal to airfoil chord, positive upward 

0 Lead angle of first harmonic of oscillatory response fimetion 

to Circular frequency of sinusoidal flap oscillation 

DEFINITIONS 


Hcduced fretiucncy: 
k = toC/2r 

CIO 

Flap Motion; 

6(t) i6, sintot 

Trailing edge do\ui considered as |KisItIve deflection 
'lypical first harmonic of a res|»nsc function: 
f(l) f sin (u,'t ' 0 ) 

0 |iositive, res|ionse leads flap motion 
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COMPUTATIONAL NOTES 


The computer program used In the present work Is closely related to the program 
used in work on the NACA 64A410, reference 11. 

Equations. Scheme 

The coupled system of four unsteacty Euler equations in conservation form is 
solved numerically using a two-step, Lax-Wendroff, explicit, finite-difference scheme. 
DifiUsion was added to suppress ra^ed overshoot in the calculated ou^mt near shocks 
and was also found to be needed to control short-wavelength oscillations in parts of 
the flow which were near- sonic* 

Mesh Arrangements, Ssmehronization 

On the order of 5000 mesh nodes arranged in several distinct grid systems were 
used to cover the field around the airfoil. Fine mesh was used around the airfoil nose; 
the basic mesh around the airfoil was 0.04 chord squares, and stretched and coarser 
meshes were used to extend the coverage to outer boundaries several chords from the 
airfoil. 

Local, airfoil-oriented, coordinate systems were used to provide mesh nodes 
along the airfoil surface. Bands of fine mesh were inserted to provide detail on the 
loci of upper and lower surface shocks. In any single mesh system a uniform explicit 
time step would be applied at each node. Many time steps would be taken in the 
finest mesh for each time step in the coarsest. Exchanges of information between the 
develoi^ng solutions in contiguous grid systems were made by assignments and Inter- 
polations. 


Boundary Conditions 

To satisfy tangency boundary conditions along the airfoil surface, the flow at nodes 
on the time-averaged airfoil surface is calculated by the method cataloged as ’’Euler 
Predictor, Simple Wave Corrector” in the survey by Abbctt, reference 12. By a 
similar process the upper and lower pressures and flow directions are matched along a 
line extending aft about 0.2 chord from the airfoil trailing edge. Further aft the wake 
discontinuity is allowed to become indistinct by numerical diflhsion. 

If the flow over the airfoil in an unrestricted stream was to be calculated the flow 
was held Invariant at the field perimeter. On the examples calculated here, the u|>- 
stream and dowustream field boundaries were about 10 chords distant from the airfoil 
midchord and the lateral boundaries were placed further out by use of stretched me.'^'h. 

A flow pattern due to a doublet and a vortex (strength commensurate with airfoil mean 
lift) plus free stream was maintained on the perimeter. 
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For the cases simulating flow over the airfoil in a slotted wall tunnel, an 
empirically determined boundary condition on perturbation velocities 

u J-0.73V ±0,17v^- 0 

was specified on horizontal lines ±1«53 chords from the airfoil* This empirical 
relation is based upon steady transonic calculations done in conjunction with tests 
which measure! pressures on the slotted tunnel walls as well as on the airfoil model. 


Computational Eitpense 

The computer programs utilized are vnidtten in FORTRAN e? .tended language and 
the calculations were run on a CDC 7600 computer. A relatively .’table solution to a 
steady flow problem would be obtained in about 2400 passes through the field requiring 
580 seconds of computation. Stationary solutions to the unsteady problems typically 
would be obtained after following the flow for about 3. 6 cycles; this would require 2400 
to 3100 seconds of computing; the expense depends inversely on the reduced frequency. 
Pressure fields at (typically) 36 steps in an oscillation cycle were recorded for forfoer 
study. 

RESUL7S 

Calculations were made at three Mach numbers; however, attention will be con- 
centrated here on computations at Mach number 0. 85. The steady free- air flow over 
the NACA 64A006 was calculated at zero angle-of- attack with the quarter-chord flap 
undeflected and deflected 1.0 degree; tiie unsteady flow for a reduced frequency of 0.179 
with 1. 0 degree flap amplitude was also calculated. These cases were also calculated 
with a slotted wind tunnel wall boundary condition included. The unsteady pressures 
were Integrated to find foe airfoil forces and moments and selected unsteady quantities 
were run through harmonic analysis to determine foe m^nitudes and phases of foe 
ftindamentals. 

The magnitude and phase of the first harmonic of foe lift due to sinusoidal flap 
deflection is shown In figure 1 , together with similar results from calculations by 
Ballhaus and Goorjian, references 5 and 6, and Ehlers, reference 2. There is relative- 
ly good agreement betv'een the four Inviscld calculations; Tijdeman’s experimental 
result (reference 10) (including viscous and wall Irtterference effects) is considerably 
separated from the calculations • 

The normalized lift due to flap deflection In steady flow (Inviscld, frec-airi Is also 
shown in figure 1 along with results from a similar calculation by Traci, reference 3. 
The magnitudes of the calculated lifts are both more than double the steady lift measured 
by Tijdeman, reference 10. 
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The disftgreement between the free-alr Invlscld calculations and Tijdeman's 
experimental results also shows In the zero-lift pressure distributions, figure 2. The 
calculated pressure distribution at Mach 0. 86 Is compared with Tljdeman's experi- 
mental data. At corresponding Mach numbers, the calculated pressures along the 

middle of the airfoil are considerably lower (and veloclttes higher) than the experiment 
al values. 

At first. It was surmised that a simple Mach number shift of about 0. 02 might 
produce better agreement between calculation and experiment; viz. , calculated Mach 
0. 86 resembles experimental Mach 0. 876 in figure 2. An inspection of the distribution 
of lift in steady flow, however. Indicates a notable difference between calculation and 
experiment; the experimental lift distributions at Mach 0. 85 and 0. 875 both show 
considerably less loading on the front part of the airfoil than the calculated, figure 3, 
The large spike in the calculated loading between 0.4 and 0.7 chord is due to move- 
ments of the upper and lower shocks upon deflecting the flap. The program calculates 
the pressure rises across shocks to be very nearly the •’normal shock" values; the 
experimental recovery is considerably less because of interaction between shocks and 
boundary layers. Hence, foe calculated loading spike due to shock migration is more 
extreme than the experimental. ♦ 

Adding a simulation of foe slotted wind tunnel wall behavior to foe calculations 
brings the calculated and experimental zero- lift pressure distributions closer to one 
another, figure 4. Similarly, the distribution of load due to deflecting foe flap is 
brought Into closer, but by no means precise, agreement with experiment by adding 
wall effect, figure 6. 

The amplitudes and phase angles of foe first harmonics of the oscillatory 
loadings In unsteady flow are shown in figures 6 and 7. Except in the region of foe 
loading spike caused by shock movement, foe inclusion of wind tunnel wall effects In 
foe calculations results in reasonably satisfactory agreement with the experimental 
values. The noted agreement may, of course, be due to fortuitous combinations of 
errors; foe agreements between calculation and experiment were not as impressive 
for steady flow, figures 4 and 5. 

Inclusion of a simple practical means of weakening the pressure rise at shocks, 
l.e. , a simple correction for the effects of shock interaction with the boundary layer, 
certainly would be desirable. This correction is needed for predictions of beliavior 
of real configurations in free-air as well as for models In wind-tunnels. 


♦Additionally, the quasl-steady experimental results are derived from tests with 1 5 
degrees flap whereas the calculations are for 1. 0 degree flap. The loading due to 
shock movement Is proportional to the area under the spike; spike height and breadth 
depend on flap amplitude in this particular style of presentation. 
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Upon integrating the pressures on the airfoil (as calculated Including wall effect) 
to obtain the lift due to flap deflection, the results (square symbols In figure 1) are In 
considerably better agreement with the data than are the free-alr results. 

CONCLUDING REMARKS 

When wind tuimel wall Influence was added to the Invlscld calculation of flic 
steady, zero- lift flow over the airfoil, the airfoil pressure distribution changed In a 
manner as If the Mach number had been lowered by about 0. 01, The wall influence 
lowers the amount of load developed on the front part of the airfoil upon deflecting the 
flap In steady flow and lessens the phase lag of pressure excursions on the front part 
of the airfoil in imsteady flow. The parts of the quasi- steady and unsteady loadings 
which are due to shock- migration are overestimated by the Invlscld calculations studied 
here; a practical first-order means for describing the weakening of shocks by Intci^ 
action wltii the boundary layer seems to be needed. 

On the unsteady thln-alrfoll problem studied, calculations done with programs 
based upon perturbation equations, references 2, 5 and 6, are In relatively good 
agreement with the results from the program based on the unstea^ Euler equations; 
the perturbation programs obtain solutions at a very small fraction of the computing 
expense needed for an unsteady solution using the explicit unsteady Euler program. 
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Figure 1.- Lift due to flap deflection for an NACA 64A006 airfoil with 
quarter-chord flap, a * 0®; M » 0.85. 
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Hgure J.- /’.eru-mt pressure distributions for an NACA 64A006 airfoil. Fi 
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Figure 3.- Chordwise load distributions in steady flow. Free-air calculation 
at M * 0,85 compared with experimental values at two Mach numbers. 
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Figure 4.- Zero-lift pressure distributions at M = 0.85. Calculations in free 
,ilr and including wind-tunnel wall Influence compared with experimental value:! 
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Figure 5.- Chordwlse load dlatributions in steady flow at M = 0.85. Calcula- 
tions in free air and including wind-tunnel wall influence compared with 
experimental values. 
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Figure 6.- (.hordwisi* load distributions in unsti’ady flow at M = O.HS lor 
amplitudes of first harm.-nies. Calculations in'lree air and includlng 
wind-tunnel wall inl lueiui' compan-d with experimental v.ilues. 
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Figure 7.- Chordwlse load distributions in unsteady flow at M » 0.85 for lead 
phase angles of first harmonics. Calculations in free air and Including 
wind-tunnel wall Influence compared w.<th experimental values. 
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ALFRED GESSOW: 


The final wrap-up session of the conference Is supposed to be one where 
wo can all offer opinions as to how useful the program has been and where It 
should go. First, I should announce that the session Is being recorded so that 
your words will become deathless prose when we print up the proceedings. I 
might just take a couple of minutes to review the objectives of the conference 
as they were laid out initially. There were three: One was to provide the 
aviation conmunlty with a review by NASA and Its supporting contractors and 
grantees of all the latest Information they've been doing on airfoil computa- 
tional and experimental work. I might say that there were forty papers of 
that type, and about twenty- four Independent papers from representatives of 
government. Industry, and educational Institutions. The second was to provide 
an opportunity for Input and exchange between the various segments of the 
aircraft industry, the academic community, and other organizations. The 
third was to conduct a forum which engaged the ccm..iun1ty and all the aviation 
interests with an opportunity to provide an evaluation of NASA thrusts, goals, 
and progress. It Is the third one, I think, that 1.> the reason for us being 
here at this session. We woul< like to get your response to the program as It 
has gone, and we would like some constructive criticism as to where It should 
go In the future. In a way, this is a second gcround of this nature, the 
first being at the 1975 Industry-NASA workshop that I alluded to earlier In the 
conference. 

I'd like to Insert a few comments of my own before we begin. I think the 

turnout, the papers, and the scope of the Information that have been presented 

here has been very gratifying to see. We have come a long way, although there 

is a lot more to go. It Is also clear to me that the computational efforts 

that have come about In the normal course of events, that weren't spurred by 
the airfoil needs, really have enabled the state of the art In airfoils to be 
where It Is today. At the same time, and I think Gary Chapman pointed It out 
to me this morning, the interest that this conference has elicited Is an 
Indication that the subject of airfoils Is the one clear area where computa- 
tions have come of age and have provided a direct application in a very visible 
way. We are beginning to see the Influence of computational aerodynamics In 
many other areas of aircraft design but It Is most clear in the airfoil 
business, and I think that Is why we have made the progress and have the Inter- 
est that we see today. 

The way we are going to conduct this roundtable, and It will have to be 
conducted on a very tight schedule. Is to have a number of Invited guests 
express their views. There are two groups; one of these which Bob Bower will 
be chairing, will consist of the workshop chairmen from the conference. They 
will be representing the viewpoints of the applications people, the analysis 
people, those dealing with experimental techniques, and the like. They are 
seated along my right. Then, starting with Mr. Abbott on my left, we have a 
group of people who were invited because they represent areas of expertise, as 
well as the views of various segments of the Industry, universities, and other 
government agenc'es. I am afraid that there Is not going to be very much time 
for a lot of extra discussion. Instead, we plan that each of our speakers 
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Iiorn at tito tablo tako no inoro than fivn minutes in turn to present thoii 
views. After they are throuqh, wo will use the reniainincj time for coiiiiiients 
and questions from the floor. 

I've probably run past my time now, so let's start with Mr. Bower and 
his group. 


UOBIUT t. BOWtR: 

Thanks A1 , I think we have quite a formidable task, particularly to start 
here with Bud and keep him to five minutes. We will get right on to it now. 
The first area that will be critiqued and summarized will be the analysis and 
design area. It will be shared by both Bud and Gary, but I think Bud will 
lead the show, and later Gary will make some additional comments. Bud. 


PERCY J. BOBBITT: 

The panel on "Analysis and Design: Needs and Issues" made a thorough 

although rapid assessment of the state of the art and identified a large num- 
ber of problem areas. It was generally agreed that methods developed over the 
past half dozen years do provide good design and analys" tools for single and 
multielement systems with attached flow. There was, however, a great deal of 
concern expressed about the accuracy, or the lack of it, of our drag predic- 
tions. No matter how hard the theoreticians try they are still unable to 
consistently agree with experimental data; all sorts of reasons are offered 
why the experimental data isn't any good. Our drag prediction methods do need 
improvement but, in defense of the theoreticians, much of our experimental data 
is deficient as well. There isn't any, what boundary-layer people call, certi- 
fied data. We need to have a series of experiments that are conducted in as 
interference-free a tunnel as we can identify, and it lias to be coupled with some 
very detailed diagnostic work. It would be desirable to have at least one set 
of data that everybody (or most everybody) could put their faith in. 

One panel member expressed concern that we may expend tO' much effort in 
trying to find improved drag-prediction methods. How accurate do we need our 
2-D predictions to be? 

Methods for calculating 2-0 flows with large separated regions are now 
appearing. The panel thought that the data base for detailed flow-fteld 
measurements on airfoils with separation was not adequate for the theoreticians 
to validate their theories and improve flow-field models. There was also the 
suggestion that the separated flow theories be modified to include the effects 
of tunnel walls. 

A whole array of suggestions were made regarding remarks on tunnel -wall 
interference. The major ones are as follows: 
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- Continued research on slotted walls st1\1 a necessity 

- Streamline-wall tunnels show potential for subsonic and low 
supercritical conditions 

- Sidewall treatnents can be Improved using thebry 

- Wall measured pressures for Interference prediction should 
be pursued 

- Wall Interference calculations for unsteady flows needed 

Streamline-wall tunnels show potential; we have a crude one working and a 
better one on the way. There are others around the world pursuing this 
approach for measurements up to the choke point. The use of pressures 
measured near the wall In combination with theory for Interference assessments 
Is another promising possibility for obtaining "Interference-free" data. 

The panel felt that we need to consider more of the 3-D environment In 
the design of airfoils, such as the effects of sweep and taper. Another fac- 
tov' of this type that should be taken Into account Is how much different the 
root env1ronn»nt Is, say, from the middle of the wing or the tip. We also 
need to pay a lot more attention to off-design problems, particularly In 
optimizers. The difficulty here Is what kind of logic to use. Clearly the 
computer program Isn't going to provide the logic; people with practical 
experience must be Involved along with those with an understanding of the flow 
physics and the rathematical processes before the proper criteria and con- 
straints can be formulated. In carrying out off-design calculations using 
optimization routines boundary-! ayer and flow stability subroutines which 
determine separation ard transition will become more Important. 

More high and low Reynolds number data at transonic speeds Is needed 
particularly for the thicker supercritical sections. Hopefully the 20- by 
6C-cm cryo tunnel can contribute significantly to these needs. 

Fi»>>rlly, a number of specific Improvements to our analysis and design 
cod '- >h. +111 required were Identified. The four which seem to be 

i ..j-s j' Indeed are being worked on to varying degrees, are listed 


• -,• -layer and trail Ing-edge Interaction "patches" 

for ly )d analysis codes 

- Separation on spoilers and trail Ing-edge flaps - theory 
and experiment 

- Improved turbulence models 

- Faster codes (should be a general goal) 


BOB BOWER: 

Thank you very much Bud. Now Gary - 
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DR. GARY T. CHAPMAN: 


Okay, I am just going to add a couple of comments to that. I would like 
to second A1 Gessow's comment. This does represent an area where computation- 
al aerodynamics has really become of age. However there are still some out- 
standing issues. Bud Bobbitt has touched on some of them already. I am not 
going to use viewgraphs, I am just going to talk about several Issues. First, 
we would all like to have one good set of transonic Interference-free data. I 
can speak pretty much from the experimental side of the house and appreciate 
that as much as anybody else. Second, I think It Is Important to clean up 
some of the areas that we are still working on, particularly for attached flows 
where I think things are really In pretty good shape, and move on to more 
important areas of work. Third, Bud Bobbitt mentioned the need for high 
Reynolds number data, and I would like to make a plea for low Reynolds number 
data. It turns out there are a couple of special applications In the area of 
high altitude drone type aircraft where Reynolds numbers tend to run below 1 
million, based on chord, and the people that do cascades and compressor type 
work indicate Reynolds numbers run from 100,000 to 1 million. The shock 
boundary-layer Interaction problems encountered at these Reynolds numbers are 
going to be a lot different than In the high Reynolds number turbulence cases. 
Fourth, there are a lot of things that need to be done to Improve the speed 
of optimization. The Idea of making many many calculations as the method of 
getting to an optimum. Is not necessarily the best way. I think there are 
some other things that can be done to Improve this. We work with a 7600 and 
I never complain about time because we can do an optimization in one or two 
minutes. That's not much time, however many of you probably have a lot small- 
er machine. Hence I believe we have to do something to speed up the optimiza- 
tion process. Fifth, I think we are moving along very well In the area of 
multielement airfoils. However, we need some data to really check them out 
Particularly for example In the work from Grumman by Volpe In transonic aero- 
dynamics^ there's really no good data for slotted transonic airfoil conflgura- 
tions. I his Is one area where the analysis tools are ahead of the experiments, 
with the exception of confluent boundaries where we still need a lot of work 
particularly on how we handle mixing. The sixth area of concern Is massive 
separation. I was very Impressed with some of the progress being made. How- 
ever, here again we need experimental data, not only to validate the 
theories, but also to provide a better physical understanding of the nature 
of the flow and hence to develop simple modeling techniques. Finally we have 
the area of unsteady aerodynamics which Is really just getting off the ground 
- It seems at this point that there are more computational tools than there 
are experimental data by an extensive margin. So here again there Is a big 
need for experimental data and with that the need to understand unsteady wall 
Interference. If we think the wind tunnel correction problems are bad for our 
static aerodynamic data, they are terrible when It comes to transonic unsteady 
aerodynamic data. I would like to conclude with a point on accuracy require- 
ments. Bub Bobbitt quoted me as saying, why the Importance on drag, and I'd 
just like to go back to Abbott's comment the other night: "We don't fly 
airfoils." So accuracy required Is to be able to separate good airfoils from 
bad airfoils, so when you go to a 3-D design you have a reasonably good start- 
ing point. So I don*t think we want to chase one drag count on drag calcula^-' 
tions for airfoils. 
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BOB BOWFR: 

Thank you Gary. The next area that we will be summarizing Is I^he 
area of facilities and techniques and Dudley Hammond will be summarizing that. 
Dudley. 


ALEXANDER D. HAMMOND: 


, The first Item listed below deals with an Issue that was brought out In 
the analysis and design discussions by Bobbitt and Chapman as well as In 
yesterday's session on facilities. In light of the rapid development of air- 
foil design and analysis codes, there Is a continuous need for us to get 
proven facilities and basic data for code verification and Improvement. 

- In light of the rapid development of airfoil design and 
analysis codes there Is a continuing need to use proven 
facilities and provide basic data for code verification 
and Improvement. 

- Considerable discussion evolved around the Issue of 
determining the validity of data. 

- Use of "standard airfoils" 

- Accurate determination of test conditions 

- Determination of the two-dimensional Ity of flow 
over airfoil at all test conditions 

- Elimination of wall Interference and determina- 
tion of residual effect on data 

- Tunnel flow quality and effect on data 

- Required accuracy of measurements and measure- 
ment requirements 

Yesterday afternoon we had a considerable amount of discussion about the Issue 
of determining the validity of the data that Is obtained In facilities, ^ong 
the subjects of that discussion was the use of standard airfoils to provide a 
basic set of data and to help evaluate the data validity as was just mentioned. 
Determination of the test conditions and just how accurate do the test condi- 
tions need to be known was part of the discussion. The determination of the 
two-dimenslonallty of the flow over the airfoils at all test conditions was 
discussed, as well as elimination of wall Interference and the effects of any 
residual wall Interference on the data. Of course, the tunnel flow quality 
and Its effect on the data Is a valid concern that we need to continue to work. 
The required accuracy of measurements and the measurement requirements are 
dependent to some extent on how wide the data base Is, and for what purposes 
the data will be used. 


As has already been mentioned and Is Indicated In the added list below, we 
need a broad matrix of data, one as broad as you can reasonably accomplish for 
any one set of data as well as documentation of the test conditions under which 
that data was run, and som assessment of the effects of any known deficiencies 
In that data. 
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- As broad a matrix of data as can be reasonably accomplished for 
any one set of data is needed as well as documented test condi- 
tions and/or limitations. 

- Development of adequate sidewall treatment for existing 2-D 
facilities needs to be expedited. 

- Requirement for increase in research on adaptive walls to 
shorten the time before they become operationally usable. 

- Assessment of tunnel flow quality by specifically designed 
programs for adequate documentation. 

- Continue effort on nonintrusive flow-field measuring tech- 
niques - LDV and holography techniques and apply to 2-D 
facilities. 

And of course, we need to develop an adequate sidewall treatment for existing 
2-D facilities and this needs to be expedited. We are working the problem; 
however, it is a requirement that is going to develop further as we go along, 
and it certainly needs to be expedited. We need to have an increase of 
research on adaptive walls, to shorten the time before they become operational 
and useable. We have a research program going on at Langley, and others are 
doing research on various adaptive wall techniques. It all seems to move at a 
very slow rate, when one is anxiously waiting for the results of the tests 
before it can be put into operational use. The assessment of flow quality by 
specifically designed programs for adequate documentation refers back to the 
point that I mentioned a moment ago. We need to have a continued effort on 
nonintrusive flow- field measuring techniques; that is, LDV and holographic 
techniques that apply to 2-D facilities, and we need to work that problem more 
vigorously probably than we have in the past. It has been shown on other 
papers, for example, that intrusive measuring devices can give problems with 
regard to accurately determining flow- field characteristics. But that's not 
to say that some of that measuring does not need to be made even though it’s 
instrusive. Even the intrusive measurements ought to be able to determine 
some aspects of the fluid mechanics flow in these highly separated regions to 
the degree of accuracy necessary for some applications. 

This review of the facilities and techniques and the needs and issues 
involved has revealed challenging areas of research requirements. It has 
pointed out that although there are ongoing programs in all of these areas, the 
pace is slow and additional effort needs to be applied in nearly all of the 
areas covered. 


BOB BOWER: 

Thank you Dudley. Next area - Low and Medium Speed Applications - 
Roger Winblade. 
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ROGER WINBLADE: 


We Just very recently completed this session and In fact* most of the 
areas that we discussed have already been mentioned. The one area that 
received a great deal of discussion and In which I think we reached a general 
consensus Is the extremely low Reynolds number field. This Is an area In 
which we are not working but should be. Potential applications for this work 
Include wind turbines, flight at very high altitudes, sailplanes and gliders. 

It was also suggested that more work Is needed on 3-d1mens1onal effects and 
the ability to predict transition. We then considered the whole range of 
real airplane effects. Airplanes have deicing boots on them, at least those 
airplanes that use the low and medium speed airfoils. The perturbations far 
exceed any of the perturbations generated In optimizing. A question of some 
concern was raised In viewing slides over the last couple of days. The start 
point and end point In optimizing the section were very close together, or 
appeared to be. Is that In fact within the manufacturing tolerance to produce 
In the light aviation arena, or In terms of variations during life time? 

In discussing near term problem areas which we*re not addressing. It was 
generally concluded that the activities underway are the ones that we have to 
live with. There's nothing that we could Initiate now that would have any 
Impact in the next two or three years. That Is consistent with what we have 
seen In the past which Is a 6-7 year time lapse. Long term activities that we 
are to keep In mind are the potential for new configurations such as canard con- 
cepts. It was also suggested that we be more adventuresome In our thinking as 
we were a^nlshed and to not be constrained to follow the traditional line, 
but to keep the options open for different, but maybe not totally acceptable 
projects. A very interesting remark was made with respect to optimization 
techniques. The analogy was drawn with structural optimization and related 
to the very short time that Gary mentioned that It now takes to optimize a 
section. When you expand that to Include more and more variables and eventually 
a full airplane, the time becomes excessive. When you start coupling that with 
the structural design systems, the total computer time becomes essentially 
prohibitive. There Is then a high potential for either losing or backing away 
from the capability and going back to empirical testing just because of the 
complexity of the problems. The resulting suggestion was that we do as was 
previously recommended, concentrate on more efficient optimization techniques, 
rather than going through the extensive calculations. 


BOB BOWER: 

Thank you Roger. We had a classified session on supercritical appli- 
cations. Within the constraints of this unclassified discussion, Ted Ayers 
from the Flight Research Center will discuss the Issues In that area. 
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THEODORE 6. AYERS: 

although the Supercritical Applications session yesterday was class- 

MV s? « 1 kfow we <H<f "ot discuss 

any Classified material In the workshop session this mornlno. I think that 

thfVSl/ people here Indicates that^almost everyone In 

thlf ^^^st started to talk, I thought 

! °rr morning. There 'vere a couple of ^ 

areas that have not been mentioned and I will bring those up. One of the 

things brought out was the need for the people working the codes to have better 
the user and to Include some of the pitfalls that are In 
ilon« it 900 d things they do. Show some oJ thrbad cases 

tations of the programs are. Also, there was a plea to simplify the codes 
Almost everyone cited a need for an experimental data base to validate the com- 
chifii^k® Includes wind tunnel as well as flight data that NASA 

providing; data for high Reynolds number, thick boundary layers, 

lulll effects, and turbulence modeling. i 

rJrfi 4 +^^ S was a consensus of support for the National Aerodynamic 
haJi III® ^ the specific concerns for the experimental data 

base was the c^^^ for supercritical airfoils. There Is Insufficient data In 

‘^‘max two-dimensional airfoils to allow people to apply these 

in the three-dimensional sense. Another Is the lack of data for 
spoilers which has already been mentioned. This was not only the analytical 
experimental area and pertains to the Inability to predict 
effectiveness because of the unsteady nature of the flow for deflected 

around quite a bit to the three-dimensional 
talking about supercritical applications. It Is rather 
difficult to divorce supercritical applications from the 3-D case where peoDle 
have to bul d an airplane with real wings. There were a number of a^as 

i^^® workshop which make It difficult to take 2-D data and tran- 
sition It into a 3-p configuration. For example, an extension of the two- 

?ItiSris®ieede2 fnl coefficients and lower thickness 

«nk+f.«'f ”®®^®‘^ for those people who deal in the low aspect ratio, swept wing 

fMnk®!^®!®®* configurations Indicative of supercrufse type airplanes I 

S Sthe?'5eS°pir' ^'’® f’®''® be^cussed 


BOB BOWER: 

Th.nks Ted. The next eree Is rotorcraft applications - Bill Ballhaus. 
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WILLIAM F. BALLHAUS: 


The statement that we have heard over the past 2 1/2 days that "we don't 
fly airfoils" Is certainly an understatement for helicopter rotors, A heli- 
copter rotor operates In a highly three-dimensional, unsteady, viscous, and 
dynamic environment. Hence, new rotor designs based on two-dimensional, 
steady analyses will not necessarily result In better helicopter performance. 

The design problem Is so complicated that no one seems to be able to define 
specifically how aerodynamic and/or aeroelastic phenomena Influence rotor 
performance. These pehnomena are usually described In general terms, such as 
"adverse transonic effects," or In terms of the Integrated effect they produce,^ 
such as "excessive pitch link loads." 

What Is needed Is a coorolnated computational -experimental program to 
develop a clear understanding cf the phenomena that limit rotor performance. 
Designers could then focus on these phenomena In efforts to extend present 
performance boundaries. In such Investigations it is Important to consider 
different operating environments corresponding to different helicopter missions:. 
Body Interference and the coupling between blade configuration and Inflow 
condition. Including wake Interaction, should be addressed. Emphasis should 
be on advanced configurations designed to meet present and future requirements. 

In the computational area, the technology Is now available to begin 
detailed qualitative studies of the Interacting effects of nonlinearity, 
three-dimensionality, viscosity, and dynamics. F. X, Caradonna of the U.S. 

Army AVRADCOM and researchers at 0NEFJ\ are Investigating nonlinear, three- 
dimensional, unsteady rotor flows. Attempts are also being made to Include 
viscous effects. These studies have provided Information concerning the motion 
of shock waves on the surface of rotors near the tip and how these shock wave 
motions Influence rotor load variations. Ultimately, quantltat'.'^ely good 
unsteady force and especially pitching moment predictions will uz required 
In the design of high performance rotors. Caradonna and his associates have 
already demonstrated the capability to predict advancing blade pressures that 
are In good agreement with experiment for simple nonlifting rotors. 

In the experimental *"ea, there Is a concern over uncertainties In test 
data resulting from a fall -re to sort out facility effects. For example, 
there have been a large number of unsteady airfoil tests, but there appear to 
be few data points for parameters that overlap. A workshop participant 
suggests that the facilities be calibrated for the NACA 0012 airfoil In the 
following way: first, document the static stall point for M<„ = 0.3 and 0.5 
(this should be done before any unsteady testing). Then, for calibrating 
the unsteady tests, data should be correlated for reduced frequencies (based 
on the half chord) of 0.1 and 0,2 for the following cases: (1) M«> = 0,3, 

a = 10° + 10° sin t and o = 15° + 5° sin t (two different approaches to 
the same maximum angle), (2) Me. = 0.5, a = 10° + 5° sin t. The same standard 
conditions should also be used to calibrate facilities for other airfoils. 
Similar types of standards should also be established for low lift, high Mach 
number tests. 


it 
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In the next few years there will be opportunities to construct and test 
advanced technology rotor configurations. For example, Hughes Is building five 
rotors to be tested, under a contract with Ames Research Center, when the 
modification to the 40' x 80' wind tunnel Is completed. Some of these conflgu- 
rat ons remain to be defined. It was suggested In the workshop session that: 
(1) the multi-point designs proposed by Blackwell, Bingham, and Hicks, all of 
Which are based on two-dimensional, steady analyses, be considered as candi- 
dates for evaluation In this contracted effort, (2) In such design and design 
evaluation efforts, parameters be limited to focus on a few aspects of the 
rotor performance, and (3) a cojrmittee be established to coordinate NASA 
(Langley and Ames) and Army activity In this Area. In the future, more sophls- 
icated designs. Including planform and spanwise section variations, should be 
attempted. These designs would be based on three-dimensional, unsteady, 
viscous, and dynamic analyses, taking advantage of advanced computational 
methods as they become available. 


Finally, I would like to Include In this discussion a wish list ("what I 
would like to see accomplished In the next year") compiled by one of the work- 
shop participants: (1) experiments with good flow visualization (especially 
of transition and vortex shedding phenomena), (2) static stall hysteresis loops 
reoorted In all airfoil stall experiments, (3) experimental data on Mach tuck, 
(4) Investigations of the sensitivity of airfoil designs and Its effect on 
Made tolerance requirements, (5) Information on the extent to which laminar 
flow exists 1m flight, (6) results of wind tunnel tests on airfoil designs 
proposed by Btackwell, Dadone, and Hicks, (7) completion of a report on the 
AH-lG program. Including power requirements and pitch link loads, and (8) better 
communication of technology among the NASA centers. 


BOB BOWER; 


Thanks Bill. The final area Is perhaps a catchall 
Joe Stickle will summarize the Issues and needs there. 


- other applications. 


JOSEPH W. STICKLE; 

Being In the "other" category this committee covered such things as sall- 
ing ships, hang gliders, RPV's, and things the other groups would like to but 
probably didn t talk about. The low Reynolds number requirement for data 
consistently came out throughout the conversation. One of the unique missions 
d.scussed for an RPV that I had not heard before was that of using an RPV as a 
Naval ship decoy; that Is you send an RPV away from the ship at high speed and 
then loiter or cruise at a speed matching the ship cruise speed. The RPV acts 
as a radar reflector or decoy for the ship. It has to fly along about 40 
knots. So there we have the requirement for low Reynolds number data that has 
to range from low speed all the way to high subsonic speed. It also means 
that It probably will be variable geometry. So there Is a need for low 
Reynolds number data covering a wide range of parameters. In the way of wind 
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turbines, one of the requirementi was looking at a wider range of useable angie 
of attack, and looking at such things as flat top lift curves. In the com- 
pressor, a need for cascade testing was brought out. ^nother area that was 
discussed is the radial flow effects on the boundary layer. The thickening 
of the boundary layer due to radial flow and its effect on the lift and drag 
for turbines and compressors need research. Multielement airfoil require- 
ments have already been described. The fourth area we talked about was 
unsteady aerodynamics. The effects of the propeller wake on the oerformance 
0 -^ the wing airfoil sections needs to be understood. We go to a lot of trouble 
to design the airfoil in low turbulence pressure tunnels, fix the transition, 
and stick it behind the propeller and expect it to perform like the airfoil we 
designed it to be. I think we need '•ome testing in the facilities to jive us 
airfoils that will work better in the propeller wakes and they ma:- or may not 
be the same airfoils that came up in steady wakes. Compressors have to be 
concerned with stall tolerance. Bill just mentioned about the rotorc'aft air- 
foil being a very dynamic environment or critical environment to operate in, 
but people from Pratt and Whitney told us that the 4th compressor in a turbo 
engine has probably as harsh an environment to operate in as anything you can 
imagine. The committee ended up with another "other” category and that is a 
requirement not for the airfoil developers but rather for us in flight area 
or design mission area to do. If we are going to look at specific aircraft 
like an agricultural airplane or an RPV, we need to provide better mission 
criteria for the airplane from which then we can decide what this airfoil 
should be doing in terms that the designer can use. 


BOB BOWER: 

Thank you Joe. Al, I think we finished on schedule and I will turn the 
session over to you now. 


AL GESSOW: 

We sure did! I hope my team can do as well. I think we will start with 
"S" Abbott on my left. We are very happy to have him here and gain his 
insight and v'ewpoints on the work we have been doing. 


IRA H. ABBOn: 

Well, let me say first I didn't realize until a few minutes ago that I 
was supposed to get up here to say something. I had presumed that other people 
knew as well as I did that I was no longer technically competent to deal out 
any words cf wisdom about this research that you are doing. 

I would like to say that I am favorably impressed by it. A lot of it I 
just didn't understand - not in detail anyway, i would like to emphasize that 
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Impressed by the whole thing, because 1 may end up making 
comments that may give you the Impression that I am not. and that Is not so, 

Zt ^ to leave that Impression with anybody. I am very much 

it should be continued, and I don't want any misunderstanding on that, 

^ “'’easy feeling 

ba®'"' said many times that we never fly 

airfoil sections. We have probably beaten that to death, but It Is something 
that can stand being beaten to death. Even In the old days 40 years ago this 

beaten Into people all the time, and I think that wnen I 
sections I must have spent a third of my time 
;!lil Yes. we're doing a lot of work on the airfoil sections, but 

kid yourself, they are not wings. The miracle, of course. Is that alr- 
^5Hf slg^flcance whatsoever. The concept dates back 

1912, and that Is a long while ago. The assumptions 
that wei^ made in the theory which led to the concept of airfoil sections are 
TJ.V^ reality that it Is a wonder that there Is any reason why any^dj 
should consider airfoil sections Important. The wonder is not that the airfoil 

practical results on three-dimensional forms 
validity that Is worthwhile. Of course. It does. You all 
know that too. We also knew 30 years ago that even If you completely disre- 
gard the viscous and compressibility effects, which become overwhelning under 
many circumstances, that on planforms that we are commonly using today the 
curvature of the flow is such as to Induce In an airfoil In many sections of 
"®9ative camber which may greatly exceed the total amount of camber 
put In the wing section as designed for two-dimensional flow. So, all I am 

® ^ast as Is feasible with three-dimensional cases. 

I know that everybody wants to. I know that the work here nas to be done first. 
In all probability we have to learn to walk before we can run. 

enamoured by It. I have an uneasy feeling that at 
least some of the work is being done not because of perceived needs by the 
aircraft designers, but because a lot of people playing with computers have 
discovered problems that they can do and they are having one very good 
time. Now, I am not against research people having a good time. In fact I 
am In favor of it, but don't forget that while you're having a good time that 
ther^ are more important problems to solve. So, instead of thinking up little 

9 lve you a little better result on 
pis, that and the other thing, do a little thinking about how you are going 
to up thep fine computers to solve some real problems. I don't know how you 

do that, pu know more about It than I do. I just ask you to think about It 

once in awhile, when you are no longer enjoying yourself quite so much. 

I hpe also noticed another thing - we're gradually changing our whole 
design phipsophy, not only in aircraft but on other things - to optimization 
of thp, pat and the other thing instead of just selecting things and putting 

them together and coming up with a pretty good thing. This Is good. It is 

nppsary. But again don't forget that when you start optimizing things, 
mii^h depends on what you optimize and how you do It, If you put too 

much faith In the optimization of things which are after all not going to work 
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the way you think they are going to, you may end up with something that over- 
all Is very very far from the optimum. Above all look at your optimum; if it 
Is a nice smooth curve that Is concave downward you can pick a maximum point 
without trouble. But If It has a sharp peak, forget It. Fair It off some- 
place otherwise you have nothing. Now, I have faith In designers, having 
talked to a great many of them over the years, so I think they are not going 
to be bothered too much by this, I am talking more to the research people now 
than I am to the designers. 

One other thing I am Impelled to say. I have heard a great deal about 
computer work, I have heard a great deal about wind tunnel work - you know I 
'am an old wind tunnel man - I heard a little about flight work. But I want to 
make a plug for It. There Is a great deal of this work that I think can only 
be properly done In flight. The question of transition Is something that 
there Is no other way to solve, or get basic Information, except In flight. 
Now, I thl^ I know something about low turbulence wind tunnels, and If you 
want a low Reynolds number low turbulence wind tunnel, I think you can get one 

that will be pretty good. But If you want to get Into what we think of as 

real airplane Reynolds and Mach numbers at altitude - I will leave some of the 

little airplanes out of It for the moment - I think It Is going to have to 

be done In flight. I am afraid there Is some thought of flight research as 
flying an airplane and see what happens to that airplanet or getting some data 
on It. This Is not the kind of flight research I mean. We are going to have 
to go out and modify airplanes to get the detailed data that would be obtained 
In wind tunnels If suitable tunnels were available. In the real world of 
flight you cannot only find out how or whether some definite aerodynamic 
promise can be achieved, but also how It has to be done so that when things 
go wrong you are not going to have penalties that are going to be too severe. 
So I would like to make a real plug for flight research. 

I have said too much already. I thank you. 


AL 6ESS0W; 

Thank you very much “S". We certainly weren't disappointed In expecting 
and getting very meaningful comments. They were just the kind we are looking 
for. I think we will be hearing more about the extension from 2-D to 3-D as 
we move down the table. Our next speaker will to Helmut Sobieczky who is 
from Goettingen and Is visiting at the University of Arizona. 


DR. HELMUT SOBIECZKY: 

The success of a conference like the one we have had the privilege of 
attending here depends critically on the amount of stimulation It provides 
through reports on major achievements and through work that we see remains to 
be done. This goal was achieved admirably here. The clear demonstration of 
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f”-- »'>-f“"5. “"‘lor- 

iire needed in applyina airfoil result? tn”th^^ clearly, where advances 
the whole aircrSn; To Mhieve ?hil "'"9= 

both the theoretical and tte“e*peH«“a? r^s’Sl??‘'c?“Jully?""‘' 

conpe^S?s'’f«“^alSf^i^;“l?,SeS^?„'^Kr'^ “"? by 

estihlishing realistic LIu of ?el? ‘ *,‘'’0 

explanations for observed phenomena alth^thL?! 5^ 9(*tng physically sound 
their value in airfoil resMroh^Sd'>i!?‘5.f?®°’^^ 2"? experiment have proven 
In three-dimensional ap?ll“?flln?? -lelineated its value and limitations 

cal rl??l?s'^? ‘bat it seems worth noting here Numeri- 

sSch eJjerliSentriJe nwded“SS®l5pJiy^or^5lnd1?at^ 

interaction problems Parole! flow 

ties based oS sJnS^itic cJi™ioJ^?^nTfesulST 1 ^‘^ 


AL 6ESS0W: 


had hlH’Jivfhis 2Se^ls‘■?^^75'■‘'?11S^ o?‘SS?"m*"3 ' ™S‘ 

the 3-D L>. 1 th?:rSe’;,m:'h.l nrcl^?o‘rtrete"mSf?s" 5? tt^‘ ?ul5:??. 


LUIS R. , .,<AsNDA: 

by saying that ! am*"very^impresLd°with^t^^ Abbott did, 

e«o"u“ge3?‘'SSt°a"ft;? llZo Lt'? *"12 «frkTh*oildTe”S2"ti2"ufd afd" 

a"?I"lw2*typ« ?^hS?e'ihSlnil^"2!;?lEfa^i‘2S'‘°'*""^ lhallhere 

p*s??Lt*"t-*sim:j??irci2^.isrI>iii* S 
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we are In terms of 3-D wing design technology and what 

?IseS?ch irrequlred to bring this technology to the same I®*®’ 

and completeness of 2-D airfoil design J« Ja« seen 

r4l"eJ?:nrSlrfili "rS«arKlf"c5nfS?e r 

tte proper place to become the focal and lead ng P®'"* JfX 
S? r?hSold acknowledge that, from the airplane ™"V!!®‘"5"4tS?lrin 
9 n airfoil research may be reaching the point of diminishing 
oj- a?wal deHortShnology. I am folly aware that my opinions are biased 

point of view of what large segment of the Industry which worKS in mgmy 
three-dimensional environments. Thank you. 


AL GESSOW* 

experienced In experimental techniques as well as other areas. 


LARS H. OHMAN: 

I'm very glad to see that NASA Is making a quantum jump 

First 5^hlve°thrs1 dewall boundary -layer problem, which has 
"urSo'real suggestion ha^been .ade on how o^^^^ r^ood^rpr^ c^aJd^o^r 
system of suction o^'er an area around the airfoil is a^goo^^ J^methlng could 
experimental data seems bo support this. Howev . boundary -layer 

be done here in the form of analysis. I m that the bounw^^^ 

concept Is all fine, but the overwhelmi ng pro blem Is s*<l] ?SjS?4?eMe*Sffects. 
sidewall boundary layer. Oust ? short co^nt on the of 

Wall pressure measurements should be mandatory on all d oi 
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course* so that sufficient data are at hand to employ some sort of wall inter- 
ference correction. The trend seems to be for everybody to use the 0012 air- 
foil to evaluate new concepts. I think there is a very great danger in doing 
this, because that's an airfoil that's not particularly sensitive to wall 
effects. A supercritical airfoil is much more sensitive and you may find 
yourself entering into a completely new ball game if you start evaluating the 
adaptive wall concepts with supercritical airfoils. One area of interest is 
z-D high lift at low speed and I just want to make a short comment here. A 

successfully tried In Canada is that with a 
solid floor and a ventilated ceiling. This seems to give almost an ideal wall 
configuration for high lift testing, and this was successfully demonstrated 
at the University of British Columbia. Some further exploring of this concept 
would perhaps be well worthwhile. There is also the question of data accuracy. 
As I mentioned yesterday, I feel that in one way we are working in the dark 
because we have not established a target accuracy that we should strive for. 

We want to compare new calculation procedures and wind tunnel data and produce 
Cn plots, which may look beautiful, but very little attention is paid to 
what accuracy is actually required in such a comparison. I think that we 
snould establish a target accuracy for what we want to accomplish instead of 
just carrying on ad infinitum Improving things without even knowing if it is 
needed. The same applies, as was touched upon by one of the other gentlemen 
he»^, to model accuracy. This has not really been discussed at the conference, 
out we know, and it has been proven, that the supercritical airfoil is 
extremely sensitive to model inaccuracies. That's all I wanted to say. 

Thank you. 


AL GESSOW: 

Thank you. Our next speaker is Andy Lemnios from Kaman Aircraft. 


ANDREW Z. LEMNIOS; 

During this past two-and-one-half days, we have heard a great deal about 
new, pow®rful analytical techniques for the design of airfoils and comparisons 
of these advanced analyses with pressure distributions on newly developed 
airfons. We have also heard, in the past two-and-one-half days, that you 
cannot fly an airfoil. I would like to expand that statement to say that you 
cannot fly a computer either. 

With regard to rotorcraft applications, which was the subject of one 
specialist session this morning, Dr. Ballhaus very succinctly sunmarized some 
of the concerns and some of the "wish list" items for airfoils that are pre- 
valent throughout the rotorcraft community. I support and strongly reinforce 
the conments by Dr. Ballhaus. Furthermore, I want to point out that in no 
other application do we have such strong interactions among the aerodynamics, 
dynamics. Inertial, and control characteristics of an aerodynamic structure 
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as we do In the rotor blade. The aeroelastic effects In rotary wings can be 
much more significant and much more Important than any airfoil Improvements 
which we may be seeking* 

In trying to develop advanced airfoils with Improved performance, compari- 
sons between two-dimensional airfoil tests and theory have been directed 
primarily at matching pressure coefficients and pressure distributions at 
constant or matched lift coefficients. As Dr. Ballhaus brought out earlier, 
rotor blade designers are equally Interested In pitching nmment coefficients 
and their variations with Mach number and angle of attack. In order to 
develop analytical tools that are useful to the rotorcraft aerodynamicist, 
added emphasis should be placed on matching pitching moment coefficients. In 
addition to lift coefficients. 

Another point, which Is Important from a designer standpoint and from a 
fabrication standpoint. Is that minor changes In surface contours on these 
new supercritical airfoils with tralling-edge cusps result In strong Inter- 
actions and strong variations In the airfoil pressure distribution. Typically, 
rotor blades have chords that are on the order of one to two feet. Required 
contour tolerances on blades of this size are very difficult to achieve 
practically, thereby creating a problem for production. As a recommended 
approach to this production problem, we should determine the Influence on 
pressure distribution of minor variations and tolerance differences over the 
leading-edge portion and the trail Ing-^edge cusp of the blade. What would these 
contour differences do to supercritical flow? Would they. In fact, create 
more problems than they solve? 

A final point Is that many of these new airfoil shapes for fixed wing 
applications are based on a single point design condition, such as cruise* 

In rotorcraft applications, we have a multi-point design problem. We have to 
design for a cruise condition which has a wide variation In flow over one 
cycle and we have to design for maneuvers. What happens when you go Into a 
1.5 g pull-up or a 0.5 g pushover? What does this do to the air loads? Per- 
haps we can design an efficient airfoil section and blade for a particular 
cruise condition, but It may have a hard stall. The point there is that we 
need to look at a forgiving airfoil - one with a gentle stall. This point 
was brought out yesterday In one of the papers. Thank you. 


AL GESSOW: 

Thank you very much, Andy. Our next speaker represents another type of 
rotating airfoil application, Mr. Carl Rohrbach from Hamilton Standard who 
represents propellers, wind generators, and other types of rotating machines. 
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CARL ROHRBACH: 


tho participate In the roundtable discussion of 

the material pi^sented at this outstanding conference* I would like to conment 

?s »“*■• ‘"an for wligHld JJtorTaSdttn 

that ^5® '■•*"*'■ f»l“tf««ent In this case Is 

that the airfoil sections range from two to ten percent thick over the Imoortant 

«ntal*dSn~ *nrf llitU destgraJd eij^l?'"* 

mental data are available on advanced airfoils for this thickness range. 

A1rfol?r satisfied with the NACA Series 16 

hair Hflh airfoil family to 

have high critical Nach numbers and good llft-to-drag ratios over falrlv wide 

M^ril^Ir f*"l; ®!?.”®®‘' now. With the rSJer?H^^^^^ 

k 4 looks like we may be able to achieve some Improvement In 

performance characteristics of the propeller along with Improve- 
n^^"^eslgn operation. However, our concern Is that all the work done 
to date has been on the relatively thick airfoils, l.e., 10-percent thick and 
above. We knw that NASA Langley have plans to Include the design and test 
of thin airfoils and urge that they expedite this work. 

use WhUcomb Indicated that we could 

l!f®* ?? *^^nd tunnel testing, to confidently deslon new 

we need alrfoll^datH^lir®'^* performance prediction, 

ijf J?®®® over a broad range of Mach numbers and angle of attack 

we require such data above the section critical Mach number as well 
as below. Thus, we would need to use the codes to develop these data where 
perhaps they may not yet be applicable. Finally, I am concerned about the 
large performance sensitivity the supercritical airfoils have to small con- 

propeller blades are probably within manu- 
!r ®”?®*4 i®®?»'d1ngly. we would like to see som^ test date Sn 
prevlously!*"^^^*^*^ airfoils covering the broad operating ranges I mentioned 

deslan”and"Sr?*«LiISii^ urge that NASA expedite and expand their plan to 

whetSr Jho! supercritical airfoils and establish 

h2f« ^^®* pei^ormance advantage over conventional airfoils as has 

been demonstrated for the thlckr- versions. a^rrons as nas 


AL 6ESS0W: 


Thank you. The next speaker Is Mr. 
represents the RPV conmunlty. 


Harry James of Teledyne-Ryan who 


HARRY A. JAMES; 


The trend I see for the high altitude RPV vehicles Is the need for new 
data In the extremely low Reynolds number regime. The regime between 1 and 3 
n ^ ‘think Is now pretty well documented and we are very grateful to 

ir* his help In the application of the supercritical airfoils for 

the high subsonic regime. When we got Into Mlnl-RPV's. particularly the low 
speed propeller driven RPV's for high endurance » I am starting to see Reynolds 
numbers as low as 20,000 on propeller balde sections. We have very little 
substantiation of the efficiency of the propellers that would operate In that 
regime. On the wing sections of the Mlnl-RPV that Mr. Reed of Dryden Is plan- 
ning to launch from a high balloon at about 100,000 feet, the Reynolds numbers 
can be as low as 100,000. This area was not addressed at this conference. 

little comment about the other end of the Reynolds number regime. 

TRA Is now Involved In a Mach 2 RPV for low altitude which will operate In an 
extremely high Reynolds number area. During the preliminary design test 
process development, I noted that we had very Httle confidence In the avail- 
able analytical procedures for the selection of the airfoils so we used the 
old cut and try reliable wind tunnel tests Involving many configurations to 
solve your problems. I think It's a cumbersome way. I didn't see anything on 
supersonic airfoils at this conference. Perhaps It was Irrelevant. 


AL GESSOW: 

Thank you very much. Our. next speaker Is Mr. Bruce Carmichael. 


BRUCE H. CARMICHAEL: 

First, I think we all owe Mr. Ken Plerpont and his wonderful staff a 
vote of thanks for an extremely well organized. Informative, and exciting 
symposium. 

My own Interests cover viscous fluid phenomena over the entire Reynolds 
number range. As the late Dr. Sighard Hoerner so beautifully put It— "From 
butterflies at 1000 to bombers at 100,000,000." Perhaps we should expand 
that to. From the fractional Reynolds numbers of falling dust motes to that 
of nuclear submarines at 1,000,000,000." I would like to address three 
regimes In this expanse for about one minute each. 

First the regime from 10,000 to 100,000. You may think this will mainly 
help your children as they strive to win model airplane contests and perhaps 
this Is not a bad reason for research, but I assure you that In addition 

® project of major national and International Importance. You 
win be hearing about this mind boggling project next year. Scientifically, 
the problem Is that of the long laminar separation bubble on wing upper 
surfaces. To define the separation bubble. Its separation point, Its transi- 
tion, Its turbulent reattachment point, the pressure distribution under and 
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well established. Force and pressure d1str1butlon'd»5» <lesl9n Is 
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and 60*8 by our group working under Dr. Werner Pfenntnger at Northrop Aircraft. 
Following those exciting days, I experimented with laminar flow under water 
for over 15 years. The Navy, who backed some of this work, was Interested In 
laminar torpedoes, but It applies equally well to laminar aircraft fuselages. 
We have worked both with and without distributed suction an(l also surface 
temperature differential. The research and development Is largely done. It 
Is now time to reap the benefits with useful vehicles. 


AL 6ESS0W: 

Thank you Mr. Carmichael. Last, but certainly not least. Is my good 
friend Barney McCormick, Chairman of the Aerospace Department at Penn State. 
He represents the university viewpoint In rotorcraft and a number of other 
good areas. 


DR. BARNES W. McCORMICK, JR. : 

I am pleased to represent the acadmnic community. As we have seen at 
this meeting, the academic community has certainly contributed to the area of • 
computational fluid mechanics. Particularly noteworthy is the work of Lee 
Carlson at Texas A&M, Joe Thompson at Mississippi State, the group at Wichita 
State University, and the airfoil center which has now been established at 
Ohio State. I would also like to mention some research at Penn State In an 
area which has not been touched on here. This concerns the field of hydro- 
dynamics. We have a large water tunnel at Penn State and have had an extensive 
research project ongoing for many years concerned with underwater propulsion, 
particularly torpedo propellers and pun^) jets. Through the application of 
aerospace technology we have been able to significantly Improve the cavitation 
performance of torpedo propellers. Classified at the time, we succeeded with 
the application of series 16 airfoils In doubling the silent speed of torpedos. 
I certainly think that the supercritical technology has some place In hydro- 
foil development for fully wetted sections. One Is concerned there with the 
same limitations which are encountered In compressibility effects; that Is, 
how high can you get the value of the minimum pressure coefficient. 

I might add something other than an aerodynamic consideration concerning 
the design of propellers, some things one would normally not foresee. Just 
last week I visited a propeller manufacturer In connection with a research 
project we have with NASA Langley on propellers. I asked them If they were 
considering going to supercritical sections of some kind for their propellers. 
The response was not too enthusiastic bec?Mse they were concerned about the 
fact that. In operation, a propeller becomes nicked on the edges. The edges 
look to them to be far from the neutral axis of the section for a supercritical 
airfoil and they are concerned about fatigue problems. 
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AL GESSOW: 

tako speakers here at the table. Now Bob Bower will 

floor**'^^^ utilize the remainder of the time to get some comments from the 


BOB BOWER: 

♦w ^ so obedient in staying within our 

time that we now have about 15 minutes left before we* re scheduled to complete 

^ '*® ^0 open up the discussion now^to 

hi?nrh comments that you may have. Well, a real bashful 

bunch. One down here please, David C. Ives. 


DAVID C. IVES; 

*" engine coiiipany. Pratt-Whitney. If you count the number of 
them have engines, and a lot of these engines are jet 
iSrtA”®*c yoo count the number of blades In the engines it's 500 to a 

I? ^i® ”“"*er of blades have the number of wings outnumbered by 

at least 2 orders of magnitude* We need some basic research on Reynolds num- 

"'1110". with levels of turbulence from nothing to 
5 percent. That s what we really need. The rest of the tools you people have 

U?klSg! "" 


BOB BOWER: 


Thank you. Behind you over there. Dennis Bushnell. 


DENNIS M. BUSHNELL: 

I'd like to conment on a question of the calculation of turbulence 
’"® 9 lons. I think first of all we'll have calculation methods for 
turbulent separated regions way before we have accurate calculations. The 
comment by Bradshaw in a 1978 AGARDograph is we know virtually nothing about 

turbulent separated flows. The fundamental 
^®®^^ threefold: There is breathing of the separated 
regions which superimposes on the turbulence, there is generation of 3- 
dimensional vortex structures, and there is interaction between these 3- 

nl!!!®?«^2II® yc’^'tex structures. As a final comment, the actual work going on 
now in the structure of turbulence and separated flows is not very much; in 
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fact. It's extremely small. At the present rate of progress, we may never get 
there within Al's lifetime, as he commented before. 


BOB BOWER: 

Thank you. Now Ted Ayers. 


THEODORE 6. AYERS: 

I'd like to make a couple of comments. Bob. One of the things that I 
didn't mention earlier was that people are saying that NASA ought to have 
faster dissoni nation of data. The manpower situation Is what dictates when 
we can get the data out, generally. The other thing I'd like to say Is “thank 
you" to Ira Abbott for saying what I would like to have said about flight. 


BOB BOWER: 

I thought you'd like those words Ted. Are there any other comments? 
Before I officially close the conference I think I would like to say a few 
words. We've all been looking forward to this conference for a long time - 
In fact some people for ten years • we've been looking forward to the time 
when we would be able to have the data as well as our theoretical tools to 
talk about and to exchange. So It has been a very gratifying experience to 
see It come off with such attendance. A conference Is not an end In Itself. 

As had already been mentioned, one of the primary purposes of a conference Is 
to stimulate thinking, stimulate exchange of Ideas, and I think with an 
attendance here of 450 people we have done that. I want to thank all of the 
parti cl pr. its, not only In this session this morning, but also for the past 
three days. And again, I'd like to ask the Individual who Is really the one 
most responsible for the conference, Ken Plerpont, to stand up. He has 
finally come In here and joined us. As Ken Indicated to me at coffee break, 
he organized this conference to be one that he would like to go to. I think 
that he has been extrenmly successful In putting on one tr^ndous conference. 
With that, I would like to officially close our conference on airfoil research 
and thank you again for coming. 
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